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The data for this report were obtained as a result 
of investigations carried out under the Environmental- 
Social Program, Northern Pipelines, of the Task Force 
on Northern Oil Development, Government of Canada. 
While the studies and investigations were initiated to 
provide information necessary for the assessment of 
pipeline proposals, the knowledge gained is equally 
useful in planning and assessing highways and other 
development projects. 


Résumé en francais 


Le présent résumé indique le degré d'avancement des études 
hydrologiques effectuées en 1973 par la Division de la 
glaciologie de la Direction des ressources en eau, dans le cadre 
du Programme socio-écologique des pipelines du Nord. Certains 
des projets mentionnés ont débuté en 1973, tandis que d'autres 
représentent la poursuite et l'extension du programme de 1972, 
Seize rapports de recherche contenant les données et les résultats 
d'étude accumulés jusau'a la fin de 1972 ont été publiés en un 
volume dont on peut se procurer un exemplaire 4 Information 
Canada. 

En 1973, nous avons surtout visé a établir les données 
hydrologigues gui n'ont pas fait l'objet de nos études de 1972. 
L'objectif secondaire a 6té de poursuivre les travaux grdce 
auxquels de nouvelles données pourraient contribuer a améliorer 
considérablement la base statistique de certaines conclusions. 


Au cours de 1973, nous avons poursuivi les travaux entrepris 
dans les bassins de recherche du ruisseau Twisty, du ruisseau 
Boot et du lac Peter, ainsi que les études de la débacle, de 
Pembacle et de l*affouillement par la glace dans le chenal 
principal du Mackenzie. En collaboration avec le Groupe d‘étude 
Ecologique de la route du Mackenzie, nous avons recueilli des 
Gonnees sur Ve debit de 5 petits afiluents du Mackenzie, afin de 
verifier I“exactitude dés critéres de calcul des buses de la route. 
Dans le delta du Mackenzie, nous nous sommes attachés 4a 
déterminer les taux de sédimentation au moyen de techniques 
dendrochronologigques, les taux d'érosion des berges des chenaux et 
la répartition saisonniére du débit entre les chenaux du delta. A 
l'ouest du delta, nous avons entrepris dans les monts Richardson, 
l'examen des caractéristiques physiques des bancs de neige semi- 
permanents. 


En plus des recherches in situ, nous avons continué 
1'étude des caractéristiques géomorphologiques et hydrologiques 
des bassins versants des affluents du Mackenzie. Parmi les 
travaux en cours, le rapport mentionne pour la premiére fois, 
les expériences en laboratoire relatives 4 la pollution par le 
pétrole des riviéres englacées. Au cours de l'année, nous 
avons terminé une nouvelle étude globale concernant la 
Simulation des effets résultant de la modification du mollisol 
dans les régions de pergélisol du Nord canadien. 
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SUMMARY REPORT ON 1973 HYDROLOGIC STUDIES 
RELEVANT TO NORTHERN PIPELINES 


BY 


D. K. MACKAY 


Introduction 

This report summarizes the progress of hydrology studies carried 
out in 1973 by Glaciology Division, Water Resources Branch, within the 
Environmental-Social Program, Northern Pipelines. Some of the research 
projects discussed began in 1973, while others reflect the continuation 
and extension of the 1972 program. Sixteen research reports containing data 
and results of studies to the end of 1972 have been published in one volume 
available through Information Canada (9). 

Our principal objective in 1973 was to provide information on hydro- 
logical concerns not actively studied or covered by our 1972 studies. A 
secondary objective was to continue those projects in which the statistical 
foundation of some conclusions would be markedly improved by additions to the 
data base. 

During 1973 projects in Twisty Creek, Boot Creek and Peter Lake 
research basins were continued as were break-up, ice jamming and ice scour 
studies on the Mackenzie main stem. In cooperation with the Mackenzie Highway 
Environmental Working Group, discharge data were collected on 5 small Mackenzie 
tributaries to check on the adequacy of highway culvert design criteria. In 
the Mackenzie Delta, emphasis was placed on determination of sedimentation 


rates using dendrochronology techniques, on erosion rates of channel banks, 


and on the seasonal flow distribution among delta channels. To the west 
of the Delta in the Richardson Mountains, an examination of the physical 
characteristics of semi-permanent snow banks was undertaken. 

In addition to field research, analyses of both geomorphic and 
hydrologic characteristics of Mackenzie tributary basins were extended. 
Another continuing project, involving laboratory experiments into oil 
pollution of ice covered rivers, is reported herein for the first eine. 

A broad, new study concerning Simulationof active layer modification effects 


in permafrost regions of Northern Canada was completed during the year. 


Work Carried Out 


ihe Mackenzie River and Tributaries 
1.1 Break-up, ice jams and river scour. 

Reconnaissance of Mackenzie break-up continued with special attention 
being paid to logging ice jam locations and to the fluviogeomorphic effects of 
ice action (8). The 1973 break-up was relatively placid with numerous minor 
or partial ice jams occurring in the main stem down to Arctic Red River. A 
fairly large jam below the Arctic Red River junction was observed blocking the 
channel for some time. 

Inclement weather over some reaches of river during break-up affected 
the river scour project. Sections of interest to the river scour study were not 
photographed at the time appropriate for comparison purposes, thereby limiting 


the success of the project. Our intention is to continue the study into the 


1974 break-up period. 


ee Olen ce; covercd rTivers. 

Studies of 011 pollution in ice covered rivers continued with 
present work concentrated on a literature search and on laboratory experiments 
(7). Hot crude 011 released in a 1.5 m. diameter ice-covered basin was monitored 
to determinespreading rate and related behavioural characteristics of the oil 
over varying periods of time. 
ioe, Discharee and fish passage through culverts. 

A project was carried out in co-operation with the Mackenzie Highway 
Environmental Working Group to assess design criteria of 7-day delay discharges 
at culvert sites (12). Surveys of 5 small basins with culvert sites on the high- 
way route south of Norman Wells, indicated that icings on stream beds and flood 
plains are regular features of their hydrologic regimes. 

1.4 Morphometric analyses 

Work involving morphometric data collection in tributary basins, 

especially on the eastern side of the Mackenzie River, was continued (13). 

Flood characteristics of a number of basins were estimated using the relation- 

ships established in an earlier report (9). Among the larger basins investigated 

on the eastern side, 100-year floods for the Blackwater River at the Mackenzie 
(drainage area 3157 (Sq. miles) as well as the rivers Ochre (drainage area 445 sq.miles) 
and Donnelly (drainage area 614 sq. miles) were estimated at 10.37, 12.24 and 14.82 
cfs/sq. mile respectively. Among basins on the western side, 100-year floods for 

the Ramparts River at the Mackenzie (drainage area 2908 sq. miles) and the Hume 
(drainage area 1754 sq. miles) and the Little Bear rivers at the Mackenzie 

(drainage area 877 sq. miles) were estimated at 12. J iseie 44 and 11 93ecrs7sa. 


mile respectively. 


ma Her 


1.5 Research basins 


Data collection continued on water budget elements in 3 small 
research basins. For the first time spring snowmelt runoff was measured 
in Boot Creek; and in Peter Lake basin, the amount of water withdrawn for 
use in drilling operations was monitored (1). In addition, an experiment 
was commenced at Peter Lake to study the undercatch of precipitation which 


is attributable to rain gauge design. In the subarctic alpine basin of 
Twisty Creek, work on water budget components was expanded and new rainfall/ 


runoff peaks were recorded (6). 


va Mackenzie Delta 
2.1 Seasonal flow distribution 
Miscellaneous discharge measurements were made over the 1973 
summer season but not in sufficient quantity to describe the delta flow 
distribution. Beginning in the summer of 1974 the Water Survey of Canada 
will install gauges in a number of delta channels to monitor sediment and 
discharge on a continuous basis in conjunction with Beaufort Sea studies. 
Mackenzie Delta discharge and ice thicknesses were measured in 
March and miscellaneous hydrometric measurements were taken during the open 
water season (2). March 1974 flow distribution in the outer delta differed 


from that of March 1972 and 1973. The delta ice cover was generally thinner 
in 1973 than in the other 2 years. 


2.2 Channel bank erosion 


Channel banks in the south end of the Mackenzie Delta were surveyed 
to estimate erosion rates of a variety of bank types (10). The significant 


factors involved in bank recession were analyzed and discussed for this portion 


of the delta. 


The results of the study suggest it should be expanded over the 
delta to improve the data base and to determine its utility in other areas. 
so ocdimentation 

Sampling of trees in the Mackenzie Delta continued with a view to 

establishing sedimentation rates by analyses of adventitious root growth (4). 
In addition, an enquiry into 3 pingo-like humps on the Peel River flood plain 
beyan with flood level surveys, determination of active layer depths and 
establishment of vegetational zoning. A detailed study of the features including 
microtopography, active layer thickness, ground ice content and plant communities 
iSs-contemplated an 1974. 
Ce Hydroclimatology-permafrost relations 
3.1 Uses of dendrochronology 

The uses of dendrochronology in permafrost environments are discussed 
by Henoch (5). Some of the clues to past climate rendered by dendrochronologic 
techniques are important in the analyses of geomorphic processes and forms. 
3.2 Active layer simulation 

The relations between active layer thickness and regional climate were 
assessed using computer simulation techniques (3). The effects of surface melting 
and drying as well as destruction or removal of the organic surface mat on active 
layer thickness were also investigated by computer simulation. The results are 
presented in the form of isoline maps showing the geographical distribution of 
active layer thickness and the degree to which it is changed by surface modification. 
3.3  Semi-permanent snowbeds 

Snowbeds in the northern Richardson Mountains and adjacent foothills 


were examined to associate their distribution and physical characteristics with 


localized slope stability and biological activity (11). These features, 
which may exist for more than one year, provide a moisture supply in summer 
affecting active layer development, erosion and mass wasting of slopes. 
4. Projects not reported. 

Other projects which are still in progress such as lake temperature 
regimes, Mackenzie river heat budget and some snow -cover/ground temperature 


relations will be reported on at a later date. 


Significant Results 


The general conclusion is that the most favourable Mackenzie 
pipeline route from a hydrological point of view remains as in 1972 
CErgure sl): 
1 Ice Jams - Partial ice jamming in the main Mackenzie may be as significant 
from a river scour point of view as complete cross-sectional jamming. When a 
partial ice jam develops, such as occurred directly below Fort Simpson in May 
1973, the channelling of flow can cause considerable scour (8). This has been 
indicated by sediment-laden plumes observed streaming below ice-constricted 
channels. 
i Mackenzie Delta Flow - Three year's of March discharge measurements show 
that more than 80 per cent of total Delta inflow is carried by Middle Channel 
north into the central part of the Delta (Station 18). Farther downstream the 
pattern of flow dispersion is unstable from year to year as is indicated by the 
range in discharge from 57 to 91 per cent in Middle Channel (Station 20) opposite 
Reindeer Station (2). This indicates that water Supply in outer distributary 
channels could be a critical factor in winter oil drilling operations. 


ae Ice Covered Rivers and Oil - In ice covered rivers the release of crude oil 
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Figure 1. Mackenzie pipeline routes from a hydrological point of view. 
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can cause severe environmental problems. Our laboratory experiments indicate 
that hot crude oil released under an ice cover separates into small globules or 
particles 0.1 to 2 cm in diameter. These particles form an 011 layer about 

1 cm thick under smooth fresh water ice. The particles do not adhere to the 
ice under-surface but become sandwiched if the ice is thickening. A period 

of 1 to 3 days is needed for the oil to become sandwiched; variations are 
dependent on ice growth rates and current velocities (7). 

4. Delta Bank Recession - An examination of bank erosion rates in the 
Mackenzie Delta shows that these rates may vary from negligible ones to 30 
metres per year or more depending on shape (10). Before structures or other 
facilities are located adjacent to banks, an examination of bank shape and 
erosion rate should be made and Outhet's study (10) should be taken into account. 
Di Design Floods - New flood peaks which should be considered in establishing 
regional design criteria for stream crossings by pipeline and highways were 
measured in our research basins. In Boot Creek, draining a 31] ee area of the 
Panipou Hihlse sa vysnowmelt ron0tt peaksof 0.22 a cee fan was reached on May 27, 
1973 (1); and in Twisty Creek, a small basin on the north slope of the Mackenzie 
Mountains, a new summer rainstorm flood peak of 2.45 ih seeacne was established 
on July 6, 1973 (6). Estimation of 100-year floods using Thakur's model based 
on morphometric parameters has been improved thereby enhancing its predictive 
value (13). 

6. Use of Dendrochronology - The integrity of an oil or gas pipeline in a 
permafrost region may depend on the evolution of landforms by active geomorphic 
processes. Some insight into the age and development of particular forms and 


processes can be gained by use of dendrochronologic methods (5). 


= 12> = 


rie Semi-Permanent Snow Beds - Where proposed routing of pipelines is adjacent 
to semi-permanent snow beds, care should be taken not to affect the soil environ- 
ment. Soils in areas around the beds often approach an unstable condition due 
to snowbed runoff. Moreover, mass wasting may readily occur when summer rainfall 
events are prolonged and/or intensive (11). 
8: Active Layer Simulation - A simple site energy balance climatic simulation 
model may be used to estimate changes in active layer development and thickness 
(3). Such changes are caused by inadvertent or planned surface modification on 
highway and pipeline routes and may be foreseen for relatively little cost by 
use of this simple computer model. 

From the hydrological point of view, there should be further work 
on comparative studies of small coastal basins as an adjunct to the Beaufort 
Sea Program, there should be extensions of the semi-permanent snow-bed and the 
delta bank erosion projects, and there should be further examination of lakes 
and bogs which could be affected by the construction and Operation of oil and 
gas plants and pipelines. In our view a greater understanding of northern 


hydrology and maximum protection for the environment are concomitant goals. 
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ABSTRACT 


The effect of organic layer removal and drying of 
the soil surface on active layer depths was investigated 
by computer simulation, Simulations were performed with 
a standard set of site conditions for 51 locations in the 
permafrost region of Canada, using monthly mean 1971 wea- 
ther data. Contour maps prepared from the simulation re- 
sults illustrate the relationship between regional climate 
and the magnitude of active layer depth change due to sur- 
face modification, The effect of a surface wetness change 
is seen to vary with the climate of a region to a larger 
extent than does the effect of organic layer removal, due 
to the relatively high sensitivity of the surface energy 
balance to the wetness parameter change, The experiments 
illustrate the utility of one dimensional geographic point 


simulators in the consideration of spatial problems, 


I, Introduction, In cold regions, soil can exist ina 
frozen state from year to year, This condition is called 
permafrost, and its existence is defined precisely by the 
absence of soil temperatures of 0° Centigrade or higher 
throughout the year, Permafrost is widespread in northern 
regions, and, with the current interest in northern de- 
velopment, studies of the climatological, geomorphic and 
biological processes operative in permafrost regions are 
of increasing importance. The possibility of resource 
extraction on a large scale has focused attention on human 
perturbations in local and regional ecosystems. In addi- 
tion to ecosystem effects, construction associated with 
resource extraction is susceptible to damaging feedback 
from the disturbed systems. It is important that the na- 
ture of the interaction between climate, site, and eco- 
system be documented before decisions are made which pro- 
duce irreversible long term environmental effects and/or 
expensive engineering failures, 

One phenomenon of importance to local and regional 
ecosystems and to construction activities in permafrost 
areas is the annual evolution of the active layer, a thin 
surface layer which melts in response to energy input by 
solar radiation during the northern summer. Energy ab- 
sorbed at the ground surface is dissipated into the atmos- 
phere via conduction and evaporation, and infrared re- 


radiation, © A small heat. flux. into 2the: soil) ie generated 


by conduction and water transport in vapor and fluid forms, 
The amount of energy available to the soil anc the annual re- 
gime determine the depth to which the active layer: extends 
for a given soil, The available soil heat is in turn de- 
pendent on the surface energy balance in toto; and that 
energy balance is a complex function of the location's 
climate and site conditions (e.g. soil thermal properties, 
vegetative characteristics such as albedo and thickness of 
the organic mat), and their evolutionary histories, 
Construction related activities may alter the surface 
energy balance, changing the course of active layer melting. 
For example, removal of an insulative layer of surface vege- 
tation may increase summer soil melting, releasing soil 
water with subsequent thermal erosion (thermokarst) produc- 
ing differential settling. Brown et al (1969) reported 
doubling of active layer depths due to removal of the sur- 
face organic layer, The reader is referred to Gold and 
Lachenbruch (1973) for a review of the literature on thermal 
conditions in permafrost, These authors have pointed out 
the lack of quantitative knowledge both of actual tempera- 
ture structures in permafrost regions and of the influence 
of climatic and site conditions in determining those struc- 
tures, The importance of theoretical knowledge of the 
operation of natural systems has been stressed by Khilmi 


(1968), 


Most studies of physical-environmental effects on 


Oa 


active layers have been restricted to a limited number of 
sites, with restricted temporal climatic data, Most pre- 
vious modeling efforts have been of limited scope, due to 
the complexity of the processes modeled, the expense of 
modeling, and data processing and availability constraints, 
Models typically require time series data for several cli- 
matic variables, and the specification of dozens of site 
parameters (Goodwin and Outcalt, 1974; Nakano and Brown, 
1973). Sensitivity testing over such a multivariate space 
is quite expensive even when the voluminous observational 
data is available, When data requirements are reduced by 
structural changes within the model, "accuracy" generally 
suffers (Outcalt et al., 1974). In addition, lack of an 
adequate soil temperature observation net, particularly 
with uniform site conditions, makes model verification ex- 
tremely difficult. 

It was with the realization of these pragmatic con- 
straints that Goodwin (1972) developed a simple model of 
active layer formation which, although it is not capable 
of precise prediction, nevertheless embodies sufficient 
deterministic theory to allow reasonable estimation of the 
effect of climate and site perturbations on active layer 
formation over regions of interest. This paper reports 
the application of that model to the question of regional 
variance of process sensitivity to site parameter change. 


The work is essentially geographic and exploratory of a 


kt 


technique. It by no means sweeps the whole space of climate 
and site parameter variation, The study is limited to two 
site parameter changes and a temporal scope of one year, 
Nevertheless, it is felt that an experiment in geographic 
climatonomy will demonstrate the potential power of that 


methodology. 


II, Methodology. It is not the purpose of this study to 
predict precisely active layer depths at specific sites, 
but rather to assess in a semi-quantitative manner the spa- 
tial variability of the effects of selected surface modifi- 
cations on active layer melting, In order to do this, sur- 
face parameter values which represent reasonable end-point 
values of the natural continuous range of the parameters 
were used for simulation, 

Two site parameters were selected for testing: sur- 
face organic layer thickness and summer surface wetness, 
Of the many parameters that could have been selected for 
study, active layer evolution was previously demonstrated 
to be highly sensitive to these, and it is anticipated that 
they are significantly susceptible to modification by human 
activities, 

Four simulations were performed for each of 51 sites 
over the continuous and discontinuous permafrost region of 
Canada, The following parameter values were chosen for the 


Simulations: organic layer thickness -- 30 cm and 0 cm, 


summer surface wetness -- 1. and .5 (referred to a totally 
wet, freely evaporating surface with wetness = 1.). Contour 
maps of maximum active layer depth are presented ina later 
section for each parameter value combination, In addition, 
contour maps are presented which illustrate the effects of 
the parameter changes in terms of absolute depth differences, 
(Absolute depth differences are mapped since it is probable 
that the deleterious effects of increasing active-layer melt- 
ing increase in proportion to actual depth increase, rather 
than percentage increase, ) 

Simulations were performed using climatic data for 
1971 only, with all site parameters except wetness and 
organic layer thickness held constant over the area, The 
latter constraint isolates the effects under examination 
and simplifies analysis; the former is necessary for prac- 
tical reasons. Both constraints give rise to problems of 
interpretation of the simulation results. 

Since initial soil temperature profiles are not avail- 
able, the soil is isolated from its history. 1In particular, 
all sites are assumed frozen for 2 meters depth from the 
beginning of the run until soil thaw begins, a condition 
which is not, in general, true at the southern limits of the 
study, In addition, the constant parameter set chosen will 
not be valid at every point over the area, 

These facts constrain interpretation of the results, 


and it must be emphasized that the values presented on the 
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maps apply only for sites where the above conditions are 
realized, The representativencss of the results will an- 
crease with the extent to which the specified site parameter 
set is representative of the area to which it is applied, 
which will, in general, increase with latitude, The prob- 
lem of "representativeness" is little understood and is a 


fertile area for future geographic research, 


III, The Model, The computer simulation model used in the 
present study is an extension of a diurnal digital computer 
model developed by Outcalt (1972), which was inspired by 
the analog computer modeling experiments of Myrup (1969). 
Both models are based on the existence of an equilibrium 
surface temperature (discussed at length by Outed li 29739, 
which brings the summation of all surface energy transfer 
components (radiation, soil, sensible and latent heat 
fluxes) to zero, as required by the principle of energy 
conservation, 

The relationships which estimate the surface energy 
fluxes are non-linear and involve climatic variables, site 
parameters, and surface temperature, With all variables 
and parameters supplied or computed for a given time step, 
in part conditioned by information from the preceding time 
step, the equilibrium surface temperature is determined 
iteratively, thus providing estimates of the mean surface 


energy fluxes over the time interval. During soil thaw 
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refreezing, the value of soil heat flux is partitioned into 
components of heat content change and phase change of soil 
water (ice), from which active layer melting is computed, 

A detailed description of the model, with FORTRAN code, is 
found in Goodwin (1972). 

The model requires mean monthly values over the model 
year of the climatic variables listed in Table I, Curves 
are interpolated to these variables by the Fourier synthesis 
of monthly data points and sampled at 5 day intervals, the 
basic time step of the model, (January through December 
values are input, but the November and December values are 
used to precondition the soil temperature vector and to 
attenuate the edge effects of Fourier synthesis, The model 
thus operates from 1 November.) With the site parameters 
listed inwTable Il specified initially, simulations are 
performed for each combination of values of the parameters 
selected for testing. 

The primary advantage of the present model is its low 
operating cost, approximately 40 cents per year per simula- 
tion on the Michigan Terminal System, In addition, climatic 
data requirements are minimal, These properties enable ex- 
tensive simulation to be performed over a large area, The 
inclusion in the model of the set of site parameters impor- 
tant to physical processes allows extensive sensitivity 
testing on the effects of parameter change. The use of the 


equilibrium temperature strategy avoids the limiting require- 


es a 


ment of either collecting surface temperatures or assuming 
temperature curves, as did Gold et al. (1972). It must be 
realized that surface temperature is a product of site pro- 
perty variation and thus surface modification effects cannot 
be explored reliably with air temperature as the site forcing 
function, The temperature of the surface is part of the 
answer, not part of the question? 

The above results are achieved with obvious sacrifice, 
Some processes of undoubted importance are not modeled, In 
particular, snow melt is assumed instantaneous when modeled 
snow surface temperature exceeds zero degrees Centigrade, 
Previous work (Goodwin and Outcalt, 1974) has shown, however, 
that the timing of snow melt has little impact on ultimate 
active layer depths for small (<50 cm) initial snow depths, 
It is only when a snow cover extends far into the summer that 
it affects appreciably the total quantity of heat available 
for soil thaw, 

In addition to model specification errors, the use of 
long term mean values of meteorological variables introduces 
uncertainty in the estimate of atmospheric fluxes, due to 
the non-linear nature of the flux equations and the frequency 
dependence of energy transfer in a turbulent medium. 

Figure 1 compares modeled active layer melting with 
real depths computed from temperature observations reported 
by Weller (1973) for a site at Barrow, Alaska, The model 


underestimates freeze-back from the base of the active layer, 


although the course of thaw and the maximum depth attained 
are reasonably close, These effects are more accurately 
mimicked by more detailed simulators at the expense of at 
least a two order of magnitude increase of input information 


and simulation time (Outcalt and Goodwin, 1974), 


IV, Data, Climatic data for 51 stations over the permafrost 
area of Canada was obtained from the Atmospheric Environment 
Service, Canada, Monthly radiation data for 1971 was available 
for 14 of these stations, These 14 Seth one will be called 
primary stations, and the remaining 37, secondary stations. 

Radiation values for all stations were computed by 
polynomial regression equations derived from the primary set 
for each month, These equations involved the first three 
powers of latitude and longitude, Table III compares monthly 
values computed by the interpolation equations with avail- 
able values for two secondary sites, 

For the sake of consistency, all simulations reported 
were performed using the computed radiation values, Asa 
check, simulations were run for the primary stations using 
original data. These revealed no significant difference in 
maximum active layer depth due to the use of computed radia- 
tion values at those sites. Since the primary sites are 
well distributed over the entire region, it is expected that 
error due to this method of interpolation will be of impor- 


tance only at particular sites and will not appreciably al- 
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ter the final maps. 


V. The Maps. It is apparent from Figure 2 that active layer 
melting over the permafrost region of Canada is closely re- 
lated to the pattern of regional climate, which is greatly 
affected by regional air mass climatology as influenced by 
topography and the annual circulation regime, 

Of particular interest is the northerly shift of the 
contour lines in the area of the proposed Mackenzie Valley 
Pipeline Corridor. This shift is probably due to a combina- 
tion of factors such as the large quantity of heat trans- 
ferred northward by the Mackenzie River, the location of 
weather stations close to the river (although most stations 
across the map are located quite close to either lakes or 
rivers); the continentality of the region, and the movement 
of tropical maritime air. ‘For ‘comparison, consider Hays 
(1971) map of mean monthly precipitable water for July. 

Active laver melting maintains the pattern of Figure 2 
upon either removal of the organic mat or drying of the sur- 
face (Figures 3 and 4), Both changes operating together 
amplify the pattern (Figure 5), When the active layer depth 
maps: are subtracted. one from the other, spatial variation of 
process Sensitivity to parameter change is clearly seen, 

In Figure 6, the effect of organic layer removal is 
isolated, and it is seen to be relatively independent of 


spatial position, Over the entire area south of the 30 cm 
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contour line there is little systematic variation, Only 
south of 65° N latitude is a small latitudinal dependence 
found, 

A more pronounced spatial dependence is seen in Figure 
Veewnneh tsOlaves the ertect of lowering summer surface wet- 
ness, he ettecr of the “two parameters operating in tandem 
is depicted in Figure 8, 

These results are quite understandable, Both para- 
meter changes are effective primarily in the summer. Soil 
surface wetness is important only when the soil is snowfree, 
The organic layer acts as a barrier or low frequency by-pass 
filter to the flow of heat from the soil surface, Removal 
of the layer allows more of the summer heat pulse to pene- 
trate the deep soil, whereas in winter the effect on energy 
transfer is masked by the presence of the snow blanket, 

Soil heat flux is a relatively small component in the 
summer surface energy balance, and removal of the organic 
layer, maintaining the soil in a wet condition, has little 
effect on the surface equilibrium temperature, Thus, an 
increase in active layer depth due to organic layer removal 
is tied relatively weakly to the climatic determinants of 
the surface energy balance, 

Changing surface wetness, on the other hand, affects 
active layer melting by changing the rate of evaporation, 
and therefore the surface equilibrium temperature, There 


is then, through the linking mechanism of the surface energy 
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balance, interaction between a change of surface wetness and 
such climatic elements as relative humidity and wind speed. 
Since these elements have spatial variance, and since the 
interaction between them and surface wetness is non-linear, 
the effect of a surface wetness change is non-linear as well, 
and exhibits strong spatial variation, 

Table IV displays the results of the simulations along 
a line from Inuvik in the Mackenzie Delta area to Ft. Simpson. 
It is apparent that artificial crainage. of -areéac oispeaty, 
soils along the line would have a considerable effect on ac- 
tive layer melting, although the effect of drainage on the 
thermal properties of the peat layer might lessen the total 
depth increase. (This relationship is not modeled presently 
but will be explored in future work with a simulator contain- 
ing a more refined treatment of soil thermal properties. ) 
Complete removal of the organic mat produces an equivalent 
active layer depth increase, and both removal and drainage 
(a probable result of pipeline construction) nearly doubles 


the maximum depth of the active layer, 


VI. Summary and Conclusions, In summary, the research re- 


ported in this paper has demonstrated the feasibility of 
employing a simple site energy balance climatic simulator to 
a spatial problem, The effects on active layer melting depth 
of surface modification have been determined quantitatively 


for a region where the active layer is an important factor 


in human settlement and resource extraction, 

The simulations have shown that stripping the organic 
layer from mineral soil produces an appreciable increase in 
the depth of active layer melting, and that this increase 
is little affected by regional variation in climate, Change 
in surface wetness, on the other hand, produces an effect 
which is sensitive to regional climatic variations, 

Interpretation of these results is necessarily con- 
strained by the theoretical nature of the work, In parti- 
cular, the problem of "representativeness" must be kept in 
mind, 

We have in effect constrained several axes of a multi- 
dimensional variate space in order to isolate the effect of 
a few parameters, This methodology is by no means novel, 
and is in fact the classical method of sensitivity testing. 
The importance of the work stems from the fact that the 
study has been conducted at relatively low cost by the use 
of a simple computer model and readily available weather 
data, Much future work on permafrost and its interaction 
with site and climate will no doubt involve more sophisti- 
cated models, and these will require extensive data collec- 
tion, The climate and site-parameter space will be investi- 
gated by many workers in great detail and over large time 
spans, This will be an expensive and difficult undertaking, 
and it is important that preliminary work be done with simple 


models and available data, so as to approach the multi- 
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variate problem with an abstract knowledge of basic system 
sensitivities. It is felt that the methodology reported 
here, and similar modeling activities, Gaal be of some 
assistance to that end, 

Planned future work at the University of Michigan will 
address the problem of temporal variance of the spatial pat- 
tern of active layer response to surface modification, It 
is hoped that this effort will lead to a better understanding 
of both the effect of past climate on present conditions and 
the temporal decay of the effects of various surface modifi- 
cations. In addition a model is being developed which will 
produce time-distance maps of surface modification effects 
for any given pipeline route, This model may be useful in 
pipeline routing decisions, 

The authors wish at this point to stress the need for 
the establishment of a net of soil temperature observation 
Sites over the permafrost areas of North America, Soil tem- 
peratures and active layer depths should be monitored at 
each location at sites with soil and surface conditions com- 
parable over the net, These observations would be of great 
value in model building, and are necessary for verification 


of the theoretical work in progress, 
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TNPUT DATA: CLIMATIC VARIABLES 
(12 Monthly Mean Values: Jan - Dec) 
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Total Solar Radiation 

Air Temperature at Instrument Height 
Wind Velocity at Instrument Height 
Barometric Pressure 

Relative Humidity 

Cloud €over Fraction 
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INPUT DATA; SITE PARAMETERS 


PARAMETER VALUE 


a 


Thermal Properties (From Nakano and Brown (1972)) 


Material Conductivity | Das Se 
(cal/em sec °C) (em yao 

Snow 0,00020 0, 00120 
Frozen Organic Layer 0,00028 0.00046 
Thawed Organic Layer 0,00067 0.00070 
Frozen Mineral Soil 0700220 0.00420 
Thawed Mineral Soil 0, 00206 0,00370 
Other Parameters 
Organic Layer Porosity 05.50 
Mineral Layer Porosity 0,42 
Snow Depth ian to me it. 40: em 


Surface refreeze to 1 Jan: 0 to 40 cm 
(linear increase) 


Snow Surface Wetness 0 
Instrument Height 150 ecm 
Snow albedo 0,80 
Summer Surface Albedo OVALE 
Snow Surface Roughness Length 0, 0 cm 
Summer Surface Roughness length PAOLO e Toiir 
Organic Layer Thickness 30 or 0 cm 
Summer Surface Wetness 0 eC re Ono 
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TAS LE er 
MONTHLY MEAN TOTAL SOLAR RADIATION: 


COMPARISON OF INTERPOLATED VALUES WITH OBSERVATION 
(Observation source: DOEAES - Monthly Record ) 


LY /OAY 

MONTH PT, NEGSON Ay BoC, EUREKA, N.W.T, 

( ) = Observation value 
Apr ite 282 267) 
May SLO. 500) 527 (540) 
Jun 469 (451) 562 (5e)) 
Jal 482 (481) 42s (449) 
Aug 72} Gao 222 (2h) 
Sep POI 226) LOS eo 2s) 
Oct Lee) S. (iy) 
Nov 47 (8) 5 FO) 
Dec 2 28) rare 
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TABLE LV 


RANGE OF ACTIVE LAYER DEPTH AND INCREASE 


ALONG A LINE FROM INUVIK TO FT. SIMPSON 


_ CC CC CC rcC1V3O3O i a ee 
ee 


ACPIVE AYER DEPTH Cem) 


— 


Summer Surface Wetness 
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Organic 30 cm Pie LS 1O5 = TH6 
Layer 
Thickness O cm Ore di 7 Whe eras be lal 


ACTIVE LAYER INCREASE (cm and %) 
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CAUSE ABSOLUTE INCREASE % INCREASE 
Organic layer removal B36 30. eS 
Decreasing Surface wetness 31 - 35 20-4 
Both C7 bay 70 xO), ee vel 
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Abstract 


Applications of dendrochronology are discussed with reference 
to dating geomorphic events which are caused or accelerated by 
ice in permafrost and to reconstructing the past climate and pre- 
dicting future climatic trends. It is suggested that the recon- 
struction of climatic variations calculated from ground temperature 
profiles can be corroborated by interpretation of temperature 


anomalies calculated from tree-ring data. 


Rees 


is Introduction 


With the accelerated rate of development in permafrost 
regions and the anticipated construction of pipelines there is a 
growing need for thorough knowledge of permafrost. This cannot 
be achieved by one discipline, but by a concerted effort of 
multidisciplinary research including dendrochronology, dendro- 
climatology, and dendrohydrology. 

Considerable progress has been achieved in many fields 
of permafrost research in the last decade but there are 
several where knowledge is inadequate. Two are: dating of 
geomorphic processes in permafrost (Brown and Péwé, 1973) includ- 
ing cryogenic processes which control hydrology and create 
landforms characteristic of permafrost regions; and studies of 
past climate which can be used to predict the trend of future 
climate. In both these fields, which are of primary interest to 
pipeline construction, dendrochronology can make a significant 
contribution. 

Recent new techniques in using a computerized system 
produce a variety of tree-ring data with speed and accuracy 
which could not be attained some years ago. In Canada, tree-ring 
research is carried out at the Western Forest Products Laboratory, 
Vancouver, B.C.(Parker, 1972; Parker and Jozsa, 1973) and at 


Brock University, St. Catherines, Ontario (McNeely et al., 1973). 
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The Glaciology yeesha of the Inland Waters Directorate 
applied dendrochronological data to the study of floods and ice 
jamming on Liard and South Nahanni Rivers and part of the 
Mackenzie River (Henoch, 1973; Parker and Jozsa, 1973). In this 
report application of tree-ring data to permafrost studies is 


discussed. 


Le The Influence of Vegetation on Permafrost and of 


Permafrost on Vegetation 


The mutual relation of vegetation and permafrost is 
complex. Briefly, it can be said that vegetation protects perma- 
frost from thawing and that permafrost inhibits the growth of 
vegetation although examples of many exceptions are quoted by 
Brown (1963) and Péwé (1966). 

Low temperature and the permafrost's impermeability 
retard the growth of the root system and force roots to grow 
laterally. Benninghoff (1952) observed that the side of the 
root in touch with permafrost has restricted growth and that 
this is reflected in ovoid shape of the cross-profile of the 
KOOU. 

If, however, trees and brush are removed, causing a 
disturbance to the surrounding moss and peat, degradation of the 
underlying permafrost follows. These changes may have significant 
local effects on the hydrology of the region and thus on the 


proposed construction of pipelines. 
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Study of anomalies of tree-ring growth in permafrost 
can be a sensitive method of determining the past climate of the 
area as any interference with growth of the tree can usually be 
dated by analysis of its growth rings. Tree stems and root 
samples have been collected in the Mackenzie Delta to determine 
whether X-ray of the samples can provide information on the 
influence of permafrost on tree growth. Preliminary results 
(Parker et al., 1973) show a reduction of growth in tree-rings 
of roots of an excavated tree stump, probably owing to the rising 
level of permafrost. As the influence of other factors has to 


be taken into consideration further studies are needed. 


oe Geomorphic Processes Caused or Accelerated by the 


Presence of Ice in Permafrost 


Dendrochronology is well suited for dating landforms 
and events originating from geomorphic processes in permafrost. 
The application of some of the techniques used in the field 


has been described by Henoch (1973) and Parker and Jozsa (1973). 


36 Cryogenic processes 


Cryogenic processes may lead to an increase or 
decrease of volume of ground ice. As the result of direct or 
indirect action of these processes, trees growing on permafrost 


landforms may be injured, destroyed, or have growth aberrations 
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that can be used for dating and other studies based on tree-ring 


chronologies, 
3.1.1 Cryogenic processes leading to the increase of ground 
ice volume 


Sit Pei ek Pingos 


Pingos overgrown with coniferous trees are found at 
the foot of Black Mountain (Aklavik Range, Richardson Mountains). 
Some aspects of their formation could be documented and interpreted 
by the systematic study of the annual growth of tree-rings. 
Homes et al., (1968), for example, estimated the age of the 
youngest pingos in Alaska using the age of trees and changes in 
the tree communities. 
3.1.1.2 Ice wedges 

The age and the formation of ice wedges can be 
determined if they are found to have affected the trees. A 
buried tree, split by an ice wedge, was found near Point Separation 
in the Mackenzie Delta in 1973. After the determination of the 
age of this tree by tree-ring and radiocarbon dating the minimum 
age of the ice wedges can be established. 
Seel lies Aggradational ice 

Aggradational ice forms when the upper permafrost 
surface gradually rises over a period of years and thereby 
incorporates the ice lenses formed at the base of the active 


layer (Mackay, 1972). Its origin depends on type of soil, 
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moisture content, vegetal cover, and temperature profile of the 
active layer. 

It is suggested here, that the age and the rate of 
growth of this type of ice can be studied, in some cases, by 
dendrochronological methods. Two examples where tree-ring 
studies can contribute to understanding of the formation of 
aggradational ice are: | 
3.1.1.3.1 Mounds of aggradational ice - Ice lenses of perennially 
frozen ground in Mount McKinley, National Park in Alaska, studied 
by Viereck (1965) appear to be lenses of aggradational ice. They 
are mounds on which white spruce grow. Viereck suggests that two 
factors contribute to aggradation of ice under the trees: deep 
frost penetration, through comparatively thin snow cover under 
the canopies of trees; inhibition of heat exchange during the 
Summer by a thick mat of moss on the mounds.The growth and 
collapse of the mounds causes trees to lean and topple over. 

These events could be dated by tree-ring analysis. 

3.1.1.3.2 Humps of aggradational ice - Humped topography has 

been observed in the Mackenzie River Delta, on the terrace of 

Peel River near Indian Village. Preliminary survey of these 

humps, which have the appearance of incipient pingos, or palsas has been 
carried out during the summer of 1973 (Henoch, Parker, and Outhet, 

in preparation). 


The humps are oval but irregular in shape, Their 


2550.5 


diameters are 40m, 45m, and 80m, respectively. The largest one is 
4m high above summer river level. They are skirted by a drained 
lake and flood channels, The southern side of humps is undercut by 
the Peel River. A cut-bank of the Peel River exposes the cross 
profile of the humps displaying fine silts with high ice content, 
horizontally stratified under the channel. However, the stratifi- 
cation under the humps appears to have been dislodged from its 
horizontal position and assumes an anticlinal structure under each 
hump. 

The origin of these humps is not certain. It is 
believed that aggradation of ground ice governed by the interplay 
of flooding and zonal distribution of vegetation are important 
factors in their formation. 

Intensive sampling of trees growing on the humps was 
carried out for statistical analysis of tree-ring data, to 
determine the influence of the climate, depth of active layer, 
and flooding on the rate of growth of these trees. 

Further research should be directed toward investigation 
ohiee 

(1) Aggradation of permafrost due to insulating and moisture 
retaining influence of caribou moss. 
(2) Ground ice under the humps, its growth and heaving 


effect. 
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3.1.1.4  Palsa and string bogs 


Palsa and string bogs occur in the southern fringe of 
permafrost and trees seldom grow on them, Palsa, in many ways 
Similar to pingos, are low knolls of frozen peat usually no more 
than 2m high. They are formed by combined action of peat 
accumulation and ice segregation in the underlying mineral soil 
(Brown, 1970). | 

String bogs are strips of vegetation in shallow 
depressions generally filled with water. Their origin is not 
certain but most researchers think that cryogenic processes play 
an important part in their formation. 

Dendrochronological methods have been applied to study 
the formation of string bogs in Finland (Yli-Vakkuri, 1960; 
Eurola, 1962). Aartolahti (1965) studied leaning trees growing 
on string bogs and determined annual growth of tree diameter and 
its eccentricity as a proportion of radial growth in the direction 
of lean. As the lean of trees was against the prevailing winds 
it was apparent that the tilting of the trees was caused by 
heaving of the ground and the heaving was caused by ice pressure. 

These data were thought to represent changes in the 
intensity of cryogenic processes. An attempt to correlate the 
data with climate was unsuccessful partly due to the paucity of 


information, 
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Sv lage Cryogenic processes leading to the decrease of ground ice 


volume. 

Examples of landforms originating due to cryogenic 
processes and leading to a decrease of the volume of ground ice 
are thermokarst and mass movement. Both can be studied by 
dendrochronological methods. 
3.162.) ~ Thermokarst 

This topography owes its name to its similarity to karst 
topography. Melting of ground ice causes land subsidence and 
produces hummocky terrain with generally level surface and steep- 
sided depressions which are filled with silt and meltwater. Lakes 
and ponds, with some dead trees in the water and some living trees 
leaning toward the water, are found frequently in the forested 
regions of permafrost. The trees indicate the recession of the 
shores of thermokarst ponds and can be used to determine the rate 
of recession of the shores. For example, Wallace (1948) has 
applied dendrochronological methods to study the formation of 
thermokarst lakes (cave-in lakes) in Alaska. 
3.1.2.2 Mass movement 

Mass movement in permafrost regions manifested by 
rotational slips, slides, slumps and solifluction is often caused 
or accelerated by thawing of ground ice. A number of landslides 
in the permafrost region of the Mackenzie and Liard River valleys 


have been studied by using dendrochronological methods. The 
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results have shown the usefulness of these methods in studying 

(1) the minimum age of the landslides; (2) periodicity of poly- 

genetic causes of movement; (3) Synchronous occurrence of many 

of these movements in various landslides (Parker, Jozsa, 1973). 
Determination of the age of the landslides provides 

information on how long the landslides induced by activities of 

man would likely be active. As a significant number of land- 

Slides has been reactivated synchronously in different localities, 

it may be possible to correlate landslide occurrence with climate. 


More data, however, are needed. 
4. Dendroclimatology and Permafrost 


Rings of trees can provide a unique source of information 
on past climate as far back in time as tree-ring chronologies can 
be established. Although old trees near the northern tree line 
are rare, nevertheless the author has found trees over 500 years 
old in the Mackenzie Delta. Giddings (1962) has reported trees 
about 1000 years old growing above the Arctic Circle. 

Tree-ring width patterns over large regions can be 
calibrated with the climatic data and the established relationship 
can be used to extrapolate back into the climatic past (Fritts et 
al., 1971). The studies have demonstrated that tree ring data 


can be used to reconstruct past surface pressure anomalies, storm 


hoe 


tracks, temperature, precipitation or any other climatically 
related data that are available for calibration with tree-ring 
data. Fritts (1973) carries out investigation to reconstruct 
yearly seasonal anomaly patterns of climate (including temperature) 
and tree ring indices back to 1500 A.D. 

Besides the variations in width of tree rings, the 
density of dated tree-ring series can also be used to obtain past 
climatic records (Parker, Henoch, 1971). For the climate studies 
of permafrost regions the variations in tree-ring density might 
prove to be a more sensitive parameter than the tree-ring width. 

While in the arid zone precipitation is the main factor 
affecting the growth of a tree and the temperature plays a 
secondary role, in the Subarctic the effect of temperature is 
more significant than that of precipitation. 

It is known, from the studies of ground temperature, 
that any change of temperature at ground surface is propagated 
downwards. Cermak (1971) has shown that a detailed record of 
ground temperature for varying depth can be used to reconstruct 
past climatic history, 

The anomalies of ground temperature profiles could be 
correlated with temperature fluctuations as derived from 
dendroclimatic studies. To achieve this, however, the temperature 
profiles from deep drillings and more good tree-ring chronologies 


from the permafrost area or vicinity are needed, 
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5. Conclusions 


The review of literature and the author's experience 
indicate that application of dendrochronological methods to the 
study of landforms, hydrology, and climate can contribute to 
solving various problems arising during the planning of pipeline 
construction and maintenance. 

A systematic effort should be made to collect tree-ring 
samples along the northern tree line to build dendrochronologies 
which will provide the basis for dendroclimatic studies of perma- 
frost regions. The progress made in the last decade in the 
development of new techniques to obtain tree-ring data and 
methods to interpret them justify a wider application of 


dendrochronology to studies of permafrost. 
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ABSTRACT 


Aerial surveillance and photography during the 1972 and 
1973 Mackenzie River break-up periods were carried out to monitor 
channel ice effects and to map water surface velocities. Some 
visible effects of ice action are discussed and illustrated with 
reference to the problems of navigation, highway and pipeline crossings, 
and settlement location. Surface water velocities are mapped for 
some sections to illustrate the possibilities of the method and its 
relevance to river scour ;,to the movement of spilled oil and to general 


hydrographic work. 
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Channel ice effects and surface water 
velocities from aerial photography of 
Mackenzie River break-up 
by 
Di Ke MacKay, D.A. Sherstone and K.C. Arnold 


Introduction 

Aerial photographs of Mackenzie River spring break-up were 
taken in 1972 and in 1973 to expand the information base available on 
northern river ice disintegration processes and to test methods of mapping 
surface current velocities. Previous to 1972, only scattered air and ground 
observations were made with consequent limitations on the interpretation 
of ice break-up. 

A continued exan*nation of Mackenzie break-up (see selected 
reference list) is warran ed for a number of reasons. From the econanic 
point of view, the river forms a vital transportation link between north 
and south. Its use as a waterway, however, is severely limited by the 
shortness of the open season. More information on break-up could conceivably 
improve its utilization b: ietailing progress of break-up and by revealing 
changes in the deep water . annel caused by ice movement. 

The economic argument may become very important in the near future 
if northern oil and gas exploration activities are successful and production 
decisions are reached. Pipelines may then be built with consequent heavy 
demands on river cargo space that could tax the capacities of freight 
handling facilities on the Mackenzie. A complicating factor is the 
possible increase in movement of road—-building equipment and supplies for 
construction of the Dempster and Mackenzie Highways. Eventually, these 


roads will provide all-weather, year-round alternatives to shipping via 


pe 


water during the short season on the Mackenzie but, in the interim, it 
seems clear that full use of the river is important for northern econanic 
development. 

In an engineering sense, break-up surveillance is important 
for the information it provides for highway and pipeline planning, 
construction, and maintenance. Bank erosion and bed scour caused by 
ice movement and ice jamming can affect location and design of pipeline 
and highway river crossings as well as river docking facilities. Ice 
jams can cause localized high water on the main Mackenzie and backwater on 
its tributaries. This, in turn, may induce washouts at pipeline or road 
crossings of rivers and streams. It seems, therefore, that for both 
economic and engineering reasons, continued scrutiny of break-up processes 
on the Mackenzie River is warranted. 

The general pattern of Mackenzie River break-up is well known 
(17, 15, 12). In summary, radiational melting of snow begins along banks 
and on slopes in early April causing a rise in river stage and development - 
of shore leads. As runoff from tributaries flowing into the Mackenzie 
main stem increases, transverse leads open across the full width of the 
river (Figure 1). A further development of this stageis the intermittent 
downstream movement of sections of the main ice cover. The movement of the 
main ice cover is partly dependent on break-up of the Liard River, the 
major Mackenzie tributary. 

When Liard ice moves, in April or early May , its streamflow can 
exceed the Mackenzie's which it joins directly above Fort Simpson. 
Normally the rising Liard flow pushes ice wreckage (Figure 2) high 
onto channel banks with remnant ice entering the Mackenzie and stopping in 
the river stretch opposite Fort Simpson and below. If the stage continues 
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to rise, ice in the Mackenzie moves downstream leaving the river open at 
the settlement. Often the Mackenzie ice cover above the Liard junction 
(Figure 3) remains intact for some days thereafter. 

During the 1972 break-up and through the early summer season, 
extremely high water conditions were experienced on the Mackenzie River and 
in its Delta. With these high water levels, a few large ice jams developed, 
two of which were of some consequence. One jam, below Fort Simpson, caused 
flooding in part of the settlement; and the other jam, at Fort Norman was 
sufficiently high to push or spill ice into the ee rim of the village. 

In the Delta, masses of ice were deposited on levees and some central-outer 
portions below Reindeer Station were covered by floodwaters. 

In contrast to 1972 conditions, water levels in 1973 were reasonably 
normal with many small or partial jams developing (Figures 4 and 5). Only one 
fairly large size jam occurred, that being near the proposed Canadian Arctic Gas 
Pipeline crossing above Point Separation. In other respects, the progress of break- 


up was normal. 


Channel Ice Effects 
Ice moving in a river during break-up erodes channel islands, beds, 
and banks. Many of the erosional effects on these features can be documented, 
at least in part, by aerial surveillance. 
The heads of many channel islands are subject to strong ice push 
(Figure 6) and those of alluvial origin are probably subject to annual modification 
in shape. Others of bedrock or till are less affected, the latter by forming 
a protective boulder veneer on their upstream ends. 
Channel beds are open to abrasion or corrasion by the tumbling action 


of ice wreckage and by ice grounding in shallows. Along some river stretches, the 
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planing action of ice has formed pouider pavements (Figure 7) visible at low 

water. Above the pavements and along the banks of many river reaches, there is 
evidence of gouging and redistribution of bank materials by ice shove action. 
Often, banks are topped by ice shove ridges. Such features are generally 

developed by initial onshore ice movement. Any subsequent ice shove or pressure 

is attenuated by the facade of stranded ice as implied in Figures 8 and 9. 

Some of these features my not be detectable on aerial photos but their information 
can be insinuated at locations where ice stranding or jamming is observed. 

The most important of channel erosion effects are due to ice jamming. 

Ice jams affect channel morphology in a number of ways. The increase in head 

with flow moving under the jam can cause appreciable bed scour. In some instances, 
a jam in the navigation channel will direct flow into an ancillary channel, 
deepening it and changing the distribution of flow, (Figure 10). The redistribution 
of flow may be accentuated by jamming in subsequent years until the ancillary 
becomes the main channel. 

Partial jams can also cause marked changes in channel cross-sections | 
through flow redirection, erosion and deposition. During partial ice jamming, 
flow is concentrated in a particular part of the cross-section. This flow con- 
centration could cause more scour than expected by normal evaluation methods. 

If this occurred at underwater pipeline crossings or bridge locations, both 
piers and pipeline could be undermined or exposed. 

Backwater on the Mackenzie's tributaries and the main stem caused by 
ice jamming is a factor in the Planningof highway and pipeline routes. Moreover, 
it is of practical importance in settlement location and in design of water 
supply and sewage systems. Aerial photography is certainly the best method of 


documenting the effects of backwater including the flooding of settlements. The 
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settlements of Fort Simpson, Fort Norman, old Fort Good Hope and Aklavik on 
the Mackenzie have all been subject to floods in which ice played an important 
role (9, 13). Aerial photography following the 1972 flood recession is 
available for Fort Simpson. At Fort Norman, there are photographs of both the 
early and late stages of the 1972 ice jam in which some ice entered the lower 


village. 


Surface Current Velocities 

The analysis of surface flow over a period of time may aid in 
distinguishing scour and fill in critical river sections of interest to pipelines 
and highways. The reasons for this possibility are that changes in channel 
roughness, shape and gradient should affect water surface velocities and, Secondly. 
these velocities can be mapped from air photos. Two ways of estimating surface 
current velocities from air photos are examined and discussed in this report. 

In the first instance, discrete blocks in moving ice wreckage are 
located with respect to channel shore control points on consecutive air photos. 

If the aircraft speed and the interval between exposures is known, then the velocity 
of ice movement can be calculated as illustrated in Figure 11. One of the 
limitations of this method is the large number of plotted movements that must be 
measured to produce flow lines (Figure 13) over a given stretch of river. 

An alternative method is to contour the apparent stereo parallax (Figure 12 ) 
that results from water surface movement over the interval between photo exposures. 
The important criteria, as indicated in previous studies (4,5,6,7,19), are as 
follows: 

1) Visible foaming, sedimentation or artificially introduced discrete control 
points in the water to permit measurement of the moving surface water; 


2) Shore control points to permit " tying in " of water data to a rigidly controlled 
shore datum; 
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3) Rigid control and recording of the interval between photographic exposures 
to permit precise determination of displacement. 

In the Mackenzie study ice wreckage combined with substantial 
volumes of sediment and flood debris satisfied the first criterion. Adequate 
shore and water control, however, was often lacking due to the mismatch of 
flight line direction with river orientation. Moreover, excessive pitch and 
yaw was often discernible when photo transparencies were examined in a precision 
stereo plotter. This was particularly evident on photos taken during times when 
the aircraft was following bends in the river. The initial 1972 study, therefore, 
was restricted to two test sites near a possible pipeline river crossing and 
where, fortunately, high quality photography was available. The river at these 
sites also contained a considerable amount of suspended sediment, discrete ice 
blocks, and some foam from turbulent flow. 

Contour plots of two 1972 test sites, located at miles 607 and 612, 
approximately 20 miles above the Sans Sault Rapids, are shown in Figures 14 and 
15. These two sites represented the first experimental tests to determine the 
usefulness of this technique and were the only test sites which could be completed 
in the plotting time available. A visual comparison of the contour and discrete 
ice block measurement plots shows a reasonably good correlation between the two 
techniques. <A further comparison with Canadian Hydrographic Service Chart 6420 
suggests close agreement of the high flow velocity paths with the two alternate 
navigation channels. 

In view of the encouraging results obtained with the samples tested in 
1972 it was decided to examine further 1972 aerial photography in this manner and 
to then examine these same areas during similar 1973 break-up stages. Photography ¢ 
the 1973 break-up is limited as ice clearance in the upper reaches of the river 


(upstream of Wrigley) occurred prior to the arrival of the photo aircraft. 
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Imagery of moving ice debris in the areas of miles 607 and 612 was 
again obtained, but the relatively low concentrations of moving ice and lack 
of turbidity prevented continuous contouring of the surface (Figure 16). The 
area in which contouring was completed proved too fragmentary to permit a full 
comparison of 1972 and 1973 data. Values derived from discrete ice block 
measurements indicate that velocities were generally similar in both years; the 
average derived flow velocities using this method were 4.07 and 3.45 Me. He LO 
the 19/2 and 1973 tests respectively. These values may be suspect due to a 
decreased density of ice " wreckage " in 1973, but the stage-discharge throughout 
the 1972 break-up season exceeded that in 1973 which tends to support the reported 
yelocities. 

Im order to provide sites for future comparison of surface velocity 
patterns, additional 1972 photography was examined and stretches were selected 


near Wrigley (2 sites) and Camsell Bend. These stretches (Figures 17, 18 and 19) 


provide ideal test sites for contouring since the river is both narrow and 
straight for a distance of roughly 20 miles. Of the Wrigley sites the southern 
stretch (Figure 17 is the most interesting since it covers an uninterrupted 
distance of ten miles and includes a channeling of the flow due to the occurrence 
of fast ice (Mile 364) and a resultant " compression " of contour spacing. This 
indicates a rapid rate of velocity change due to channel constrictions (this 


effect also occurs at Mile 370). The section covered by the test plot also 
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includes a zone of back eddies (indicated by the " negative " or " depression 
contours) and a large area of significant surface shear in the Y plotting direction, 
which causes increased residual parallax and the resultant loss of a contourable 
surface. 


Over the 10-mile Wrigley stretch (Figure 16), ice wreckage is evenly 


distributed with block sizes grading from large at the down-stream end to small 
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at the upstream one. The only area of open water is just below the constriction 

at Mile 364 which poses some difficulty in carrying out the contouring. The 
Canadian Hydrographic Service Chart 6414 covering the plotted section shows sand 
deposits and shoals at the ice covered constrictions. These are normally exposed 
during periods of low summer flow. Maximum surface flow velocities at Miles 364 and 
370 (Figure 18) coincide with the deep water navigation channels shown on the 
hydrographic charts. 

The contour velocity values in the 1973 test sites were derived from 
the discrete ice block measurement system (4,5,6,7). In 1972 the calculation of 
velocity contour values was initially made using a similar triangle solution and 
the discrete ice block measurement system as a check. 

The results to date indicate the following: 

1) Maximum ice wreckage concentrations occur in the areas of maximum velocity. 
Examination of this data using a Spatial Data color display densitometer lends 
eredence to this evidence (20). 

2) It is possible to contour the velocities of a complete water surface using 
photogrammetry. This suggests the possibility of relating changes in surface 
patterns to changes in bed conditions. 

3) Surface water velocity plots could be a significant tool in assessing 
distribution and movement of oil following spills, and 

4) Such plots may be beneficial in staging hydrographic work on rivers by 
indicating changes in the navigation channel. 

In summary the use of aerial photography during particular stages 
of spring break-up can provide useful data on changes in channel morphology 
caused by ice action and on the opening of sections of the river for navigation 
purposes. Possibly, of greater importance is the application of such photography 


to oil spill, bed scour and other hydrographic problems through the provision of 
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a graphic portrayal of surface water conditions. 
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Figure 1 
Mackenzie River at Wrigley showing transverse cracks 
across the river at Hodson Creek and Wrigley River. 


Figure 2 
Liard River ferry crossing site on the Fort Simpson side, 
May 5, 1973. Note the ice wreckage stacked approximately 
10 metres above the river level. 


ie 


Figure 3 
Junction of the Liard and Mackenzie Rivers on May 6, 1973. 
The ice cover on the Mackenzie above the JUNC oLOn. 1S Stall 


intact. In the foreground, silt-covered ice wreckage from 
the Liard litters the water approaches to Fort Simpson. 
be ull —— 


Figure 4 


Minor ice jamming adjacent to Normal Wells settlement 
on Way 13; 1973, 
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Figure 5 
Ice concentration in the deep water channel below 


Fort simpson, May 7; 19732 The dark sediment=laden 
1cE sis sirom= the Liard River. 


Figure 6 
Ice push on the head of an island in the Mackenzie 
near Mile 225. In the background, an ice ridge may 


be seen traversing the main channel. 
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Figure 7 
Boulder paving on the Mackenzie River. The planing action 


of moving ice during break-up is partly responsible for 
development of these forms. 


Figure 8 
At Fort Simpson (May 6, 1973) stranded shore ice acts as a 


buffer limiting erosion by the mass of ice wreckage moving 
downriver. The dark area of jumbled ice, encompassed by 


Shear lines, is from the Liard River entering the Mackenzie 
from the right background. 
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Figure 9 
Shear lines caused by moving ice near the cliffs of 
the Remparts, May 16, 1973. 


Figure 
Ice jamming the navigation 
some deflection of flow to 
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channel at Mile 470 causes 
the left around the islands. 
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Figure 12 
Stereo pair illustrating false parallax caused by ice 
movement, Mile 610, Mackenzie River, May 19, 1972. 
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Figure 13 


Discrete ice block movement test plot, Mile 612, 
Sans Sault area (1972 photography). Vector arrows 
indicate downstream movement of each block between 
photo exposures. 


MACKENZIE RIVER A-7 TEST SITE ¥i58 
MILE 612 (CARCAJOU RIDGE AREA) 
HORIZONTAL SCALE: APPROX. |:8,850 
VERTICAL SCALE: APPROX. 1:8,333 
CONTOUR INTERVAL: REAL: 1.0mm. 
APPARENT: 65 FT. 
PHOTO INTERVAL: 23 SEC. 
DASHED ‘CONTOUR’ LINES INDICATE AREAS IN 


WHICH RESIDUAL PARALLAX REMAINS AFTER 
CLEARING REMAINDER OF MODEL. 


HORIZONTAL SCALE 116,850 


FEET 600 5001000 __ 1500 2000 FEET 
METRES 100 oO 200 °«3300«=~=«00=S*=«iDS:«OD METRES. 
VERTICAL SCALE 1:8,333 
FEET $00, 500 2000 FEET 
METRES 100 0 100 200 300 ~ 600 METRES 


=o 


Figure 14 


Surface contouring, Mile 607 Sans Sault area, 1972 
photography. Values shown for each contour level 
are machine intervals and not velocity values. 
Velocity values derived as described in text. Dashed 
lines indicate area of slight plotting difficulty due 
to small amounts of residual parallax. 
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Figure 15 


Mile 606, Sans Sault area, 1972 photography. The 
first test site attempted using photogrammetric method 
of stream velocity contouring. Excessive parallax and 
stereo gaps prevented further work at this location. 
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Figure 16 


Surface contouring, Miles 608-610, 1973 photography, 
Sans Sault area. Contouring is fragmentary due to 
lack of turbidity. Photography of this stretch of the 
river was acquired later in the break-up than previous 
years’ photography. Contour values are shown in m/sec. 
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Figure 17 


Surface contouring, Mile 375, Wrigley area (northern 
site), 1973 photography. 


SS MACKENZIE RIVER A-7 TEST SITE 73-1 
MILE 375 (OCHRE RIVER AREA) 
CONTOUR INTERVAL: APPROX. 


NOTE: ‘DASHED! CONTOUR LINES INDICATE AREAS 


OF SLIGHT RESIDUAL PARALLAX PREVENTING 


\cé COMPLETE CERTAINTY IN LOCATING CONTOURS 


HORIZONTAL SCALE: APPROX. 1:18,000, 
FEET 2000 1000 ° 2000 2000, 


4000 FEET 


- =] 
METRES 500 250 o 300 1000 1300 METRES 


OLaKE 


AS Ma 
Mac, 
Sey2ie —HIGHWay 


\ 


= 1'035— 


Figure 18 


Surface contouring, Mile 360-370, Wrigley area (southern 
site) 1973 photography. 
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Figure 19 


Surface contouring, Mile 280, Camsell Bend area, 
1973 photography. 
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ABS TRACT 


Physical characteristics of snowbeds in the Richardson Mountains 
are investigated. The study conducted in the vicinity of Canoe 
Lake (23 miles west of Aklavik) was carried out in the field 
season lasting from May 14 - August 20, 1973. It examines the 
impact of late snowbeds upon biological elements and other 
physical processes in and adjacent to snowbeds. Role of other 
factors affecting active layer development, ground temperature 
regime and slope stability, in association with snowbeds is 


discussed. 
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Introduction 


This study was undertaken to examine the physical 
characteristics of snowbeds in the northern Richardson Mountains, 
Northwest Territories, Canada (Map A). Though snowbeds do vary 
in size, in most cases snow at the beginning of the melt season 
ranges from one to six meters in depth. The last of the snow melts 
sometime in the latter part of July or in August, or during colder 
summers, can last for more than one year. For the most part, though, 
snowbeds in the northern Richardson Mountains are annual or biennial 
features. This study was conducted in the vicinity of Canoe Lake, 
N.W.T. (68 ,04'N,135°,30'W),23 air miles west of Aklavik, N.W.T. One 
snowbed, referred to as Snowbed *l, was studied intensively during a 
field season lasting from May 14 to August 20, 1973. Several other 
snowbeds in the immediate area were kept under observation throughout 
the summer, and, as time and access permitted, surveys of major 
snowbeds in the region were carried out over a period of four months. 

Work performed during the 1973 season was the continuation 
of a three year project on the impact of late snowbeds upon associated 
biological elements and upon physical processes in and adjacent to 
the snowbeds. Other than their importance in the life cycles of many 
small mammal species and insects in the mountains, as well as containing 
a very specialized flora, snowbeds can be important as a source of water 
during the summer for smaller streams in the Richardson Mountains. 
Moisture from melting snow can also have a profound effect on the stability 
and active layer development of slopes under and near the deeper snow 


accumulations. 


Materials and Methods 


During the fall of 1972, wooden poles 270 cm long and 2.5 cm 
in diameter were placed around three snowbeds in order to mark the 
borders of the snowbeds and to determine snow depth the following spring. 
Snow depth throughout the 1973 season was determined by measurement along 
these poles or by probing with a 3/8 inch diameter, 2 meter sectioned, 


aluminum pole. Snow melt was determined by either the drop in actual snow 
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depth as measured by probing or by melt along one inch diameter, 

white wooden poles placed at 5 or 10 meter intervals in three profile 

lines perpendicular to the length of the snowbed. Due to some dirricul bres 
with thick ice layers within the snowpack, (snow depth and melt), were also 
checked by digging snow pits, generally to ground level. Profiles of snow 
type, such as location of ice layers, were measured in these snow pits and 
the density of the snow in each major snow or ice layer was determined using 
aluminum cylinders 20 and 48 cm long and 4.7 cm wide. A beam balance was 
used to weigh the samples on location. Densities of the surface layer of 
snow were also measured throughout the summer. The melt rate of Point CD75, 
at the center of the snowbed, was used as the standard rate, as CD75 represented 
both average conditions on the snowbed and the longest readings for snow melt. 
To record horizontal retreat of the snow, the snow border was marked with 


coloured flagging tape at 3 to 7 day intervals. 


Albedo and incoming shortwave radiation were measured along the 
three profile lines with a portable MK6 Sol-A-Meter, produced by Matrix, 
Inc., Phoenix, Arizona. This meter, with a spectral response between 
0.35 and 1.15 microns, measures incident solar energy in BTU per hour per 
square foot, with a +5% accuracy, and albedo, in percentage of incident 
radiation, with a +2% accuracy. All values were taken perpendicular to 
the immediate slope, and those of incident solar radiation converted to 
calories per cm? per minute. Incoming shortwave radiation was also 
measured with a recording bimetallic actinograph, made by C.F. Casella & 
Co., Ltd., London. The actinograph was placed at the upper edge of 
Snowbed *1 (location point CD 20), on the same SE facing slope as Snowbed *1. 
Values, taken at 15 minute intervals, were converted to cal/cm2/15 minutes 
or 24 hours, by a computer program written by M. Smith, Geography Department, 
Carleton University, Ottawa. 

Wind profiles throughout the summer for Snowbed *1 were measured 
with a hand-held anemometer, made by C.F. Casella & Co., London, model 
204026, in meters per second or kilometers per hour. Measurements were 
taken at 1 meter above the ground and averaged over 2 minutes. Wind 


direction was determined with a compass. Total wind runs, in kilometers 
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LOCATION OF INSTRUMENTS 


anemometer at the top of Ridge 1 


upper anemometer, upper rain gauge, actinograph, 
hygrothermograph 


lower anemometer, lower rain gauge 


approximate May 15, 1973 snowborder, snowbed 1 


see ARM 


over recorded time periods, were measured with three cup counter 
anemometers placed at 1.5 meters above the ground. The anemometers, model 
W5002) madewby~ Gaselilas® Co., had an accuracy of + 1.6 kilometers per 

hour. The three anemometers were located on a transect down Ridge 1 

(see Map B), with the anemometer on the top of Ridge 1 at 1550 feet 
elevation, the upper anemometer at 1350 feet (just before the top edge 

of the gulley containing the main part of Snowbed *1), and the lower 
anemometer at 1300 feet (in the band of lower shrub willows, but not 
sheltered by adjacent willows during winds from a northerly direction). 

Local air temperature and relative humidity were determined 
with a recording hygrothermograph, model E9389, produced by R. 

Fuess, Berlin, placed on a stand 30 cm above the ground, on the 
upper edge of the snowbed. For determination of daily averages, 
chart values were taken every half hour. 

Precipitation was collected in two standard 5 inch diameter 
rain gauges, model W5002, Casella & Co., Ltd., London, set at 30 cm 
above the ground to avoid splash-up from surrounding vegetation. A 
standard snow gauge was installed in September, 1972, but unfortunately, 
it was knocked over by a snowmobile during the winter and no record 
obtained. 

Readings for total wind, precipitation, snow melt and snow depths 
were taken as early as possible in the morning. Values for daily average 
temperature and relative humidity, maximum and minimum temperature and 
relative humidity, and total shortwave radiation were taken from midnight 
to midnight. Profiles of albedo, incoming shortwave radiation, wind and 
snow density were taken as close as possible to solar noon. These profiles 
were run along three lines, perpendicular to the length of the snowbed, 
about 40 meters apart, starting with the upper May 15, LO 722and e973. 
snow border as the zero point. The East Doweling (ED), Central Doweling 
(CD), and West Doweling (WD) lines were respectively on the east end, 
central portion, and west end of the snowbed. All distances along these 
lines were measured with a 30 meter metal tape. 

Soil temperatures beneath the snow bed were sampled by 66 thermis- 


tors with a direct temperature read-out, non-recording meter burl. by the 
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Carleton University Science Workshop. The thermistors had been 
inste@bled in August) 1972, at depthsvot (0, S710, 2155530, oreo0acm 
in the ground. Those at "0"' cm were placed either next to the soil surface 
under the immediate vegetation cover over the ground, or under the 
first leaves of the litter layer. Those thermistors at lower depths 
were inserted into the ground with a 4 inch aluminum probe, thereby 
minimizing disturbance of the ground. At the end of September, 1972, the 
thermistor wires were elevated on poles 1.5 to 2 meters above the 
ground. Temperature readings taken in November were done for those 
thermistors which could be located on their poles, and, in the 1973 
season, readings started once the wires and tops of the poles had 
melted out or were dug out. Unfortunately, the weight of the winter 
snow broke some of the poles, so that the wires became deeply buried. 
Readings, in these cases, were not obtainable until much later in 
the season. For most of the thermistors, though, temperature values with 
a snow cover still present were recorded. Other gaps in the readings 
occurred when the author was in other field locations. During the 
1973 season, readings were usually taken in the morning, at noon, 
and as late as possible in the evening (generally about 6 pm). At 
the end of the summer, the thermistors were recalibrated in an ice 
bath. 

Active layer depths throughout the summer were measured with a 
% inch solid aluminum probe along profile lines, or by the actual depth 
of soil pits. The presence of large chunks of rock or of bedrock in 
areas adjacent to the snowbed made active layer determination difficult 


there. 


RESULTS 
Distribution of Snowbeds 


As seen by Pictures FAl and FA2 in the Appendix, snowbeds are 
by no means unusual features in the northern Richardson Mountains. 


Though many slopes are snow-covered at the beginning of the melt season, 
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SNOWBED LOCATION 


MAP C 


July 1, 1973 


location of snowbed, 
- name of snowbed 
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this snow cover is not long lasting, melting by the beginning to 
middle of June. The deeper accumulations of snow are found in 
depressions greater than one meter, in creek gulleys of any orien- 
tation (though the drifts are deeper on southerly facing banks, or, 
in the case of wide stream valleys, are confined along southerly to 
easterly facing banks), and on southerly to easterly (usually SE) 
facing slopes. Map C shows location of snowbeds, chiefly those that 
were watched more closely over the summer. All those labelled on 

Map C still had snow left in the second week in July, but only Snowbed 
A had snow left at the end of the summer season. The melt of the 
majority of Canoe Lake snowbeds in 19/73 was three to four weeks ahead 


Of the melt im 0971: and 972. 


Snowbeds are not features unique to the Richardson Mountains 
but are found throughout the western Arctic. As Pictures Fl to F4 
show, the Shingle Point, Yukon region has many such larger snowbeds. 
During the middle of July, 1973, these snowbeds were much larger and 
thicker than those of the Canoe Lake area. Some, such as in Picture F3, 
had as much as one meter of ice on the bottom, suggesting that these 
accumulations of snow had lasted longer than one year. In an aerial 
survey over Shingle Point on September 1, there still were large patches 
of snow along the coast and creek valleys. Location of some snowbeds in the 


Shingle Point area is indicated in Map D. 
Profile Lines and Areas 


Pictures Al to A4 show the snowbed most intensively studied 
during the summer. Though three profile lines, used for various data 
determinations, were run across the snowbed perpendicular to the length 
of the snowbed and the contours of the slope, data used in this report 
have come mainly from the center line (CD line - see Picture A3), which 
ran in a straight line, heading 32° west of true north, from the May 15: 
1972 and 1973 snowborder (Point CDO) to the easterly gentle slope of the 
next ridge. The following outline gives CD points, taken at 5 meter 


intervals along the slopes, giving slope angle, area location, and major 
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River 117A sheet, Scale 1:250,000) 


30 


type of vegetation. Plant specimens were identified by the author 
and verified by Dr. J. Lambert, Biology Department, Carleton University, 
and by Bill Cody, Plant Research Institute, Ottawa. Some of the plant 
communities have been previously described by Lambert (1968). For 
correct plant names, see Table 7. (Appendix) 
CD -20 to 0 - Picture A2 
- slope angle 5 to 10°, direction South 


- above the May 15, 1972 and 1973 snow border, 
on Ridge 1 (see Map C). Winter snow cover 
0 to 5 cm over tops of solifluction features, 
up to 60 cm in depressions. Most areas without 
winter snow cover. 


- landscape dominated by elongated solifluction 
features (lobes), lichen-moss communities inter- 
spaced with bare mud and rock patches, and 
hollows with 10 to 30 cm high shrubs and a cover 
of moss. Soil very rocky. Average growth height 
on lobe tops 1 to 5 cm. 


CD 0 - Picture A2 
- slope angle 59, slope direction S 


- May 15, 1972 and 1973 snowborder. 50 cm drop 
in elevation from the slope above. 


CDMO-=235 - Picture A2 
- slope angle 0 to 10°, slope direction S 


- early melt area (May 17-27, 1973). Snowpatch 
with 5 to 35 cm of snow. Always covered with 
shallow snow in the winter, but the snowis early in 
melting. Moderately drained - no standing water. 


- dominated by wet mossy hollows, patches of lichen 
(Cladonia’ sp.) , scattered 107 to 30 cme hieh bircr 
(Betula glandulosa) and willow shrub (Salix pulchraje 
occasional tussocks of cotton grass (Eriophorum 
vaginatum L.), patches of Cassiope tetragona and 
other heath plants (Empetrum nigrum, Vaccinium 
uliginosum, Vaccinium vitis-idaea, and Ledum 
palustre), plus other assorted vasculars including 


Arctostaphylos alpina, Carex sp., Petasites frigidus, 


Equisetum arvense, and Arctagrostis latifolia. 
Average growth height 10 to 20 cm. 


CDs 2 — 30 - location of actinograph, upper anemometer, upper 
rain gauge, and hygrothermograph 
CD 40 - 45 - Picture Al 


- slope angle 10 to 15°, slope direction SE 


- melt area May 30 - June 6, 1973. Start of snow 
over one meter in depth. Upper edge of Snowbed *1. 
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- Drop in elevation of at least one meter 
from CD 35. Well drained. 


- band of 50 to 100 cm high shrub willows 
(Salix pulchra) of almost continuous cover. 
Leaf litter over 80% of the ground surface in 
this band. 


CD. 45 - Picture Al 
- slope angle 10°,slope direction SE 
- melt area June 6, 1973. Well drained. 


- edge of upper willows and upper meadow (see CD 50) 


CD50! "95 - Pictures Al and A3 
- slope angle 20°, slope direction SE 


= melt area June 1/-27, 1973. Soil with only 
scattered large rocks. Locally, the slope up 
to 45°. Well drained. 


- upper meadow vegetation generally 20 to 40 cm 
high, ‘though itn patches up to. 60 cm in height. 
Usually an 80 to 100Z cover of moss (especially 
Polytrichum sp.) with a lush, thick growth of 
grasses, sedges, horsetail and flowering herbaceous 
species. Some areas with almost a 100% cover of 
ground willow (Salix Chamissonis) 2 to 5 cm in 
height. Meadow species including Petasites 
frigidus, Arctagrostis latifolia, Equisetum arvense, 


Saxifraga punctata, Carex podocarpa, Artemesia 
arctica, Artemesia Tilesii, Arnica lessingii, 


Antennaria a monocephala, Achillea Ilea borealis and 
Polygonum | bistorta in this area. 


CD 60 - Pictures Al and A3 
- slope angle 10°, slope direction SE 
= melt area June 27 — July 2, 1973. Well drained. 


- 40 to 70% cover of moss (especially Polytrichum sp.). 
Vegetation 20 to 30 cm in height, not as thick as CD 
45-55, with 30% cover provided by the ground willow, 
Salix polaris subsp. pseudopolaris, Ito: 2 cm in bei ght. 
Other vascular plants in this area including Carex 
podocarpa (30% cover), Saxifraga punctata, 
Arctagrostis latifolia, Equisetum arvense, Oxyria 
digyna, and Myosotis lyosotis alpestris. 


CD 65 - Pictures Al and A3 
=| «slope angle 12°, slope direction 55 


- melt area July 5-8, 1973 (late melt area). Well 
drained. 


Selo. 
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CD75 


CD 80 


CD. 65 


30 to 80Z of the area bare mud and rock. Vegetation 
sparse, low growth (average 5 to 10 cm), but species 
composition like CD 60. 


Pictures Al and A3 
slope angle 1. slope direction SE 


melt area July 8-12, 1973 (Tate melt area). 
Well drained. Bottom edge of SE facing slope 
from CD 40 to CD 70. 


102% bare mud and rock with a 40 to 70Z cover of 


moss (especially Bryum sp. and Drepanocladus sp.). 
Vascular plants sparse, with Oxyria digyna (very 
common) , Poa arctica, Equisetum arvense, and 


Ranunculus p pygmaeus. 


Pictures Al and A3 
slope angle sae slope direction S 


melt area July 12, 1973 (late melt area). Moderately 
drained. Flat bottom area of snowbed. 


patches of bare mud and rock over 30% of the area. 
Moss cover (especially Drepanocladus sp.) over 70% 

of the: area.,. Plants 2 to 20 cm in’ heieht, with 

such species as Arctagrostis latifolia, Poa arctica, 
Carex Lachenali, Eriophorum Scheuchzeri, Ranunculus 
pygmaeus, Saxifraga punctata, and Salix polaris subsp. 
pseudopolaris. 


Pietures Al and A3 
slope angle shat Slope direction §S 


melt area July 8, 1973 (late melt area). Locally 
well drained. 


100% cover of moss (especially Drepanocladus sp.) 
with ground willow (Salix polaris subsp. pseudo- 


polaris), Poa arctica, Oxyria digyna, Arctagrostis 
Jatifolia, Ranunculus pygmaeus, Petasites frigidus, 


Saxifraga punctata, unctata, and Equisetum arvense. Growth 
2 to,.20 cm iniheitent. 


Pictures Al and A3 
slope angle isa slope direction §S 


nelt area, Iulyods-<= 8, 1973 (ate melt area). 
Poorly drained. 


lower wet meadow with 100% moss cover (especially 


Drepanocladus sp. and Bryum sp.). Vascular plants 


including Ranunculus nivalis, Carex Lachenali, 
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Saxifraga punctata, Arctagrostis Fatifolia, 
Equisetum arvense, Salix polaris subsp. pseudo- 


polaris, and Oxyria digyna. Growth generally 
10 to 20 cm high. Occasional dead shrub willows 


(Salix pulchra). 

Pictures Al and A3 

slope angle as slope direction § 

Melt area June 27-July 2,°1973. Poorly drained. 
lower wet meadow. Vegetation 20 to 50 cm high 
(see CD 85). Scattered dead Salix pulchra bushes. 
Pictures Al and A4 

slope angle a, slope direction §S 


melt area June 14, 1973. Poorly drained — 
ground very soggy. 


lower wet meadow. Vegetation 10 to 20 cm high 
(see CD 85). Scattered live 10 to 20 cm high 


young Salix pulchra. 

Pictures Al and A4 

slope angle On 

melt area June 10-June 14, 1973. Poorly drained. 


main band of lower willows (Salix pulchra)- 

average height one meter,though up to two meters 

in height. 100Z moss cover (especially 
Drepanocladus sp.) with scattered elongated 
vasculars, including Equisetum arevense, 

Artemesia Tilesii, Polemonium acutiflorum, Carex 
Lachenali, Arctagrostis latifolia, and occasionally 
Salix Chamissonis. 


location of lower anemometer and lower rain gauge 


in a patch of open meadow in the lower willow band 


Pictures Al and A4 
slope angle oe 
melt area June 10-June 14, 1973. Poorly drained. 


lower willows (50 to 80 cm high) mixed with a lush 
meadow growth (20 to 35 cm high) composed of 
Artemesta Tilesii,, Poa arctica, Mpilobiumvaneust= 
ifolium, Saxifraga punctata, Equisetum arvense, and 
occasional Salix Chamissonis. 
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CD 105-110 - Ground collapsed this spring between CD105 
and CD110 due to the underground flow of water 
during spring runoff (see diagram). 


Gp 210 - Pictures Al and A4 
- slope angle 0° 
= melt area June 17, 19/3. Well drained. 


- lush, 40 to 60 cm high Calamagrostis canadensis, 
Senitiss with scattered plants of Carex podocarpa, 
Polemonium acutiflorum, Equisetum arv arvense, Salix 
Chamissonis, and Viola epipsila. 


CDMS - Pictures Al and A4 


- slope angle Or 


- bottom of creek bed at base of NW facing slope. 
Large pieces of rock over most of the area, in 
places bedrock (sandstone) exposed at the surface. 
Active layer "depths" taken on adjacent bank. 


CD20 - Pictures Al and A4 


- slope angle 45°, Slope direction 12° west of 
true north (NW). 


- melt area June 10-June 14, 1973. Well drained. 
Five meters above the creek on the NW facing slope. 


- Vegetation mainly shrub growth with Salix pulchra 
(50 to 70 em high), birch (Betula glandulosa, 20 cm higl 


and Vaccinium uliginosum (20 cm high). 
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CD 125 = Pictures -Aleand AG 
- slope angle 45°, slope direction NW 
- melt area May 27, 1973. Well drained. 


- 100% moss cover (especially Polytrichum sp., 
Dicranum sp., and the Iichen, Peltigera sp.) 
Vegetation damp, but well drained, 10 to 15 cm 
high with heath type species and shrub willow 
and birch. 
CD 7133 =, Pictures: Al and (A4 
- slope angle ees slope direction NW 
- May 17, 1973 snow border (winter maximum snow 
border) 
CD 125-140 - Pictures Al and A4 
- slope angle 40 to 45°, slope direction NW 


- well drained, vegetation 1 to 10 cm high, 
mostly lichen, moss and woody tundra species. 


CD 140 - top of NW facing slope. Ridge 2. 
- area of scattered mud patches. Vegetation lichen 
dominated. 
CD 140-160 - Pictures Al and A4 


- slope angle 2 to ee slope direction NNE to NNW 


- on slope of next ridge (Ridge 2 - see Map C). 
Main ridge slopes eastward, though these points 
locally slope gently in a northerly direction. 
Well drained, though hollows damp. 


- area of Eriophorum vaginatum Le tussocks and 


hummocks dominated by sedges and grasses. 80 to 
100Z ground cover of moss and lichen, though some 
patches of bare mud on hummock tops. Hollows 
between hummocks 10 to 50 cm deep with a 100% 
cover of moss and with 10 to 20 cm high birch 
(Betula glandulosa) and willow (Salix pulchra) 
shrubs. 


The snowbed was further divided into representative areas, chosen 
on the basis of vegetation type, date of exposure, depth of snow, 
drainage, and location in the snowbed. Some of these factors are outlined 
in Tables3A and 3B. The rest will be described here. Albedo values are 


usually the average of 22 area values, taken over snow-free ground in August. 
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Upper Lobe 


Picture Al 
slope angle 0 to phe slope direction $ 


10 to 30 meters above the May 15, 1972 and 1973 
snowborder. Snow largely confined to hollows 

5 to 30 cm (maximum 50 cm) deep, though snow 

cover over the tops of raised solifluction features 
can be up to 15 cm during some parts of the 

winter. Most lobe tops bare of snow throughout 

the year, though. Well drained on the tops of 
lobes, hollows of variable drainage (the vegetation 
in hollows is usually wet, with the surface 1 to 

3 cm dry occasionally). Albedo of lobe tops ranges 
from 21.5 to 232, withean average of 222. Ihe 
albedo of hollows ranges between 21 and 232%, with 
an average of 21.972. 


landscape dominated by flat, elongated solifluction 
features running downslope in a northern to southern 
orientation. Vegetation on lobe tops generally 
lichen dominated (Cladonia sp.,40 to 602 cover), 
with a 10 to 20% moss cover (especially Dicranum 

and Polytrichum sp.) and woody vascular plants such 
as Vaccinium uliginosum, Ledum palustre and Salix 
artica Pall. subsp. artica. Other less common 
vasculars include Arctostaphylos alpina, Pedicularis 
Kanei, Anemone narcissiflora, Lupinus arcticus and 
Vaccinium vitis-idaea. Vegetation usually dry and 
well drained. 


most hollows have a 100% moss cover (especially 
Dicranum sp. and Drepanocladus sp.). Wide variation 
in shading is proportional to the amount of shrubby 
plants in the hollow, such as Betula glandulosa and 
Salix pulchra. Some hollows deep, with 70 to 902% 
shrub cover over moss (shrub cover including Betula 
glandulosa, Salix pulchra, Vaccinium vitis-idaea, 
Ledum palustre, and Vaccinium uliginosum). Other 
hollows 5 to 10 cm deep, vegetated mainly by moss 
and Vaccinium uliginosum. Plants tend to reach the 
tops of hollows or extend to a maximum of 10 cm 
above the level of the hollow (see diagram). 
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Upper Birch 


Upper Willows 


Picture: AL 


slope angle 10 to neve. Slope direction aor 
east of south 


from above the winter snow border in some years 
(1972 covered, 1973 exposed) to a melt area of 
May 26, 1973. Starting where the slope of Ridge 1 
suddenly drops in elevation into the gulley 
containing Snowbed *l1. Well drained. Albedo from 
11 to 23%, depending on the ground cover under- 
neath the shrub layer. Average albedo 18.6%. 
Temperature of the ground under the shrub layer 
is dependent upon the depth of leaf litter, amount 
of moss cover, and amount of overhead shading by 
the birch. 


intermittent border 3 to 5 meters thich of shrub 
birch (Betula glandulosa) along the upper edge 

Or ionowbed “1, *Birch height 10 to 70 cm, 

average 41.5 cm. Birch cover usually complete, 
with a ground cover of dry, loosely packed moss 
and dead leaves, with scattered plants such as 
Vaccinium vitis-idaea (in some areas up to 80% 
ground cover), Vaccinium uliginosum, Empetrum 
nigrum, Pyrola grandiflora Radius, and occasional 


grass and sedge leaves. 


Pictures Al and A3 
slope angle mes slope direction SE to ESE 


melt area May 27 to June 7, 1973. Can be with or 
without a birch shrub border immediately above. 

Well drained. Albedo from 12 to 23%, depending on 
the ground cover under the shrub layer. Average 
albedo 18.8%. Temperature of surface proportional 

to exposure (nearness to border of the willows), type 
of ground cover and amount of overhead shading. 


band 3 to 5 meters wide of shrub willows (Salix 
pulchra) along the upper edge of the snowbed (level 
with CD 40-45), Willow height from 30°to 150 cm, 
average 79 cm. Willows not continuous-broken by 
open patches with meadow or heath vegetation. Under 
the willows proper, loosely packed leaf litter is 
common, with much dead wood and large, live, 
horizontally lying stems and trunks. Patches of 
loosely packed moss (especially Drepanocladus) 

also common. Undergrowth, for the most part, is 
sparse and spindly, with such species as Arctagrostis 
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Upper Meadow 


Upper Late Melt 


latifolia; Artemesia erctica, wtellarta 
longipes, Vaccinium uliginosum, EOawartica, 
Petasites frigidus, Viola epipsila la and Achillea 
borealis. 


Pictures Al, A3 and A4 


slope angle 20 slope direction SE to locally 
E. ‘or S facing 


melt area June 7-June 27, 1973. Well drained once 
the snow border passed the area. Albedo from 24 to 
35k, WLEn an averace Of 20,97. 


lush, mostly herbaceous flowering plant growth. 

On level with CD 45-60. 0 to 15 meters from the 
upper willow band. Most values used for this 
report were recorded at the level of CD52, 9 meters 
from the upper willows, in a melt area of June 10 - 
June 17, 1973. The vegetation here is dominated 

by Salix Chamissonis (up to 90% cover), which has 
an overall height of 2 to 10 cm Other meadow areas, 
with overall growth heights of 5 to 30 cm contain 
variable amounts of Salix Chamissonis along with 
the more common herbaceous species of Achillea 
borealis, Dodocatheon frigidus, Equisetum arvense, 
Carex podocarpa, Artemesia arctica, Aconitum 


delphinifolium, Arnica lessingii, Anemone Richardsonii, 
Anemone narcissiflora, Polygonum bistorta, Polygonum 


viviparun, m, Arctagrostis latifolia, , Castilleja Raupii, 
and Pedicularis sudetica. The Sanne of moss cover 
is variable, with pockets of 70 to 100% cover 


(especially Polytrichum sp. and Drepanocladus sp.) 


over about 50% of the meadow area. 


Pictures Al, A3 and A4 
slope angle 15°, slope direction SE. tos 


melt area June 2/-July 11, 1973. Level with CD 
65-70. Well drained once the snow border has 
retreated past this point. Local drainage pathways 
usually filled with flowing water, though. Most 
areas with a hard, tannish surface mud, wet below 

2 to 5 cm. Albedo from 17.5 to 24% with an average 
Of 20542. 


band of sparsely vegetated mud and rock 20 meters 
from the upper willows. 40 to 100% of the surface 
hard mud, pebbles, and/or sandstone boulders up to 
50 cm by 30 cm by 20 cm in volume. Vegetation 
present (growth from 2 to 4 cm, maximum 20 cm) 
includes clumps of moss (especially Drepanocladus sp. 
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and Bryum sp.), and scattered plants of 
Arctagrostis latifolia, Ranunculus nivalis, 


Carex podocarpa, Arnica lessingii, Oxyria — 
digyna, Ranunculus pygm: pygmaeus, Equisetum arvense, 


and Saxifraga punc punctata. Most growth is below 
10 cm in height and scattered between the rock 


and mud patches. 


Lower Late Melt 


Pictures Al; A3vand A4 
- slope angle 0 to ae slope direction § 


= melt area July Sto Jast melt inv July 21973. 
Poorly drained, ground usually saturated to the 
surface. Flowing water in 1 to 3 cm deep drainage 
pathways. Debris from melted snow common on the 
surface of the exposed ground. 


- lower flat slope of the snowbed. Next to and 
extending 10 meters south from the SE facing upper 
slope. On level with CD 75-80. Patches of bare 
mud and rock over 30 to 100% of the surface in 
areas. Moss cover important - usually 70% cover 
(especially Drepanocladus). Plants 2 to 20 cm 
high. Sparse. Common species include Oxyria 
digyna, Ranunculus pygmaeus, Carex Lachenalii, Poa 


arctica, Equisetum arvense, Arctagrostis (aici 
Salix olerie subsp. pseudopolaris, Petasites 
frigidus, Saxifraga punctata, and occasionally 
Eriophorum scheuchzeri. 


Pictures Al, A3 and A4 


Lower Meadow 
- slope angle 0 to a slope direction S 


= melt area June /-June 30, 1973. Poorly drained, 
ground usually saturated to within a few cm of 
the surface. Standing water in places. Albedo 
from 21 to 29%, with an average of 26.82. 


- wet meadow on the lower side of the snowbed. 0 to 
8 meters from the lower willow band. Vegetation 
not high, 5 to 20 cm, with a 100% ground cover of 
moss (especially Drepanocladus sp. and Bryum sp). 


Common vascular plants include Carex podocarpa, 
Carex lachenali, Carex aquilinus, Arctagrostis 
latifolia, Salix polaris subsp. pseudopolaris, 
Equisetum arvense, Petasites frigidus, Saxifraga 
punctata, Ranunculus pygmaeu: pygmaeus, and Poa arctica. 
Sedges (Carex sp.) usually form the largest 
proportion on of the plant cover. 
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Pictures Al and A 3 


Lower Salix Cham. 
- slope angle he slope direction S 


- melt area June 14, 1973. Locally better drained 
than most of the lower snowbed. Saturated during 
times of high runoff from the snowbed or slopes 
above. No standing water. Albedo 27 to 29.52, 
average 28.32. 


- on the lower side of the snowbed; better drained 
area with a 60 to 80% ground cover of the ground 
willow, Salix Chamissonis, common also in upper 
meadow areas. Height usually 5 cm. Scattered 
individuals of the same vascular plants as the 
lower meadow (growth up to 30 cm) mixed through 
the Salix Chamissonis, with occasional young 
shrub willow stems 10 cm high (Salix pulchra), 
plus Vaccinium uliginosum (5 cm), Artemesia 
Tilesii (10 cm), and Pedicularis sudetica leaves 
(38cm). 


Lower Willows = Pictures Al,,A3 and A4 
- slope angle 0 to yaa slope direction S 


- melt area June 7-14, 1973, though as late as 
June 30 in localized, sheltered patches. Ground 
usually saturated within 10 to 20 cm of the surface, 
and often saturated within 1 to 5 cm of the surface. 
Albedo from 10 to 19%, average 15.1%. Light reaching 
the ground cover under the willows patchy due to 
shading by the willows. 


- band 10 to 20 meters thick of 1 to 2 meter high 
shrub willow (Salix pulchra). Located south and 
below the SE facing, steep slope of the snowbed 
and the lower wet meadow, but north of the creek 
at the bottom of the snowbed and the steep, NW 
facing creek bank. Includes CD 100-105. Ground 
cover, under the almost continuous canopy of 
willows, generally 100% cover of moss (especially 
Drepanocladus sp.), with elongated vegetative 
shoots of such green vascular plants as Viola 
epipsila, Arctagrostis latifolia, Petasites frigidus, 
Equisetum arvense, Artemesia Tilesii, and 
Polemonium acuteflorum. Leaf litter is sparse. 


Formation and Ablation of Snowbeds 


Formation 


As seen from Tables 3 and 4, (Appendix), the total snow depth in 


this region over a winter is unevenly distributed, even early in the winter, 
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due to drifting by the wind. Over the tops of ridges, snow deposition 
is minimal, though local depressions can be filled to a depth of one 
meter. On the sides of ridges and in flat bottom lands, the snow is 
generally 0 to 25 cm deep. Slight depressions or drops in topography, 
hollows between solifluction features and vegetation tussocks, and most 
mudslide scars have average snow depths of 25 to 70 cm. These patches 
of snow are penetrable by light when the snow is dry (Sellers, 1968) 
and generally melt early in the season (from the middle of May to the 
first week in June ). Snowbeds, on the other hand, with depths of snow 
in excess of one meter, melt out from the middle of June to the end of 
the summer season. The larger snowbeds usually ablate by the end 

of August. However, in some years they may not melt completely. 
Snowbeds are confined to the deeper depressions, creek gulleys, and 
sudden steep drops in topography along southerly facing slopes. 

The occurrence of larger drifts on SE facing slopes is an expression of 
the prevailing northwest wind in the winter. 

Snowbeds in the Richardson Mountains occur year after year in 
the same location, as manifested by airphotos (1954, series A14133-93, 
Air Photo Division, Department of Energy, Mines, and Resources, Canada), 
personal communication with Dr. J. Lambert, who also worked in this area 
in the mid-1960's, the author's three summers in the field in the Canoe 
Lake area (1971-1973), and the distinct vegetation patterns associated 
with the snowbeds. From the author's experience, though, the snowbeds 
did vary in depth and in the exact location of the face of the drift 
each year, but covered approximately the same area each season, with a similar 
proportion of snow distribution throughout the snowbed. Time of snow 
disappearance did vary considerably, though, with a three to four week 
range over the three summers. 

The characteristics of snowbeds are influenced by inter- 
relationshipSamong topography, wind and amount of snowfall. Table 2 
(Appendix) illustrates the relationship between these variables. The 
location of the anemometers in Table 2 is shown on Map B. As seen in 
this table, the drop of wind down the slope of Ridge 1 is consistent, 


and is especially noticeable during storms from the northwest 


Soe 


or west (such as August 13-the wind varied from 20 km/hr average 

at the top of the ridge, to 14 km/hr just above the snowbed with 

a drop of about 45 meters in elevation, to 4 km/hr at the bottom 

of the snowbed with a drop of about 15 meters in elevation). The 
August 1 wind profile, Graph 1E, better illustrates this wind drop 
down the slope. When the wind passes over the sudden drop in topo- 
graphy from CD35 to CD50 (a gradient of 5° at CD35 to a gradient 

of 45° at CD50, with a change of 2.5 meters in elevation), wind 

speed decreases drastically, with the effect that, if the wind is 
carrying snow, much of its load is dropped, and a drift develops 
downslope from the drop in topography. The wind picks up again 

along the lower slopes, but is reduced in the bank of one to two 

meter high willows at the bottom of the snowbed area (CD95 to COLES). 
again dropping snow. Once going up the northwest facing slope along 
the 3 to 5° easterly facing slope of the next ridge (CD120-160), the 
wind rapidly picks up speed, blowing snow away from these areas. This 
pattern of snow deposition is shown in the snow depth recordings 
listed for Tables 3 and 4. The snow depths above the snowbed on 
September 19, 1972, were negligible on the surface of solifluction 
features, though in the birch filled depressions, snow was up to 60 cm 
deep. Snow patch areas (such as CDO-35) had an intermittent cover of 
up to 20 cm. With the sudden drop in topography at the upper willows 
(CD40-45), the snow depth rose to 120 cm, with a drift extending about 
five meters into the upper meadow. Below the drift, the snow depth 
over the upper and lower slopes ranged between 10 and 20 cm, up to 

the area of the lower willows. Snow depths listed for the lower 
willows (10 to 20 cm) are misleading, as much of the deposited snow 
actually rested on the branches of the willows. The northwest facing 
slope was clear of snow, but hollows in the hummock area above were 
filled to a depth of 30 cm. By November 24, 1973, the snow deposition 
pattern had changed. Areas immediately above the snowbed were covered 
with snow from 0 to 15 cm. Snow in the region of the upper willows 
had been redeposited further downslope, decreasing the snow depth here 
from 120 to 75 cm. The main drift of Snowbed *1, now much further 


extended, covered the upper slopes of the meadow with depths around 140 cm. 
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Even the lower flat meadow areas had depths from 60 to 160 cm, 

with the greater depths being recorded in areas where the upper 
slope ended and the lower flat slope began (160 cm of snow). The 
lower willows, with a snow ground cover of 75 ecm, still had much 
snow caught up on branches and leaves. As can be seen by these 
data, the snowbed changed shape from a steep faced drift at the 

top of the slope in September, to a much flatter, deeper drift 

by November. By May, 1973, though, as Tables 3 and 4 show, the areas 
above the snowbed were mostly free of snow, the shallow snow patch 
areas had snow from 8 to 33 cm, and the snowbed proper had measured 
depths from 84 to 510 cm. Picture C(Appendix) shows the lip.of 

the drift, which by this time had extended at a 10° angle as far as 
CD90-95. The actual drift itself had a face 2 to 24 meters high. 
Upper willow areas in early May were covered with snow over 2 meters 
deep. This deep drift, extending out at a 10 to 20° angle from 

the drop in topography, was deepest in the flat meadow area adjacent 
to the bottom of the SE facing slope (CD65-75). Once past the face 
of the drift, the snow covering the lower meadow and lower willows 
was 150 to 200 cm deep, with a deeper area at CD115 (the location 

of the creek bed). Snow extended up the northwest facing slope 
about 10 meters, but had been swept from the upper parts of the hill. 
Only hollows between hummocks on the slope behind were filled with 
snow, along with shallow drainage pathways. 

Other snowbeds examined in the region also showed this 
development of drifts along southerly facing slopes, due to sudden 
drops in topography. The other common mode of snowbed formation was 
the filling in of local, deep depressions that were narrow enough to 


be totally covered, such as creek gulleys. 


Ablation 


The factors affecting melt rates of snowbeds are many and inter- 
acting. It was found, for the particular snowbed studied, that incoming 
solar radiation was by far the most important factor determining melt rate. 


First, though, the other factors affecting ablation will be considered. 
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Wind 

Wind has an input into the melt rate of snow by controlling 
the rate of horizontal water vapour transfer in the air above the snow. 
Wind moving across snow surfaces will create a water vapour 
gradient above the surface proportional to the rate of removal 
of saturated air. Sublimation of snow and evaporation of water from 
melted snow is faster into air of lower water content. 

As seen from Graphs 1A to 1G, the movement of wind across 
the snowbed changed considerably through the season as the snow 
melted. On May 26 the snow cover was at a slope of 10 to 205 
compared to a slope of 10° for the area above the snowbed. Snow 
melt had only begun, with a drop of about 30 cm since May 14 in 
the central portion of the snowbed. Snow cover was continuous, 
except for CD30. The face of the main drift forming the snow cover 
was at CD90, with CD95 half way down the face of the drift. CD100- 
105 were located in the flat creek gulley which was still covered by 
115 and 163 cm of snow respectively. CD110 and 115 were unreachable - 
spring runoff in the creek had begun the night of May 25 and water 
was flowing over the snow in this area. Wind flow, from the WNW, 
across the snowbed (Graph 1A) did decrease slightly down the slope 
but its effect upon water vapour levels must have been relatively 
uniform across the slope and must have been of some significance 
since the wind only dropped slightly in speed down the snowbed. 
Table 2 (Appendix), also shows this effect; e.g. on May 21, the 
daily average wind speed on top of the ridge behind the snowbed was 
15 km/hr. At the start of the main part of the snowbed (CD25), it 
was 13.8 km/hr, and at the bottom of the slope (CD120), it was 11.5 km/hr. 
The decrease in wind velocity was less than a third across the snow- 
bed. On May 26, the averages were 9.5, 8.5, and 6.7 km/hr respectively. 
This even effect of the wind was noticed also with snow piled up at 
the sides of snowpits. This snow soon melted (in less than a week) 
to the level of the rest of the snowbed from a height of 40 to 70 cm 


above the general surface of the snowpack. 
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From May 18 to June 2, the snow depth thinned an average of 
of 80 cm. Though underlying topography had only just begun to 
be exposed, the actual snowbed itself was at a greater angle 
(slope of 20 to 25°) to the slopes above (slope of 1020 Snow 
started at CD45, meaning that from CD40-45, the upper band of 
shrub willows had melted out, with a drop in elevation of at least 
50 cm. Graph 1B shows a noticable drop in wind speed in this 
area. CD45-120 at this time were snow covered, with CD90 at the 
face of the drift and CD95 half way down the drift. CD100-115, 
in the flat creek region had snow depths of 53, 90, 139 and 265 cm 
respectively at this time. The edge of the snow with the NW facing 
wall reached CD120, leaving CD125-140 above the snowbed and CD140-155 
on the snow-free,flat hummock area above. There was a noticeable 
decrease in wind speed across the snowbed, with the largest drop in 
the region of the creek gulley. Speed rapidly picked up again along 
the NW facing wall and hummock areas above. During the wind storm of 
June 2, predominately from the NW, the wind on the NE facing slope on the 
other side of Ridge 1 was 25 km/hr, on the top of Ridge 1, 23 km/hr, 
and on the general slope above the snowbed averaged 13 to 15 km/hr. 
The rapid decrease in wind speed once the wind entered the snowbed 
area should mean that the wind would have a decreased effectiveness 
in removing water vapour from above the snow compared to earlier in 
the season. Further measurements done as the snow continued to melt 
back closer to the underlying topography show that the wind was 
increasingly affected by the steeper angles of the snowbed. By June 9, 
the immediate drop in wind speed, upon passing into the snowbed area 
at CD40 is even more apparent than June 2 (see Graph 1C). The snow, 
by this time, had dropped a further 44.5 cm since June 2, though the 
border was still at CD45. Further changes had occurred on the lower 
slopes, with the snow border having retreated to CD115, leaving CD120- 
140 on a snowfree NW facing slope. The snowbed at this time was at 
a 25° angle, compared to 10° for the slopes above. As can be seen 


from Graph 1C, the wind did pick up again slightly across the snowbed, 
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but dropped again going over the lip of the drift, and remained 
low across the flat creek gulley. It immediately increased 

in speed again going up the NW facing slope and over the 

hummock area to the south. Even in this condition of higher 

wind speeds (20 to 25 km/hr from the NW), the actual wind passing 
over the snow surface was a half or less of the wind over exposed 
areas. 

By July 5, the snow had retreated to CD65-85 with a 
considerable drop in volume. The angle with the slope above (now 
the slope of the upper meadow SE facing sidehill with a slope 
of 20°) was 20°. With this drop in snow, the upper willows (CD40) 
and upper meadow (CD45-60) were completely exposed, as well as 
the areas from CD90-160. The wind, from the NW, averaged 33 km/hr 
with gusts over 50 km/hr on the slopes above the snowbed. As 
Graph 1D shows, there was an immediate drop in wind speed past the 
upper willows, with the speed on the upper side (CD50) at 11 km/hr. 
The wind speed remained very low (5 to 7 km/hr) down the SE facing 
slope, picking up slightly in the lower meadow (13 km/hr maximum), 
but dropping again in the lower willows (4 to 5 km/hr.) Across 
the entire extent of the remaining snowpatch, the wind remained 
at speeds less than half that for the tundra above. Its effect 
on the snowbed was most likely minimal, though the drying out of 
upper tundra areas was readily apparent. 

Table 2 (Appendix) shows how the effectiveness of the wind 
over the snowbed decreased as the snow melted back, both horizontaliy 
and volumetrically. It also shows most major wind storms throughout 
the summer coming from a northwesterly direction (27 out of 41 days 
where the average wind speed was greater than 10 km/hr on top of 
Ridge 1). Even during winds from a southerly direction, this drop of 
wind speed across the snowbed was apparent (see Table 2, June 11, June 
13 and June 23). These values can be compared to data from earlier 
in the season (see Table 2, May 21, May 28 and June 2). The decline in 


speed is less than 50% over the snowbed while the snow was of 


considerable volume. 
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The changing correlation of wind and melt rate was not shown 
in melt rate data. The effect of high wind speeds tended to be 
minimal when they coincided with periods of low temperature, low 
solar radiation, and high relative humidity (Table 6, Appendix). 
During May and June, high winds did occur on sunny days, but the 
air temperature on these days was low. As seenin Graph 1G and 
Table 1, there is no readily apparent correlation between total 
wind run or daily average wind speed and daily average snow melt 
or total snow melt. As Graph 1G shows, solar radiation seemed 
to be a more determining factor in snow melt. During the time 
period June 22-25, the wind speed was high (average 14 km/hr) but 
the incoming radiation was low (340 eaten jaan) along with 
the melt rate (2.7 cm daily average). In contrast, during the 
period June 15-21, with a high average wind speed (14 km/hr) and 
a high average incoming radiation (600 elfen jase the melt 
rate was high (14.3 cm daily average). Of course, the process 
was accelerated by ambient temperature, albedo, and snow density. 
During the period May 18-22, radiation was high (560 ea lvena/ ie 
wind was moderate (6 km/hr), but temperature was low (3.6°C 
average), albedo was high (602%) and density was low (0.54 en 


The resultant daily average snow melt was 1.4 cm. 


Rain 


Rain can affect the melt rate of snow by providing an energy 
input in the form of unfrozen water into the snow system, forming 
ice layers in the snowpack and increasing the snow density. As seen 
by Table 1 and Graph 2A, periods of heavy rainfall did not correspond 
to periods of high melt rates. Days with large amounts of precipi- 
tation usually corresponded to days of high relative humidity, low 
temperature and low solar input. The actual energy input of rain 
was not significant in melting the snow, but the rain augmented the 
snowpack water equivalent. This water was stored for later release 
during snowmelt, especially during days of significant incoming solar 


radiation. Days of high solar radiation, as Graph 2A shows, tended 
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to correspond to days of lowest total rainfall. The direct 
relationship between amount of total rainfall and daily average 
melt rate appears to be a reverse one, because of the factors 
outlined above. One important effect of rainfall, though, is the 
movement of runoff waters from upslope areas, under the snowbed, 
next to the frozen ground, creating air spaces between the snow 
and ground which allow movement of warmer air under the snowpack. 
Some of the spaces created are as deep as 10 cm in localized areas. 

As Table 2 shows, there was no observable gradient of 
rainfall between the upper snowbed (CD30) and the lower snowbed 
(CDLIO)= 


Density 


Densities listed in Table 1 are averages of snowpit 
profiles. Actual surface densities, in areas with 10 to 15 cm of 
snow/ice left, were recorded as high as 0.90 on one Snow on 
top of ice layers next to the ground (this type of ice layer 
was usually 5 to 15 cm thick) was very wet, its density averaging 
0.80. As Table 1 shows, there was an overall increase in density 
throughout the summer, due mainly to physical settling of the 


ripened snow pack. 


Albedo 


Albedo, the percentage reflection of incoming short wave 
radiation, decreased over the summer (see Table 1), corresponding 
to the increasing daily average melt rate shown on Graph 4A. This 
decrease in albedo was due to the melting out of plant debris and 
soil particles that had been deposited in the snowpack along with 
the snow during the winter. The snow changed from white (albedo of 
70%), with most of the incoming light being reflected at the start 
of the melt season, to blackish in appearance (albedo 24%) by July. 
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Graphs 3A to 3G show this change in albedo over the season. 
Graph 3A, the May 26, 1973 albedo profile, shows the 
contrast between the high albedo of snow covered areas early 
in the melt season and the low albedos of adjacent tundra areas 
with little or no snow cover. From CDO-35, the snow was from 0 
to 25 cm deep, with a 10 cm thick ice layer under the snow and patches 
of exposed vegetation. CD15 was covered with yellowish ice, CD20 
contained clear ice in a depression, CD30 was a mossy area free 
of snow, and CD35 was covered with greyish ice. The average 
albedo for this area of shallow snow and ice was 21%, meaning that 
at least 75% of the incoming short wave radiation was absorbed. 
The resultant rapid melt of snow is manifested by the fact that 
all these points were free of snow by May 29, yet the May 19 snow 
depths had ranged from 8 to 33 cm (see Table 4) and the air 
temperatures during this period were relatively low (see Table 6). 
An early input of solar radiation into upper tundra areas is indicated. 
In contrast, the average albedo of the snowbed May 26 was 
59% (CD40-125), with values as high as 70%. Lower albedos (53 to 512) 
were recorded in the creek gulley (CD95-100) where accumulated 
debris was scattered over the snow surface. Around shrub willows 


protruding through the snow in this area, the albedo was 50Z. 


By June 10, the snow extended from CD50-115. CD110 contained very 


muddy snow due to runoff by the adjacent creek. CD105 was a locally 
exposed patch of mossy wet ground (albedo 37%). Melted tundra areas 
of the S-facing slope had a fairly uniform albedo of 24 to 272, 
averaging 23%.The snow, in contrast, still had a high average albedo 
of 46%, while the average albedo of the NW-facing slope and south 
hummock areas was 24%. Exposed areas near the snowbed were absorbing 
about 25% more of the incoming short wave radiation than the snow 


surface, which was still reflecting about half of the incoming light. 


This distribution of albedos changed little by June 13, when 
the snow, from CD50-95 had an average albedo of 44% and adjacent 


areas kept within the same albedo value ranges that they had June 10. 
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By June 29, though, the snow, from CD60-90, had decreased 
in albedo to 27Z, an albedo almost equal to the average albedos for 
the NW facing slope and adjacent tundra areas (26%) and the now 
exposed lower meadow and willows (21% average). The snow now was 
reflecting only about one-quarter of the incoming radiation. 

The July 4 profile, the last profile done while snow 
remained (see Graph 3E), showed that the snow had kept its 272 
average albedo. Debris accumulating on the snow surface coloured 
the snow blackish in appearance. Readings taken July 8 on the last 
remaining patch of snow gave values between 19.5 and 25% albedo. 

The decrease in albedo during the summer, with the correspon- 
ding increase in absorption of incoming solar radiation, is one of 
the several factors which accderates the daily average melt rate 


over the summer (see Graph 4A). 


Ambient Air Temperature 


As Graph 4B shows, there seems to be a definite correlation 
between an increase in air temperature and an increase in daily 
average snow melt. Low temperature can override high incoming 
solar radiation, such as the period May 22-28, when the solar input 
was high (556 eaou days but the ambient average air temperature 
was low (3.6°C) and the melt rate was low (1.4 cm/day). Warmer air 
temperatures, not always corresponding to locally high incoming solar 
radiation (see Table 6), accelerate melt rates by the transfer of 
sensible heat to the snow surface. As explained before, though, 
the combined influence of wind, albedo and snow density also 


have an effect here. 
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Relative Humidity 


As can be expected, Graph 4C shows a trend towards increased 
daily average snow melt as daily average relative humidity decreases. 
A decrease in water vapour over the snowbed meant an increase in 
the gradient of water vapour between the air next to the surface 
of the snow (this air tends to be saturated with water vapour) and 
the air a few centimeters above. This gradient produces an upward 
flux of water vapour. The rate of sublimation over a snowpack 
and the rate of evaporation of melt waters are increased as ambient 
relative humidity decreases. Another obvious factor producing 
the relationship between relative humidity and melt rate is the 
tendency of days of lower solar input to be days of high relative 


humidity and low snow melt and vice versa (see Tables 1 and 6). 


Creek and Spring Runoff 


Pictures Dl and D2 show the sudden start in 
the spring flow of the creek at the base of Snowbed *1 the night 
of May 25, 1973. Water first flowed over the surface of the snow 
in the creek gulley and eroded away enough snow to form a trench 
one to three meters deep through the snowpack. This gulley 
continued to increase in width until such time as the surrounding 
snow melted away and/or the creek had worked its way under the snow 
to its natural bed. Erosion of the snow by the creek cut the snowbed 
into two unequal pieces on the north and south sides of the creek, 
and opened up caverns up to 60 cm high under the snow. Exposure 
of the sides of the gulley and holes under the snow to the ambient 
air and to sunlight most likely increased the snow melt rates in 
these areas. 

Water from snowmelt on surrotnding hillsides early in the 
season flowed in drainage pathways or, during periods of particularly 
vigorous melt, over the tundra vegetation and under the snowbed, 


wearing drainage channels between the snow and the ground. This, as for 
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the creek, created passageways (up to 10 cm high) for warmer air 

to move under the snowpack. The Significance of the air movement 
probably wasn't too great, as usually the bottom of the snow formed 
a dense ice layer, preventing upward movement of air. In many spots, 
also, the snow was frozen right to the ground, with an ice layer 5 


to 25 cm thick (density of 0.67 to 0.73 onene 


Incoming Short Wave Radiation 


As Graph 4B and Table 1 show, periods of high daily average 
incoming short wave radiation tended to correspond to periods of 
high daily average snow melt. Radiation values in excess of 700 cal/ 
one aay were not reached during the melt period due to the local 
weather of the region which had many days of fog or high percentage 


cloud cover (see Table 6). 


Local variations in incoming short wave radiation, can be 
seen in Graphs 5A to 5F. Values were obtained with the Mk-6 Sol-A- 


Meter, sensitive to the range of 0.35 to 1.15 microns. 


Graph 5A, the May 26 profile, was drawn from data obtained on 
a clear day from 1 to 2 pm, and shows higher incoming short wave 
values for the slope of the snowbed (slope angle 10 to DO) than for 
the snowpatch area, CD0-35, with a 3 to 5° slope angle. This increased 
solar input for CD40-105 did not increase the melt rate of this area 


over CD0-35 because of the difference in albedo (see Table 1 and Graph 3A). 


The June 10 profile, Graph 5B, shows a significant difference in 
incoming solar radiation around solar noon for the NW-facing slope 
(CD120-140) compared to the remainder of the snowbed. This lower 
input of energy contributed to the observed slowness of snow border 


retreat along the N facing slope. 


The July 4 profile, Graph 5D, again shows this decrease in 
radiation reaching the NW facing slope. Both the July 4 and August 1 


(see Graph 5E) profiles indicate that the radiation hitting the steep 
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SE~facing slope of the snowbed is higher that that striking the 
flatter tundra areas adjacent to the snowbed (higher around solar 
noon, the time of highest intensities of incoming radiation). This 
implies that as the snow decreases in volume and increases in 
slope, the amount of radiation reaching it should be higher than 
while it was at a lesser angle. This might be another factor 


affecting the increase in daily average melt rate shown on Graph 4A. 


Daily Average Melt Rate 


In conclusion, though the melt rates varied locally 
depending upon albedo, snow density, solar input, wind and runoff 
patterns, the actual melt rates were remarkably consistent across 
the snowbed (see Table 1). The greatest variations occurred at 
CD100-115, on the flat creek area, where the snow, from the beginning 
of the melt season and throughout the summer, had a lower albedo 
than the rest of the snowbed. 

Of all the factors considered, it appears that incoming short 
wave radiation, ambient air temperature, ambient relative humidity, 
and decreasing albedo had the most obvious effects on the snow melt 
rate. Graph 4A indicates that the actual daily average melt rate 
had a tendency to increase through the summer, probably due to changes 


in albedo, air temperature, and incoming short wave radiation. 


Ground Temperature Regimes and Active Layer Development 
Active Layer Development 


The effect of the snowbed on active layer development is readily 
apparent on Graphs 6A and 6B and Table 3. Beyond CDO and CD115, in 
tundra and lobe areas, active layer development ranged from 20 to 60 cm in 
hummock areas and hollows between solifluction features and to 75 cm on 
the bare mud and rock surfaces of some of the lobes or frost boils. 
Snow cover did exist in the hollows but did not usually exist over 


winter for the lobe tops, the upper three-quarters of the NW-facing slope, 
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GRAPH 6A 
ACTIVE LAYER DEPTH vs SNOW DEPTH “A” 
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and the tops of most hummocks and tussocks in the region south 
of the snowbed. True active layer depths were not obtained for the 
NW-facing slope due to the presence of bedrock at or near the 
surface. In the area of CDO to CD25, where winter snow cover ranged 
from 24 to 33 cm, active layer development was between 45 to 56 cm. 
CD30 and CD35, which had 8 and 30 cm of snow respectively, had 
significantly deeper active layer development (65 and 74 cm) than 
CDO-25. The border of deep snow was at CD40. There appears to be 
a perimeter effect by the snowbed on the area within 10 meters of 
the deeper snow. From CD40 through to CD105, the May depth of snow 
ranged from 84 to 509 cm, yet active layer development seems consistent 
between 110 and 130 cm. Discrepancies exist in rocky areas (CD50, 110 
and 115) where the true active layer depth was difficult to obtain, 
and in the area near the creek (CD105) where active layer development 
was comparatively very deep (168 cm). Data from soil pits dug to 
permafrost, listed in Table 3, can be used for comparison. The 
shallower depths obtained with soil pits for the upper willows is most 
likely a localized effect, since other active layer profiles gave depths 
of 98 to 100 cm for upper willow sites. 

The consistency in active layer development in the snowbed 
(see Graphs 6A and 6B) is difficult to explain as the drainage, water 
regime, ground temperature regimes, incoming short wave radiation, 
albedo, vegetation cover and snow depth vary considerably across the 
snowbed. One possible explanation is that the presence of snow depths 
greater than one meter gives enough of a uniformity in insulation from 
winter cold that the annual temperature regime of the entire snowbed 
area is relatively consistent and that the other factors modifying active 
layer development are either cancelling each other out or of lesser 


importance. 


Ground Temperature Regimes 


Ground temperature regimes, though an expression of regional 
climatic conditions, are affected on a microscale by; local topography 
(slope exposure to incoming radiation and drainage); time, depth and 


duration of snow cover; water regime; plant cover; and type of soil. 
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The snowbed studied was divided into nine nlajor areas, as outlined 

earlier and the above mentioned factors are outlined for each 

area in this section and in Table 4. Temperature regimes for 

the 1973 season with corresponding data for September and November 

1972 are shown in Table 4B. Since the soil temperature meters 

were non-recording, gaps in Graphs 7A to 7I and 8A to 8F, including 
June 18-20, July 12-21 and July 25-27 resulted when the author was 


at other field locations. 


Upper Lobe Area 


Located above Snowbed *1 on the southerly facing slope of 
Ridge 1, this area had 0 to 15 cm snow cover in the winter (except 
for hollows between the solifluction features). On May 14, 1973, 
snow cover was lacking over those probes used for upper lobe area 
temperature determinations. The top 5 cm at this time showed 
quick freezing and thawing, even this early in the season, but at 
50 cm, the ground temperature was -3°C. November data suggests 
that this area is even colder during most of the winter. 

Throughout the summer, surface temperature fluctuations 
were limited to the upper 10 cm of the soil, but these variations 
were rapid and over a wide range of temperatures. These fluctuations 
were most likely due to the reduced plant cover of the area - bare 
frost boils interspaced with low lichen cover and a very rocky ,well- 
drained soil with reduced organic content. Heat penetration below 


30 cm was slow, with soil at 50 cm depth remaining within a6 


Upper Birch Area 


The upper birch band, located at the top edge of the snowbed 
and having a thin snow cover, was also an area of rapid surface 
temperature fluctuations but lower temperatures at depth. Once free 
of a snow cover, this area had maximum temperatures recorded at 
the 5 cm depth as high as 15 to 20°C. Some of the fluctuations 
later in the summer penetrated as far as 10 cm. This moderation 


in surface temperatures in most likely due to the cover 
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oGebireh shrub sand leaf jitter. The shrub 

cover, along with the shallow snow cover during the winter, are 

two of the main factors which produce a cold temperature regime 

in this area. The soil below 30 cm remained below o°c until 

July 11, and the soil at 50 cm €2.7° on May 18,)did not ex- 

ceed. 1 © until after July 24. It fell below 1°¢ again on 

August 16. Interestingly enough, the average albedo of the 

surface under the birch shrub was 18.6%, a low albedo not reflected 


by heat penetration deeper into the soil. 


Upper Willows Area 


Not exposed until the first week in June, this area had 
a 100% cover of willows up to 13 meters high, with 5 to 10 cm of leaf 
litter. Usually the first area of the snowbed to accumulate snow 
(up to 1 meter in summer snow storms), by May, 1973, the snow cover 
averaged 2.25 meters thick. The insulating qualities of this snow 
are shown in the temperatures of the first 10 cm of the soil on May 1 
(5 cm: ae, LOvtem* Brac eKe ys Even: at) a-depth of 50 cm, the 
temperature was 1°C warmer than the birch or upper lobe areas at 
the same depth. Once free of snow, the upper willows also 
showed temperature fluctuations in the top 10 cm, not penetrating 
to 30 cm, but sometimes exceeding a range of 20°C. At 50 cms 
the sotivdid not: reach 1°c until July 9 and never went above 5C. 
This area, though somewhat warmer than the birch or upper lobe 


localities, also had a relatively cold temperature regime. 


Upper Meadow Area 


Covered with over a meter of snow in November, 1972, and 4.5 
meters in May, 1973, this area showed a considerably warmer temper- 
ature regime than the upper willows. Though the data from September, 
1972 showed the upper willows with warmer temperatures than the 
meadow, by May, 1973, this situation had reversed. The major factor 
here is most likely the greater snow cover, even as early as November, 


over the upper meadow, as both areas are well drained. 
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Temperatures for the upper meadow area at 50 cm never rer 
below -1°C and were above ec by the end of June. Surface heat 
penetrations seldom went below 10 cm but the general temperatures 
at 30 to 50 cm for most of July and August fell in the 5 to TOC 
range. Summer factors - a low ground cover of 5 to 10 cm high 
willows with an average albedo of 26.3% also contributed to the 
relatively warm annual temperature regime of the upper meadow, 
where temperatures, even to a depth of 50 cm, remained close to 


the freezing point throughout the winter. 


Upper Late Melt Area 


Readings in this area did not start until July 1, due to the 
burial of wires and supporting pole under 5 meters of snow at the 
beginning of the season. 

Temperatures on July 1, to a depth of 50 cm, were at 
freezing point. Whether these warm temperatures are due to an early 
snow cover of 160 cm by November and a deep winter snow cover of at 
least 5 meters, or due to melt waters and warm air passing under the 
snowbed during the melt period, raising the soil temperatures to a 
freezing point by July 1, is difficult to say, though the former 
seems the better explanation. Once melted out, this well drained 
area of chiefly bare mud and rock showed rapid, deep temperature 
fluctuations even to a depth of 50 cm. Once exposed from snow, the 
soil at 50 cm rapidly reached a temperature range between 5 and 152: 


This area showed the warmest temperature regime recorded for the snowbed. 


Lower Late Melt Area 


Lower late melt localities, on the lower, flatter slopes of 
the snowbed and usually saturated with water, had less snow cover than 
the upper late melt area throughout the entire winter. Not as warm as 
the upper late melt area, heat penetrations were not as deep or as high 
in temperature as the upper late melt, probably due to the shallower 
snow cover and the moderating effect of the water present in the lower 
late melt section. The 1°C mark was not passed at the 50 cm depth 
until July 10, but temperatures at this depth did remain in the 5 to 10 °¢ 


range for at least part of the summer. 
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Lower Meadow Area 


The lower wet meadow, with 1.2 meters of snow in May, 1973, 
showed a considerably colder temperature regime than the upper 
meadow area, exposed two weeks later but having 4.5 meters of snow in 
May. Situated at the lower side of the snowbed, the lower meadow 
tended to be saturated all summer with runoff water from the snow. 
Upper surface temperature fluctuations, once the snow cover had 
melted off, seldom reached below 10 cm. Surface temperatures 
tended to be cooler than the late melt areas and upper meadow area, 
probably due to the presence of standing water over much of the 
season. At depth in the soil, temperatures were much colder than 
the late melt areas or upper meadow, with the temperature at 30 cm 
being below 1°C until June 28, and below 1°c at 50 cm until ASE ba ae ae 
The thinner snow cover over the winter, compared to the other three 
areas, allowed penetration of colder temperatures into the ground 
in this area. 

In the later part of the summer, the lower meadow maintained 
a relatively warm temperature regime, with the daily average temper- 
aturesat 50: cm. from auouse loto 1/7) at ec te compared to 7.8° for 
the upper late melt, 4.9° for the upper meadow, and 5.2° for the 


lower late melt area. 


Lower Salix Cham. Area 


Though comparable in vegetation to the upper meadow area, 
this site was sometimes saturated with runoff water, being on the 
lower side of the snowbed. For most of the summer, because it was 
on slightly raised ground, the site was better drained that adjacent 
lower meadow areas. Covered with 1 meter of snow by November, 1972 
and 2 meters by May, 1973, the lower Salix Cham. area had temperatures 
around = lea on June 1 at all depths, compared to temperatures around 
freezing for the upper meadow. 1°C was not totally reached at 50 cm 
until July 4, compared to June 28 for the upper meadow. Daily 
temperature fluctuations seldom went below 5 em, and only occasionally 


reached 10 cm in depth. The surface temperature range was lower over 
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the summer for this area than the upper meadow, and the August 
1-17, 50 cm average of 30h Cis again lower than that of the 
upper meadow (4297 Oe The moderating effect on the lower Salix 
Cham.area is probably saturation of the ground with runoff water 
throughout much of the summer, compared to the well drained 


conditions of the upper meadow. 


Lower Willows Area 


The lower willow area, in spite of a snow cover of only 1.6 
meters in May, had a warmer temperature regime than the upper willow 
site, which had a May snow cover of 2.25 meters. Temperatures in 
the lower willows were not below -1°C at 50 cm throughout the winter, 
and exceeded 1°C June 28. Good cover, provided by the shrub willows, 
kept surface fluctuations in temperature usually below 20°C. Most 
fluctuations scarcely reached 5 cm in depth, and rarely reached 10 cm. 
The temperature regime at 50 cm was warmer than that of the upper 
willows - the August 1 to 17 average was 4.4°C compared to 3.2°C for 
the upper willows. Differences may be due to the presence of the 
creek within 10 meters of the lower willows and the passage of runoff 
water through or near this area throughout the melt season, or perhaps 
the influence of colder tundra areas above the snowbed upon the soil 


temperature regime of the upper willows. 


Soil Temperatures at 50 cm 


Soil temperatures at 50 cm depth do not generally show the 
rapid temperature fluctuations found in the surface 10 cm of ground, 
but tend to reflect the overall heat balance of an area. The 
magnitude of temperature penetration to greater depths, such as to 50 cm, 
has a direct relationship to the development of the active layer in 
a site. Graphs 7A to 7F show the temperature regimes, at 50 cm, of 
the major areas considered. It is apparent from the graphs that the 


areas above the snowbed or at the periphery are much colder at 50 cm 
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GRAPH 7A 
SOIL TEMPERATURES — 50 CM 
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GRAPH 7B 
SOIL TEMPERATURES — 50 CM 
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GRAPH 7C 
SOIL TEMPERATURES — 50 CM 
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GRAPH 7D 
SOIL TEMPERATURES — 50 CM 
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GRAPH 7E 
SOIL TEMPERATURES — 50 CM 
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GRAPH 7F 
SOIL TEMPERATURES — 50 CM 
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than the areas within the snowbed (i.e. the birch, upper lobe and 
hummock areas were colder than any snowbed sités). Within the 
snowbed, the upper willows were colder than the lower willows, yet 
the upper meadow was warmer than the lower Salix Cham. site. The 
upper meadow appeared to cool off faster than the wet lower meadow 
(Graph 7C), while the upper and lower late melt areas, once exposed, 
had the warmest temperature regimes at 50 cm of ail the locations. 
Factors determining these results have been discussed in the previous 
sections. One can conclude, though, that vegetation cover (which 
provides insulation from incoming solar radiation), water regime, and 
nearness to tundra areas, have the greatest effects on soil temperatures, 


at depth, within any one area. 


Ice Content and Stability Of Underlying Soils 


Ice content of soils found in snowbed areas below the active 
layer depth was not determined, as a drill was not available for 
probing into the permafrost. Due to the presence of either bedrock 
near the surface and/or the rockiness of most soils within the region, 
soil augers were ineffective in core drilling so that the structure 
and water contents of soils were examined by digging pits to the level 
of permafrost. 

The soils within and around Snowbed *1 showed many signs of 
instability or of factors leading to unstable situations. Many 
snowbeds similar to Snowbed *1, form along creek gulleys or gentler 
slopes. In the rugged topography of the Richardson Mountains, these 
gulleys and gentle slopes represent potential paths for movement of 
ground vehicles or potential areas for gas pipeline construction. 

Such areas of deep snow accumulation represent a hazard to any such 
activities. To begin with, the wet soils created by snowbed runoff, 
the deep active layer depths, and the soil temperatures near 0°C 
throughout the winter, combine to create a potential for disturbance 
by vehicular activity. 

The steep banks and runoff waters associated with most snowbeds 


give them and adjacent areas a natural tendency for unstable soils. 
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Picture El shows a mud slide which occurred on June 1, 1973, on 

the steep NW-facing bank next to the snowbed. Even more dramatic 

in 1973 was the collapse, to a depth of two meters, of ground four 
meters from the natural bed of the creek at the bottom of Snowbed *1. 
During spring runoff the creek cut an underground channel through the 
lower willows, emerging 20 meters downstream through a large hole 
about three meters from the natural creek bed, and cutting a new 
channel across to its normal gulley. This underground channel, which 
was not used in 1971 or 1972, suggests the presence of a large ice 
lens which was rapidly eroded by the pressures of spring runoff. 

The presence of springs which suddenly emerge halfway down 
the slope or at the bottom of snowbeds was common. These under- 
ground water drainage passageways suggest also a movement of water 
along weak areas in the permafrost. 

In summary, many of the snowbeds of the Canoe Lake area 
represent areas of wet, unstable soils throughout the summer, and 
areas of slow freezeup and high water content soils in the winter. 
Danger of damage to peripheral areas near deep snow accumulations is 
suggested by the deeper active layer depths here than at sites 
further away. Both the willow and birch bands, upper meadow 
and especially lower meadow areas could be damaged by removal of 
snow. This would expose the plants to colder winter temperatures 
than normally experienced under the insulation of snow and to wind 
abrasion and damage. Compaction of snow could lead to breakage of 
shrubs, an increased melt period, and increased conductivity of cold 
into the snow (Sellers, 1965). 

Snowbeds of other Arctic regions also represent potentially 
unstable areas. Picture E2 shows snowbeds found in the Shingle 
Point, Yukon area along the Arctic coast. Fresh mud slides and old 
mud slide scars are common about the snowbeds, as well as large cracks 


and sliding pieces of vegetation (see Picture E2). 
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Estimation of the Energy Flux Over Snowbed 


Valuesof wet and dry bulb temperatures were taken at the 
center of the snowpack during the melt season in June to determine 
estimates of the energy fluxes over the melting snow. On June 13, 
from 11:00 to 17:00 hours, values of wind, albedo, incoming short 
wave radiation, and wet and dry bulb temperatures were taken hourly. 
The average density of the surface snow was also determined. Table 5 
outlines the values and estimates obtained. 

The net energy balance over a snow pack can be determined by 


= 
the equation R R,% = rk y + RY = Rt e oT. 


L 
where Beg net radiation 
RoY = incoming short wave radiation 
Ee short wave albedo 
RY = incoming long wave radiation 


reflected long wave radiation 


[an 


banal: = Stefan-Boltzman equation - the emmittance of heat from 
a body of surface temperature T. 
T. = 0°C or 273°K for snow 
ea 7 emissivity = 0.90 for old snow 
June 13 = clear day, without cloud cover. 


Incoming short wave radiation between 0.35 and 1.15 microns 
was obtained with the Sol-A-Meter MK6. Incoming long wave radiation 
was estimated as 1% of the energy value of RY (Sellers, 1965). Most 
likely this is an underestimate and values presented here represent 
a minimal range. 

Using the energy balance equation R -G=LE+H 

where G = flow of heat downwards into the snow 
LE = latent heat flux 


H = convective heat flux 


estimates of the energy balance over the snowpack were obtained 


(Rouse & Kershaw, 1971). G was determined with thermistors placed 
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at 5 and 10 cm depth in the snow, and with the equation 


Be AT 
G r iG 


A = thermal conductivity = ep (Geiger, 1966) 
p = snow density = 0.499688 pee: 

c = constant= 0.0066 cal/cm/sec/degree 
AT = change in temperature 


Az = 5 em change in-depth 


Ro - G, then, was a measurable value. Using the Bowen Ratio 


method where B =H = -pc Kat Cane 
LE 
A L Aq 
Pee Srey 
eS = heat capacity of the air at constant pressure 
L = latent heat of vaporization 
AT = gradient of temperatures 
Aq = gradient of specific humidity 
Gd; 9 U022 6) 
Ep 
a 
e = water vapour pressure (millibars), determined 
from wet and dry bulb temperatures by the method 
outlined in Smithsonian Tables, 1971. 
P. = atmospheric pressure = 1010 millibars 


Both temperature and water vapour gradients were measured by sets of 
wet and dry bulb thermometers placed at 150 and 50 cm above the snow. 
As can be seen by the results, difficulties with air inversions 
over the snowbed and water vapour condensation on the surface of the 
snow (a downward flux of water vapour), were encountered on this and 
other days when attempts were made to measure energy flux values. 
According to the literature (Geiger, 1966; Munn, 1966) work done over 
glaciers, especially the Blue Glacier in the United States and the 
Fedchenko Glacier studied by Kazansky, has shown that condensation over 
snow surface is a common phenomenon and latent heat released in this 


process is used to melt snow. The results obtained at Snowbed ¥*1 suggest 
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a large input of energy from water vapour condensation and cooling 

of the air, into the snowbed system and into melting the snow, not 
evaporating water or heating the air above the snowbed. Further data 
analysis and study is required, though, before any definite conclusions 


can be reached. 


CONCLUSIONS 


During examination of the physical characteristics of snowbeds in 
the northern Richardson Mountains, N.W.T., formation of snowbeds was 
found to be mainly an expression of prevailing wind (NW) and local topo- 
graphy, with the sudden decrease in wind speed when passing over a steep 
drop in topography, leading to drift formation, especially on slopes 
facing south to east. Melt rate of the snow appeared to be most influenced 
by ambient air temperature, relative humidity, albedo and incoming solar 
radiation, though other factors such as wind, rain, density and creek 
runoff had some influence. The daily average melt rate was found to 
increase over the summer. 

Active layer development was influenced by the presence of snow 
cover greater than one meter, but active layer depths for areas with snow 
from one to five meters deep were remarkably consistent. Ground temperature 
regimes, on the other hand, more expressed the depth of snow cover and the 
water regime of each area. Important factors affecting ground temperature 
were found to be plant cover, drainage, depth of snow cover, time of snow 
cover and duration of snow cover. 

Slopes within and near the snowbeds showed signs of natural in- 
stability manifested by slides, underground springs, and slumpage of 
ground in wetter sites. Snowbed soils, with temperatures near freezing 
all winter and high moisture regimes during the summer, could be more 


easily damaged by human activity than adjacent tundra areas. 
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TABLE 2 WIND AND RAIN RECORDS, 1973 


DATE RAIN (cm) WIND km/hr 


werereee2reqgeeacec= eeeerenecereeceeeeceecesece 


OIRECTION 


May 17 NW-NE 
18 E 
19 2.32 E 
20 3.38 SE 
21 15.20 NE-NW 
22 1 0.140 | O-278 pe 5.12 NW-NE 
23 | O 4.97 NE 
e4 | 0.564] 0.561] x 8.33 SW=MW 
2S || Manos) || CQ i@Zi) se 10.88 51 | Nw 
26 | O fe) 9.46 6.71 | WNW-SE 
2A \\ Oss) O5855)) mE 15.95 8.08 | NNW-NW 
28 |- 02268 1 “0,291) x 11.09 6.15 | SE-NW 
29 | 0.756] 0.726! x Sal, 2.88 | NY 
SOM MOR SoulmmOsOL2 |G Weoe 3.97 |SE-E 
41 | -07932"|' 0,829!) Xx x 7.52 3.43 | Nw 
June 1 | 0.002! 0.006] x . 4.81 | NW 
OOO Os OO! | mac 20.65 9.98 | Nw 
D |e0¢b50 |) 6 Os205| =x. 11.40 5.85 | Nh 
4 | 06085 | 0.070) XxX 5.43 4.18 | NW-SE 
3 || © O 8.05 6.75 |SE 
6 | 2020014) 0,001) 7X 5.63 4,93 | NW=-NE 
YN © | pia Deg, 6.60 |NE-NwW 
Sure fe) 25.92 11.93 | Nw 
9/0 (@) Pea Gal 5.02 |NW-NE 
20. 16 ) 6.73 6.52 |ESE 
nal |} © ) 16.13 8.65 |SE-NE 
Wa |1 © fe) 15.49 7.76 |NW 
13) oO fe) 10.24 Zaeitinl she 
Ty || o fe) 19.56 5.78 |NW-SE 
Ls} CO) Clon || Cen) 9% 32.39 7.63 |NW 
MS |L @Oo@os | CWpKOrerxt| se 8.31 3.21 |SE-NE 
a7 
TOM OL.205') (Oseor) x 12.24 NW 
19 
20 
AN IP @a25s) || OnaQsrh|) s% 17.89 SE-NE 
22, 1 2.776 | 2.505). x 16.73 NW 
223 || Ose || WABES se x 28.59 Nw 
24% | 0.109} 0.103) x 6.99, NW-SE 
25. | 0.007 | 0,004) xX 8.51 NW-SE 
26° 0.154) | GO,227) xX 8.71 2.71 |NW-SE 
2m @) 9.12 6.50 3.20 |NE 
30 | O fe) 10.46 6.88 3.64 |KE 
sulLy oy 10.182.) “Os2901" xX 12.92 7.98 3.20 |SE-NE 
fe) fe) &.54 6.34 2.48 INW 
Sl OsOOls Oso] (x 7.23 6.60 3.63 |SE 
hl (0) fe) 16.62 11.99 4.58 |NE~SW 
Sh OsS16 0, sie) x 16.26 22657 3.80 iNW 
6 1 0.258 | “0.205! xX 6.51 5.77 2.05 |SE 
DIN RL EOS 1 250051 5x 9.54 7.46 2.15 |SE-NW 
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TABLE 2 (Continued) 


5 ee -— = = igs ee 


8 | 0.063 0.055 x 5.79), 4.77 1279) SE 
9] 0.032 0.032 x F200) eee 1.73 | NW-NE 
107° ) 17.28) 12353 3.91 | SE-NW 
Liar ial esses 6.584 x 12.0) 21337 *2.91 
blown 
over 
22 | 0.639 0.647 x O77) 7,58 1.14 | SE-NW 
en fe) SLO 7a 2.00 | SE-NE 
oR? Ne fe) 15.69 12.31 3.37 | SW 
28 | 0.884 0.876 x 11.9] 8.96 2.45 | SE 
20 | 0.691 On777 x 12.79] 9.49 2.98 | SE 
30 | O ) ao Ji 2Oce7 5.84 WNW 
425-9 fe) 11.66] 9.56 3.04 | SE 
Aug. se Naan O 23-99) 37.70 4.71 | Welw 
Bit 20 fe) 9.45) 8.71.) 2.82 | WNW-SE 
3 | 0.339 O2555 x 14.02| 8.08 312 | SE-NE 
i fe) 9,24] 6.93 2.13 | NW 
Sti cOstlos 02072 x 13.25] 8.44 228 15S 
6 | 0.038 0.049 x 30.53| 20.82 6.17 | NE-NW 
7 | 0.860 0.695 x PEL9S) 19.42 5.70 | wW-NVW 
8 | 0.006 0.006 x 10517). "Gs6e2 1.88 | HW 
So. 1 0,00! 0.001 x 28.42| 24.36 7.38 |} t 
{O°} 0,071 CO, 071 x 20.47| 16.02 4.05 | WNW 
TP OoGer 0.033 x 15.46| 10.14 ope an bet 
12 | 0.987 1.026 = 7.25) 9209 3.14 | SE=-NW 
13 | 0.024 0.032 ye 20.39| 14.06 B21 | NY 
14 0.008 0.016 X 20.23| 14.17 3.96 | NV 
45° 0.001 0.016 x 8.55} 5.95 2.47 | NW 
16.0 ) 13.92| 10.27 2.68 | SE 
py OM Wi Ge O.? x SE-WNYW 
Total 
May-Aug. 


16|20.647cm |21.116cm 8021km |22047kn| 9218km 
2158hr} 2230hr 2204hr 
12.98] 9.88 4.18 


km/7br| km/hr | km/7or 
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TABLE 3 ACTIVE LAYER DEVELO: MENT 


DEPTH OF SNOW ACTIVE LAYER DEPTH (cm), 1973 
TIME OF 4 
EXPOSURE F#- m4 —-—-------- ae el 
: : : : : ; * SOIL 
ie MAY29°0UN : ers 
A ‘ A AVERAGE = Sea 
4 
Uoper O} 15] 25 33 | May 29) 6 16 ZAL Zot Bi eeu! Lobe-~ 
lobe 25h." "| 24 | May 29] 8 17} 26] 441 48/56 | a2 
And 10) ° esa | ay eo «7 17 es. | oO ee 156. | 60 
Heath 15 4 M aL hay 29] 7. 16 | 25] 326 {| 39/45 4 50 |Heath 
20s ees Se ee ee U7} eh oS Ae. ESOn of 20 
Pos " | 321 May 29] 8 ED esi ea Aes fe 
s 
(ransi~ : 
tron oO a> 29 8 | May 29] 13 22 40 58 64 | 65 490 270 
35, 30 | Nay 29] 7 17? | 33 | 50) 74 | 74 Birch 
/ 
Upper ; 
Willows , 
4o| 70 | 50] 84] June 5 36 | 57 | 68 | 92 |132 4409] 60 
451120 97 | 174 | Junel2 32 72 Ties Ni abawss |) alex (ie 
f 
s 
Upper , 
Meadow { 
50} 120 276 June 21 45 | 62 | 82] 83 4 120 
55| 120 348 pune 25 44 | 73 |118 |120 $445] 120 
60 417 July 1 20 C7 ao a eae } BD 
f 
f 
Upper i 
Late 65/10 -1160 |496 puly 8 46 82 1122 S456 
Melt 70|20 509 puly 22 95 |130 : 
Lower 
Late 75110 - 485 uly 11 5 ayey abate). 1) 
Melt 80)20 452 Duly 11 18 84 1113 4/112 
85 408 July 8 25 pi eGeliaany 
y 
s 
Lower 4 
Meadow f 
90}10 - 358 July 1 15| 30 | 97 1127 4449 | 100 
95/20 | 63 | 249 fuly 21 38| 68 |107 |110 | 100 
s 
f 
Lower s 
haaacts ae 9 \a 4 a i 
100 ood ete ele) geese (2) 1205 4105 
105/20 |101 | 196 | June 1 a2 183 has lies li6a (237 | 220 
s 
s 
Creek s 
Area 110 234 |June 15 71} 44155 | 90 | 90 : 85 
EES 378 |June 21 a4 | 54 SOM |S Onn 
s 
4 
| | 
4 
| f 
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TABLE 4A SNOW DEPTHS SNOWBED*1, 1973 


DEPTH OF SNOW DATE OF 


(cm) EXPOSURE | 
LOCATION D = |---------------------------- ee FROM | WATER REGIME 
Desh Sept.19 |Nov. 24|May 18 | RECORDED DAY SNOW 
19722 Ve lOge lr 9798 1973} 1973 


Upper Surface 
Lobe 5 


0 15 0 0 | May 2 
18 


Well drained. 
Upper slopes 
above snowbed. 


Upper Surface} 30-60 25 20 20 May | May 26. | Well drained. 
Birch 5 0 0 18 ? Upper edge of 
10 0 0 snowbed. 
30 0 
50 0 


Upper Surface] 120 ie 225 225 May June 5 | Well drained. 
Willows 5 19 | Upper edge of 
10 snowbed. 


10-20 1 3yz 450 98 June] June 21 
up to is 
150 cm 


HSe Sia) 


Upper Surface 
Meadow 5 


Well drained 
after exposure. 
Upper slope of 
snowbed. 


Upper Surface 
Late 5 
Melt 1s) 


10-20 160 HS 50-1004) Joly | Jaky Well drained 
after exposure. 


| 
! 
! 
| 
| 
Upper slope of 
| 
| 
| 
| 
| 
| 
| 
| 
| 


snowbed. 


Lower Surface} 10-20 - 345 88 June] July 5 | Mostly saturated. 


Late 5 25 1 Lower slope of 
Melt 10 snowbed. 
30 
50 
Lower Surfacef 10-20 63 118 Sy) May June 7 Mostly saturated. 
Meadow 5 29 Lower slope of 
10 1 snowbed. 
30 
50 | 
Lower Surface] 10-20 103 204 138 June] June 14 | Sometimes satu- 
Calix 2) ie : rated. Lower 
10 | slope of snowbed. 
30 
50 | 
Lower Surface} 10-20 74 160 87 June] June 13 | Sometimes satu- 
Willows 5 3) | rated. Lower 
10 | Slope of snowbed. 
30 
50 
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TARLE 5 = BOWEN RATIO VALUFS, SNOWBED*1, 1973 


~ISHORT WAVE RADIATION BS We ee 
Bi reores tegcta qprekt yas nc 
Pee eee ee ee NLICM2MIN 
260 1.1758i1.0.7172 0.03015} 0.027010. 7443 
ES. 1.2436]0.7710 |0.03273| 0.0295/0.80048 
290 1.3114]0.8262 |(0.03544! 0.0319 /0.85809 
410 1.4019/0.8832 0.03789] 0.0341/0.91723 
295 Wag S4O OL C27 0. OF 51sO. OS16 | O.S5E70- 


260 1.1758)0.7172 0.05015] 0.0271 |0. 744-54 


0. 5659 


215 a as coe ad pacers? 0.0199 
=e pao nered = 7 x LR om ene eeaaer meetin a tae , = eer 
FLUX INTO. | TEMPERATUR ADJACENT 


VEGETATE 
SNOW (G) Lee AREA (CO 50) 


TIMED Pigs ee ey Gy nokes eiee eo Sarma at re eas face 


for G ! CaLicM2/miNn 


00 PO ipl) ae 
iB gr [tess] 19 
ee 5.8 |15.6| 45 
| 13200 eis 
\ Sees os 
fae 00 1.7 |129-6) 44 
i 4% 115.6) 44 
j, 15:00 LAs on 44 
l Ai [LGee 0 
|| 16:00 2 SONG oF 
i L [ 3.0 [17.2] 37 
te OO. | ate yene2 | oy aut ie 27) 350 
mm re Sa - wd RSS RRETS 
\WATER VAPOUR a 
PRESSURE SPECIFIC HUMIDITY wy | BOWEN 
(SO) mnillibats. We ON 00 ee oe oe LE H 
TIMES. [ie hoo eee iis ieee Pik 1 RATIO 
150 cm : 50cm 10cm, 50cm Aq a cal/emYymin calicm4n 


paler |) ale) A stal 


12:00 |11.873;10.554 | 7.312 |6.500 10.8127 |~0.396 

1.1436 |70.178 |-0.729 |+0.296 
2.2788 |-0.088 |-0.839 0.074 
1.9424 |+0.272 |-0.666 F0.181 


13400-32211 110,354: 7-520 \6.577 
14:00 |13,417|10.216 | 8.570 |6.292 
15:00 |a4.248]12.094 | €.774 |6.832 
1.1939 | 0.168 |-0.629 [o.104 
1.0516 |+0.424 + 0.377 + 0.160 


16:00 J15.094 113.945 | 9.782 18.588 


i] 

1 

1 
9.086 | 7.520 15.596 |1.9245 | -0.507 
19:00 113.841 112.129 | 8.524 17.469 
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TABLE 6 


DAY 


ne) 
w 
LJ 
PSPSPS PS OS PI OS OS PS PS PS PS POS PS 


June 


=) 
= 
PSPS PS PS PS PS PS PS PS PS PS PS PS PS OS PI PS PS PS PS 


be DSN 


ws 
a 


MO SS 


PL Dd bed Dd Od 


mS PS PS 
PS PPS PSPS 


mS PS PS 


Po PSPS PS PS 


PSPSPS PS OPS 


bd PS 


WEATHER RECORDS, 1973 


SNOwB"D*1 
) RELATIVE 
TEMPERATURE (C), HUMIDITY (%) 
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WIND (kmlhr) 


(cm) | 2 


SHORT 
WAVE 


RAIN |RADIATION 


SE 
NNW 
NNW 
6.5 |NW-E 
4.7/E 409.7 
Zao Sb 669.7 
2.G1SE S68 
13.8 |NE-N\ 628.8 
5.0 |NW-NO.14-46€0.8 
4.04 INE OKO) ||) Sexo 3 
6.5 |SW=-NWO. 58] 507.5 
Oe ME 0.10} 259.7 
&.5 NW O Wes 
12.7 INNW |0.36) 466.3 
9.0 INW 0.27] 625.3 
4.8 |NW 0.76) 303.6 
5. Sk OnlGeis252.5 
5.8 |NW 0.94] 125.6 
5.9 |NW 0.02} 297.4 
16.9 |NW tr 456.5 
| 9.2 |NW Qaals) Sy =k 
5.1 {Nw 0.09 3/5.2 
6.8 ISE O 535.9 
4.9 |NW 19% ADO o/ 
9.2 INE (0) 489.1 
20.6 IW O Obes 
6.6 INW O 409.6 
6.5/#SE |0 706.8 
ayer tay fia; 146) 681.2 
6 SW O Poa 
eH \tSia: (0) 706.0 
BLS 7) Net (0) 691.3 
26.3 |W tr 622.5 
Viste |SE or 613.0 
NW O 
9.3 0.2 |601.3 
Bais IS OneomGe ay, 
1.9 |NW Ne Qe eal 
9.6 |W OAat Boece 
5.4 |NW Sore eee 
6.0 |NW O.0OL 370.4 
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TABLE 7 SOME PLANT SPECIES ASSOCIATED WITH SNOWBEDS 


(Key: Hulten, 1968) 


Achillea borealis Bong. 

Aconitum delphinifolium DC. subsp. delphinifolium 

Anemone narcissiflora L. subsp. interior Hult. 

Anemone Richardsonii Hook. 

Antennaria monocephala DC. subsp. monocephala 

Arctagrostis Latifolia (RK.8r.) Griseb. war... datifolia 

Arctostaphylos alpina (L.) Sprong. (snowpatch or tundra) 

Arnica Lessingii Greene subsp. Lessingii 

Artemesia artica Less. subsp. arctica 

Artemesia Tilesii Ledeb. subsp. Tilesii 

Betula glandulosa Michx. 

Calamagrostis canadensis (Michx.) Beauv. subsp. canadensis 

Carex aquatilus Wahlenb. subsp. aquatilus 

Carex lachenalii Schkuhr. 

Carex podocarpa C.B. Clarke 

Cassiope tetragona (L.) D. Don subsp. tetragona (snowpatch) 

Castilleja Raupii Pennell. 

Dodocatheon frigidum Cham. & Schlecht. 

Empetrum nigrum L. subsp. hermaphroditum (Lange) Bécher 

(snowbed, snowpatch and tundra) 

Epilobium angustifolium L. subsp. angustifolium 

Equisetum arvense L. (snowbed, snowpatch and tundra) 

Eriophorum angustifolium Honck. subsp. subarcticum (Vassiljev) Hult. 

Eriophorum Scheuchzerii Hoppe var. tenuifolium Ohwi(snowpatch and snowbed) 

Festuca altaica Trin. 

Ledum palustre L. subsp. decumbens (Ait.) Hult. (snowbed, snowpatch and 
tundra) 

Loiseleuria procumbens (L.) Desv. (snowbed, snowpatch and tundra) 

Lupinus arcticus S. Wats. (snowpatch and tundra) 

Luzula confusa Lindeb. 

Myosotis alpestris F.W. Schmidt subsp. asiatica Vestergi 

Oxyria digyna CL.) Hill 

Parnassia Kotzebuei Cham. & Schlecht. 

Pedicularis Kanei Durand subsp. Kanei (snowpatch or tundra) 

Pedicularis sudetica Willd. subsp. interior Hult. 

Petasites frigidus (L.) Franch. 

Poa, aretica Rk. Br. subsp. arctica 

Polemonium acutiflorum Willd. 

Polygonum bistorta L. subsp. plumosum (Small) Hult. 

Polygonum viviparum L. 

Ranunculus nivalis L. 

Ranunculus pygmaeus Wahlenb. subsp. pygmaeus 

Ranunculus Turneri Greene 

Rubes chamaemorus L. (snowpatch or tundra) 

Rumex arcticus Trautv. 

Salix Chamissonis Anderss. 

Salix polaris Wahlenb. subsp. pseudopolaris (Flod.) Hult. 
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Salix pulchra Cham. 

Salix reticulata L. subsp. reticulata (snowpatch and tundra) 

Saussaurea angustifolia (Willd.) DC (snowpatch and tundra) 

Saxifraga punctata L. subsp. Nelsoniana (D. Don) Hult. 

Senecio lugens Richards. 

Spiraea Beauverdiana Schneid. 

Stellaria longipes Goldie 

Taraxacum alaskanum Rydb. 

Trisetum spicatum (L.) Richter subsp. spicatum 

Vaccinium uliginosum L. subsp. alpinum (Bigcel:) Hult. “(snowbed, 
snowpatch and tundra) 

Vaccinium vitis-idaea L. subsp. minus (Lodd.) Hult. (snowbed, 
snowpatch and tundra) 

Valeriana capitata Pall. (snowpatch and tundra) 

Viola epipsila Ledeb. subsp. repens (Turcz.) Becker 


* Not all these species are confined to snowbeds only. 
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PICTURE FAl - CANOE LAKE AREA - June 7, 1973 - SW OF SKYLINE RIDGE 


PICTURE FA2 - CANOE LAKE - July 11, 1973 - LICHEN RIDGE 
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PICTURE Fl - SHINGLE POINT July 17, 1973 - SNOWBEDS ALONG COAST 


PICTURE F2 -— SHINGLE POINT - July 19, 1973 - SNOWBEDS ALONG COAST 
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PICTURE F3 - SHINGLE POINT - July 18, 1973 - SNOWBED NEAR WALKING RIVER 


PICTURE F4 - SHINGLE POINT - July 19, 1973 - SNOWBED NEAR BLOW RIVER 
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PICTURE Al - SNOWBED *1 
July 28, 1973 
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PICTURE A2 - SNOWBED *1, August 13, 1973. CD PROFILE LINE. 
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PICTURE A3 - SNOWBED *1, July 28, 1973 
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PICTURE A4 - SNOWBED *1, July 28, 1973 
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PICTURE Cl - SNOWBED *1, May 24, 1973 
FEATURES OF THE DRIFT FORMING SNOWBED *1 
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PICTURE Dl - SNOWBED *1, May 30, 1973 
CREEK ERODING BOTTOM SLOPE OF SNOWBED 
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PICTURE D2 - SNOWBED *1, June 21, 1973 
GULLEY CUT THROUGH SNOWBED BY CREEK 
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PICTURE El - SNOWBED *1, June 6, 1973 
Minor mud slide (x), NW facing bank of Snowbed *1 
Exit of underground stream (y) 


PICTURE E2 - SHINGLE POINT, YUKON, July 19, 1973 
Unstable areas associated with late snowbeds 
Cracks in ground and sliding blocks of vegetation (x) 
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ABSTRACT 


Two small drainage basins in the Mackenzie Delta 
region, N.W.T. are under study to obtain hydrologic data needed 
for assessment of pipeline and highway designs, and to determine 
water supply conditions. Each basin is instrumented during the 
summer field season to measure precipitation, runoff and latent 
evaporation. Results for 1973 are presented, and comparisons with 
former year's findings are made. Precipitation is observed to be 
light, the annual normal for Inuvik being 285 mm. Maximum runoff 
typically occurs during the spring snowmelt period. A peak discharge 
of 0: 22 at eae was recorded at one of the basins in 1973. Direct 
runoff following summer storms is minimal, less than 5% of precipitation. 
Potential evapotranspiration for each basin for the period June to 
September is estimated at 225 to 275 mm. Water balance calculations 
for one of the basins have revealed a storage depletion of up to 115 mm 


of water, for the period October 1, 1972 to September 30, 1973. 
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PROGRESS OF HYDROLOGIC STUDIES AT BOOT CREEK AND PETER LAKE WATERSHEDS 
NeW... DURING, 1973 


J.C. Anderson and D.K. MacKay 


INTRODUCTION 


Hydrologic studies have been in progress at Boot Creek basin since 
1967, and at Peter Lake basin since 1970. Figure 1 shows the location of 
these basins in the Mackenzie Delta region. Progress during the 1973 field 
season is discussed in this report. Preliminary results obtained from 
former years' data, as well as a general description of the basins, are to 
be found in an earlier report (Anderson and MacKay, 1973). 

In 1973, the instrument networks were established basically as 
they had been in 1972. In addition, a second rain gauge was installed in a 
pit at the north end of Peter Lake, in order to perform an experiment on 
rain gauge catch efficiency. Also, the atmometer network was expanded 
Sslaenuly, and tie.250 one atmometers were replaced with 500 cm” imstruments., 
to overcome the difficulty of reservoirs running dry between instrument 


readings. Analysis of data collected in 1973 at the two basins has improved 


our knowlege of water budget components in the Delta region. 


PRECIPITATION 
Annual precipitation is light in the Delta region, and decreases 
northward toward the coast of the Beaufort Sea. Rain is the dominant form 


of precipitation from June through September. Snow predominates for the 


510)! Bice 


remainder of the year. Annual precipitation normals for Inuvik and 
Tuktoyaktuk are: listed in Table 1. 

During the summer field season, precipitation is measured at 
each basin in two locations, using M.S.C. tipping bucket rain gauges. 
Precipitation data collected by the Atmospheric Environment Service (AES) 
at Inuvik airport, 6 km south of Boot Creek basin, are used to supplement 
the records of the latter basin. Monthly precipitation totals obtained 
at Boot Creek have been in good agreement with those recorded at the 
airport in recent years (Anderson and MacKay, 1973, p. 48). 

No AES weather observing stations are located close enough to 
Peter Lake to be of any use in extrapolating records for that basin. 
However, the summer data that are available for Peter Lake suggest that 
precipitation during that season is normally slightly less there than 
at \Boot Greek (table 2). 

It is noteworthy that monthly precipitation can vary 
considerably about the mean (Table 2). For example, in 1971, only 
14% of the mean May and 35% of the mean June precipitation fell at 
Inuvik. The following October, 230% of the normal monthly precipitation 
occurred. In 1973 at Inuvik, 4 months were above normal in precipitation, 


8 were below normal, and the year itself was below normal by 20%. 


RUNOFF 

Peak runoff for the year at Boot Creek and Peter Lake watersheds 
occurs in the spring, with the melting of the snow pack. The pack 
normally reaches its maximum water equivalent in early May, following a 


winter snow accumulation period which can last as long as seven months. 
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Snowmelt runoff occurs in late May and early June, over a period of two to 
four weeks. Permafrost and seasonal frost beneath the snowpack act as a 
barrier to infiltration of the melt water, and a large percentage of the 
snowmelt contributes directly to runoff. 

The 1973 snowmelt flood at Boot Creek peaked on May 27, at a 
measured flow of 6.9 ne C245 ELS) AOL 1c? ear (20 Spey (Figure 2). 
As discharge was not being recorded continuously at the time, the maximum 
flow could have been still greater. The peak flood came in response to 
inputs of rain and heat. On May 25, 6 mm of rain fell at Inuvik. The 26th 
was the warmest day of the month, the maximum temperature being DDG (yale Eye 

Other studies on basins of about the same areal extent as Boot 
Creek (31 ee panticuleriy<those an the, high arctic, have reveateqetic 
spring snowmelt flood to be the major event of the year. At Weir River 
basin (area 29.4 mice on Ellesmere Island, the 1973 maximum unit area 
discharge was 0.17 Hea (Ambler, 1974). At the Mecham River basin 
(area 78 km), on Cornwallis Island, a maximum unit area discharge in 
1970 of 0.34 m?/s/km* was recorded during the spring runoff (McCann, 
Howarth and, Cogley, 1972.) ps1.09). 

On May 1, 1973, the snowpack water equivalent at the Inuvik 
airport snow course was 140 mm. By June 10, Boot Creek basin was almost 
completely free of snow, and the flood had receded to l m/s (3568S ))% 
Neglecting all water storage sources except the snow pack, and adding to 
that the 44 mm of precipitation which occurred from May 1 until June 10, 
there were approximately 185 mm of water available for runoff, 
evapotranspiration and change in storage, during the snowmelt period. 


Through June 10, 155 mm of runoff were measured, or 85% of the “available” 
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water. This high percentage is consistent with findings in other arctic 
watersheds underlain by permafrost (for example, Brown, Dingman and 
Lewellan, 1968, p. 13: concerning a small watershed on the Alaskan coastal 
plain) 

No measurements of peak snowmelt runoff at Peter Lake are available 
to date. However, the hydrographs illustrate that spring snowmelt produces 
the most significant runoff of the year (Figure 4). Maximum recorded 
discharge at. the basin 1s 2.4 aos (84ers) or 005 Te yatie a7 eae) 
on June 12,°19735, during the snowmelt flood secession. 

A secondary maximum of runoff commonly occurs at these basins in 
August (Figures 2, 3, 4). Runoff increases noticeably following storms of 
a magnitude that bring little or no runoff response in June and July. Compare, 
for example, the marked increase in runoff following the storms of August 8 
and 22, 1973 at Boot Creek basin, with the lack of response tothe storm of 
July 13, 1973 (Figure 2). In the-menath of August, there 1s normally a trend 
toward cooler and wetter weather in the Delta region, with a consequent 
decline in the rate of evapotranspiration. It is hypothesized that the increased 
runoff in August is primarily due to a rise in base flow. This is thought to 
result. from a greater ‘availability of soil moisture, which 1tselt isa és pone 
to the aecline in the rate of evapotranspiration. 

Elsewhere (Anderson and MacKay, 1973, p. 54), mention has been made 
of the meager response of these basins to storms, in the form of direct 
runoff. Ratios of ‘direct runoff to precipitation are typically 0.05 
or less. Analysis of midsummer storms is continuing. Plots of the 
logarithm of discharge versus time for storms at Boot Creek yield a 


recession limb with a changé in slope, the change likely reflecting a 
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shaft from base’ flow plus interflow (direct runoff negligible) to base 


flow only. These recessions can be defined mathematically as follows: 
Q. = QK (1) 
Qe7 mcrypt eye (2) 


Q VS Lhe dascharge at time t, oe as the discharge at time: tO, (kes 

bite recession. Constant, ,and t* 15 a decay constant. For the intertiow plus 
base flow portion of the recession, K has varied from 1.26 to 1.48, and ¢t~ 
FLOMROltO s0SenOurs,. Ol.) £00544 Frey nae for midsummer Boot Creek storms. 
For the ensuing base flow recession, K has been 1.05 to 1.10, and t*, 247 
co) 497 hours, or ‘8.0 sto’ 16 erie: At Glenn Creek, Alaska (basin area of 
1.8 ioe t* was observed to be of the same order of magnitude, varying 


between 10.9 and 42.4 ey ene (Dingman, 1973, p. 449). 


EVAPOTRANSP IRATION 

"Latent evaporation" is measured with 500 cm” black porous disc 
atmometers, in order to determine the "drying ability of the meteorological 
environment'! (Robertson and Holmes, 1958, p. 399). Analysis of atmometer 


data from several stations across Canada has revealed the relationship: 
PE = 0.086(LE) 


where PE is the potential evapotranspiration (mm) and LE is the latent 
evaporation (em”) (Roberston and Holmes, 1958, p. 406). Using this 
conversion factor, the potential evapotranspiration near Cynthia Lake in 
Boot Creek basin for the period June 19 to September 21, 1973 was 186 mn. 


Near the south end of Peter Lake, the potential evapotranspiration from 
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June 12 until September 23 totalled 238 mm. For the period June 1 to 
September 30, respective totals for the two basins were likely in the 
neighbourhood of 225 and 275 mm. 

No class A evaporation pan data are available for Inuvik for 
comparison. However, estimates of mean small lake evaporation have been 
made for Canada (Ferguson, O'Neill and Cork, 1970), and at Inuvik, the 
total for June to September inclusive is 175 mm. These estimates are 
thought to provide a "useful index of potential evapotranspiration" 
(Ferguson, O'Neill and Cork, p. 1631). Note that the mean small lake 


evaporation is only in fair agreement with the 1973 atmometer data. 


THE WATER BALANCE 

For the 12-month period October 1, 1972, to September 30, 1973, 
an estimate of the water balance at Boot Creek has been derived. 
Precipitation at the basin, as supplemented by data from Inuvik airport, 
totalled 285 mm. (equal to the mean for 1961-70 for Inuvik). Runoff, 
from the beginning of snowmelt in May, 1973 until September 21 amounted 
to 210 mm. Runoff was not measured in October, 1972, prior to freeze-up, 
or in late September, 1973, but the total for those missing periods would 
have been small by comparison, as runoff is of low magnitude late in the 
cig Es 

If the assumption of zero storage change is made for the 12-month 
interval, then actual evapotranspiration is at most 75 mm. This is 
probably too small a value. Snow sublimation could be at least 40 mm 
for the winter of 1972-73, and evapotranspiration for the summer of 1973 


might be 150 mm or more (letter dated March 5, 1974 from G. R. Holecek, Hydro- 
logy Branch, Department of Environment, Alberta, 10040 - 104 Street, Edmonton). On 
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this basis, the total runoff, snow sublimation and evapotranspiration 
would be at least 400 mm. Assuming the precipitation total to be correct, 
a storage depletion of the order of 115 mm would have to have occurred 

at. BooteCreck basa. 

It is possible, however, that the precipitation total may be 
inaccurate. Precipitation gauge deficiencies that cause undermeasurement 
are known to occur, particularly when high winds accompany the precipitation, 
and especially if the precipitation is in the form of snow. (See, for 
example, Rodda, 1971). The more important gauge ep events: turbulence 
and eddies caused by the gauge itself, can be overcome by putting the 
receptacle into a conical pit, in such a manner that the top of the 
instrument is at the same level as the ground surface. A comparison 
between the precipitation catch of a tipping bucket rain gauge located 
in a pit, and the catch of an identical gauge located on level ground 
nearby, was made at Peter Lake in the summer of 1973. The experiment 
revealed that over a 3-month period, the pit gauge caught 15% more rain- 
fall than its counterpart. It is conceivable that undermeasurement of 
precipitation of the same order of magnitude could be occurring at Inuvik 
the year round, especially when it is realized that more than halt of the 
annual precipitation normally falls as snow there. On the assumption of 
a 15% error in measurement, precipitation at Boot Creek basin for the 
year October 1,:1972, to September 30,1973 would be 328 mm, and storage 


depletion would then be 72 mm. 


CONCLUSIONS 
Precipitation in the vicinity of Inuvik is normally light, 


averaging about 285 mm annually. In 1973, the total was 20% below normal. 
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Peak runoff of 6.9 ne occurred during the snowmelt runoff 
period in 1973 at Boot Creek. This was by far the greatest TUNOLE Of 
the year. A secondary maximum occurred in August, as base flow incheased. 
Direct runoff response to midsummer storms continued to be minimal, 
constituting less than 5% of precipitation for major storms. Long 
recessions with large decay constants were noted for those storms. 
Potential evapotranspiration for the period June to September, 
1973 has been estimated at 225 to 275 mm for the basins studied. Water 
balance calculations for Boot Creek basin revealed a storage depletion 
of as much as 115 mn. for the period October 1, 1972 touSeptember 30, 2973. 
However, if precipitation was being undermeasured, this depletion would 


be considerably less. 


IMPLICATIONS FOR PIPELINE AND HIGHWAY DEVELOPMENT 

Exploration for oil and gas in the Mackenzie Delta region is 
continuing. If sufficient reserves are found, pipelines may be constructed 
to transport those commodities to southern markets. Construction has 
already begun on the Mackenzie and Dempster highways in the Delta region. 

In building these roads and pipelines, it is important that a 
proper assessment of the magnitude of hydrologic events be made. This is 
particularly so in the case of river crossings. Inadequate planning may 
result in spring wash-outs or winter icings, which would jeopardize the 
continuity of the transportation system. Analysis of data collected at 
Boot Creek and Peter Lake watersheds should prove of great value in 
predicting the magnitude of hydrologic events to be expected in the eastern 
portion of the Mackenzie Delta region. For example, the peak flow data 
collected during the 1973: snowmelt speriod are extremely. userul for the 


planning of the size of culverts to be used along the local highways. 
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PLANS FOR FUTURE WORK 

In addition to the program of data collection and analysis 
already underway, other projects are being planned. The spring snowmelt 
runoff peak has been identified as the most significant event of the 
annual hydrograph in the Mackenzie Delta region, and attention is to be 
focussed upon that event. Instrumentation has been installed in Boot 
Creek basin to study the ground and snowpack temperature regimes during 
snowmelt. Time and manpower permitting, measurement of the snowmelt 
flood will be extended to some larger drainage basins Brbeced by the Mackenzie 
Highway near Inuvik, in order to compare unit area discharges for drainage 
basins “Oi various: Sizes: 

A comparison of pit and regular tipping bucket rain gauges is 
to be continued at Peter Lake, and a similar experiment could be set up 
in Boot Creek basin. Representative values of gauge catch deficiency need 
to be established for these basins, in order to obtain better estimates of 
basin precipitation. 

Soil temperature regimes are to be monitored on a gravel spit in 
Peter Lake, and on north- and south-facing slopes at Boot Creek basin. 
Readings of water temperature at the discharge measurement location in each 


basin will also be taken. 
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Table 1 Annual precipitation normals for Inuvik and Tuktoyaktuk, 
for the peraod: 19o)-70- 


Inuvik Tuktoyaktuk 
(mm) (mm) 
Total precipitation Zoo Loy 
Rain 1235 80 
Snow (water equivalent) 162 oY, 
Snow (depth when newly fallen) 1840 [370 


Table 2. Monthly precipitation totals for Boot Creek basin and Peter Lake 
basin (where available)* 


Inuvik Airport Boot Creek and, an parenthesis, Peter Lake 


Normals (mm) 
1961-70 
(mm) 1970 1971 1972 1973 

JAN 19 Z 10 14 24 
FEB 1S S 14 12 4 
MAR 15 9 vi LS 10 
APR 15 10 9 28 9 
MAY 21 LS 3 is 36 
JUN 23 Bs 8 44 28 
JUL 40 28 Daj Pas) 12 42(34) 
AUG 54 59 40 50 (63) Soo) 
SEP 19 39 19 oy), 4 
OCT SZ 43 74 39 20 
NOV 15 15 19 wo) 14 
DEC Zi 18 14 36 6 
ANNUAL 285 281 244 309 230 
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ABSTRACT 


A hydrometric study is in progress in the Mackenzie 
Delta, N.W.T., to determine the discharge and seasonal distri- 
bution of flow through the delta, hydraulic geometry relation- 
ships in delta channels, water and ice levels, ice thicknesses, 
and the causes of delta flooding. A network of 28 stations to 
measure discharge has been established in the delta region. 
Results of the flow distribution program for three successive 
winters are presented. At the time of measurement (March), the 
Middle Channel carried most of the inflow to the delta. More 
than 80% remained in the Middle Channel as far downstream as a 
site west of Inuvik, and more than 55% at a site west of Reindeer 
Station. Below the latter section, flow was dispersed, into major 
distributaries branching away from the Middle Channel. Delta 
discharge (expressed as a percentage of total inflow to the 
delta) entered the Arctic Oceanvia the following routes: 
Shallow Bay -- 33 to 62%; Mackenzie Bay, east of Shallow Bay -- 
4 to 28%; Kugmallit Bay -- 25 to 35%. Ice thicknesses, obtained 
during the winter flow program, have varied from 0.7 to 1.6 m 
for individual cross sections. Mean ice thickness at all sites 


has been 1.0 to 1.2 m. 
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PROGRESS REPORT ON 
WINTER DISTRIBUTION OF FLOW IN 
THE MACKENZIE DELTA, N.W.T. 


by 
J.C. ANDERSON and R.J. ANDERSON 


INTRODUCTION 
In co-operation with Water Survey of Canada, a hydrometric 
Study is being conducted in the Mackenzie Delta, N.W.T., to determine: 
(1) the discharge and seasonal distribution of flow 
through the delta; 
(2) the hydraulic geometry equations which pertain to flow 
during the open water season; 
(3) variations in water levels, ice levels, and ice thicknesses 
on major delta channels; 
(4) the causes of flooding in the delta. 
Since publication of a preliminary report (Anderson and 
MacKay, 1973), research has been focused on the winter (March) dis- 
tribution of flow. Results of this program to date are presented 


in this report. 


THE FLOW MEASUREMENT NETWORK 

The network of discharge measuring stations is shown in 
Figure 1, and a list of station names appears as Appendix I. 
Table 1 outlines the format used to estimate flow at various 


sections across the delta. 


Se Aue 


The amount of water which enters the Mackenzie Delta in 
winter is considered to be the sum of contributions from the Peel, 
Arctic Red and Mackenzie Rivers. Discharge from smaller streams 
such as the Rat and Rengleng Rivers is assumed negligible during the 
winter season. Section U-U (Figure 1) represents delta inflow. 

Within the delta, discharge measurements are made on the 
main channels crossing various sections of the delta. Totals for the 
cross-sections V-V, W-W, and X-X + Y-Y are calculated as a check on 
the conservation of flow. Several smaller channels are omitted from 
the measurement program. Flow, if any, contributed by them across 
a section is estimated by subtracting measured output along the next 
downstream section from input above the section in question. 

Neglecting time lags, the flow across sections U-U, V-V, 
W-W and X-X + Y-Y at any one time should be approximately equal 
during the winter season, owing to the relatively static flow condi- 
tion of the input rivers. However, the length of the measurement 
program (usually about two weeks) introduces a time lag. Also, the 
severity of winter weather conditions can reduce the accuracy of 
discharge measurements made in winter over those taken during the 


summer. For these reasons, some discrepancies in flow can be expected 


EOuOCCIUT. 


THE MARCH FLOW DISTRIBUTION 


Three winters! data are now available for comparison. 


A full data listing is provided as Appendix II. Discharge data 


i eg 


for the three years are compared in Tables 2 and 3. 


Flow into Delta - Section U-U 

At the time of measurement, inflow to the delta was higher 
by OL seimeloyoeandaey 2 in 1974 ethane in 19722 ahs tlow variation 
is dependent upon conditions upstream in the Mackenzie basin, and 
no attempt is made here to explain it. As would be expected, the 
Mackenzie River is the largest by far of the three measured contri- 


buting sources, supplying between 97.0 and 98.4% of total flow. 


Flow at Station A 

Station A is on a delta channel that links the Peel Channel 
with a branch of the Mackenzie, near Point Separation. Reversals of 
flow are known to occur along this waterway. Discharge rates 
obtained by actual measurement at A are listed in Table 2. In 1972; 
no flow was recorded, but an overflow of water onto the ice surface 
occurred, suggesting a pressure build-up behind an ice blockage in 
the channel. If flow at A is calculated using the method of Table 1, 
negative values result (1973: -41 ae 1974: -32 oe This 


implies an inflow from the Mackenzie to the Peel Channel. 


Flow in Middle Channel at Vy 


—<—<$———$ $$ 


Discharge at this location in the Middle Channel is 
assumed to be equal to that not measured elsewhere along V-V. With 
only minor possible exception, this is true. Note that between 


96.0 and 97.5% of total inflow remains in the Middle Channel as 


Ber ACiC ieee 


far as Vy > at the time of measurement (Table 2). 


Flow Across Section V-V 
Very favourable results are observed in a comparison of 
flow across V-V with that across U-U. Over the three seasons, the 


maximum deviation at Section V-V from total inflow is only 0.4%. 


Flow at W, 
This represents the discharge which remains unmeasured 
across that portion of Section W-W lying westward of Y-Y. To obtain 
it, outflow along the West Channel (station 11) is subtracted from 
inflow of the Peel, Pokiak and Aklavik Channels (stations 9, 10 and 
12). The 1974 estimate, calculated as above, was -48.5 aie This 


result is assumed due to measurement error, and flow at W, in 1974 


has been set at 0 ee 


Flow Across Section W-W 

Discharge to the east of Y-Y across Section W-W is all 
accounted for in measurements, and flow across W-W follows as noted 
in Table 1. Comparison of total discharge between Sections V-V and 
W-W shows a loss of water in 1972 and 1973 of 15% and 8%, and 
a-daiveined 9/4 0f 57. 


Flow at Xy 


This is the residual of discharge not measured across Section 


X-X. It accounted for 3% of total inflow in 1973, and 4% in 1974. 


- 234 - 


Flow Across Section X-X 
Discharge along this section enters Shallow Bay from 
the west and south. The quantity has varied between 2 and 6% 


OF total int low. 


Flow at Y, 


This represents unmeasured flow across Y-Y. In 1972 
and 1974, it was 4 and 5% of total inflow, respectively. In 
1973, when inflow was greatest, the corresponding value for Y, 
was 26%. The smaller channels crossing Y-Y were obviously playing 
a more important role in flow distribution in 1973. Lower ice 
thicknesses (to be mentioned later) were probably an important 


cause of this occurrence, in addition to greater inflow. 


Flow Across Section Y-Y 

The Middle Channel is by far the major contributor of 
flow across this section. More than 80% of total inflow has 
been measured each year at station 18, and more than 55% at 
Station 20. 

Summation of flow across Sections X-X and Y-Y 
represents a complete flow section across the delta. Discharge 
at this point in the delta was found to be less than total 


inflow in 1972 and 1973, but more in 1974. 


Flow at Z, 


All main channels crossing ZZ are accounted for with 


= 92300. 


the exception of one (Z,). In 1973, when flow Y was significant, 
Z, was only 3% of total inflow. In 1972 and 1974, the reverse 
Situation existed, with Z, accounting for 19% and 27% of the 


total inflow. 


Flow Across Section Z-Z 

Most of the water crossing this section is supplied 
by the Middle Channel; flow in the East Channel at and above 
Station 21 is minor by comparison. Near Tununuk, the flow is 
distributed westward into Shallow Bay via Reindeer Channel, 
northwestward into Mackenzie Bay, and northeastward via the 
East Channel into Kugmallit Bay. 

Kugmallit Bay has received between 25 and 35% of total 
delta inflow during the periods of investigation in March. Discharge 
into Mackenzie Bay, eastward of Ellice Island, has accounted for 
20, 4 and 28% of the total inflow in three successive years. The 
remainder has entered Mackenzie Bay via Shallow Bay. Percentages 
of total inflow to the latter are: 1972 - 33%; 1973 - 62%; 

1974 - 43%. 


ICE. COVER 
Ice Thickness 

Ice thickness data obtained in conjunction with discharge 
measurements are listed in Table 4. As would be expected, thickest 
ice occurs at the more northerly stations. A major reason for 
this is that the snowfall and snow cover decrease northward toward 


the arctic coast. Snowfall normals for Ft. McPherson, Inuvik and 


Se igo) 


Tuktoyaktuk illustrate this (Table 5). A heavy snow cover insulates 
the ice below it from extremely low air temperatures, and consequently 
ice accretion takes place at a slower rate. 

Other important parameters that affect ice depth include 
channel width and the orientation of channels with respect to 
prevailing winds. Narrow channels are better protected from the 
wind by the proximity of their banks. This reduces the amount of 
snow redistribution, and a thick snow cover can develop. Such 
is the case at channel Sections 4 and 5 in the south end of the 
delta, and at Sections 17 and 19 near Inuvik, among others. Wider 
channels, where they are oriented parallel to the prevailing wind, 
can be swept clean of snow along some sections. Broad stretches 
of snow-free ice were observed on the East Channel below Reindeer 
Depot and on the West Channel below Aklavik, in 1974. 

Ice was generally thicker in 1972 and 1974 than in 1973. 
Temperatures were only slightly below normal through the winter 
prior to March 1972 (Table 5), but snowfall was less than normal 
at Ft. McPherson. At Inuvik, a snowfall anomally of 118 cm 
occurred in October, 1971. Much of it fell before the ice formed 
on delta channels. Thus, the October snowfall did not affect ice 
accretion later in the winter. Snowfall was close to normal 
for the remainder of the winter at Inuvik. 

The winter of 1972-1973 was above normal temperature-wise 
in the delta region, and above normal for snowfall at Inuvik and 


Tuktoyaktuk. As a result, ice was thinner, the mean depth being 


as ON teas 


0.2 m less than for the previous winter. 

Ice thicknesses in 1974 were similar to those of March, 
1972. This was apparently a result of low snowfall during previous 
months. Ft. McPherson, in particular, received only 42% of its 


normal fall for the period September to February. 


Other Ice Properties 


Overflows of water onto the ice surface, which happened 
in March, 1972, and which were caused by blockage or restriction 
of flow due to ice, did not occur in the next two seasons. No 
slush or frazil ice was found beneath the ice cover in the delta, 
and no polynas were observed. 


Rough ice was encountered on some larger channels. These 


were the last to freeze over in the fall, and they became the resting 


place for broken ice from upstream. Greatest difficulty with 
rough ice was encountered on the Middle Channel. Near Horseshoe 
Bend (site C), ice protruded at all angles to heights of two and 


three feet, and sometimes more. 


CONCLUSIONS 
Based upon 3 years of March flow data, the following 
conclusions are made: 
(1) Flow entering the delta in March can vary considerably 
in magnitude. This is dependent upon upstream conditions. 
The Mackenzie River is the most important source, 


contributing 97.0 to 98.4% of total flow. 


SS os 


(2) 


The Middle Channel of the delta carries most of the 
winter flow. More than 80% of total inflow has been 
recorded as far downstream as station 18, each year. 
At station 20, near Reindeer Station, percentages 

of total inflow for the years 1972, 1973 and 1974 
were 78, 57 and 91, respectively. 

A small amount of water appears to be entering the 
Peel Channel from the Middle Channel near Point 
Separation. This was recorded in two of the three 
years, as was flow down the East Channel at stations 
B and 16. 

Outflow to the Arctic Ocean from the Mackenzie Delta 
is distributed in the following manner (expressed as 
percent of total inflow): Shallow Bay -- 33 to 62%; 
Mackenzie Bay east of Ellice Island to Richards Island 
Bay -- 25 to 35%. 

Ice thickness in the delta has varied between 0.7 and 
1.6 m. Mean ice depth at all sites has been between 
TeGsandel 20m: 

Overflows of water onto the ice are not common during 
drilling in March. 

No slush or frazil ice has been encountered beneath the 
ice cover. Rough ice is to be expected on the larger 


channels, especially the Middle Channel. 
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Table 2. March distribution of flow in the Mackenzie 
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Table 4. March ice thickness data, Mackenzie Delta, N.W.T. 


Area Station Mean Ice Thickness (m) 

1972 1973 1974 
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Appendix I Station Locations, Mackenzie Delta, N.W.T. 


Peel River above Fort McPherson - gauging station 
Arctic Red River at Martin House - gauging station 
Mackenzie River above Arctic Red River - gauging station 


Husky Channel 

Phillips Channel 

Peel Channel 

Kalinek Channel 

East Channel (Mackenzie River) 

Peel Channel above Aklavik 

Pokiak Channel above Aklavik 

West Channel (Mackenzie River) below Aklavik Channel 
Aklavik Channel above Schooner Channel 

Nikoluk Channel 

Leland Channel 

Hvatum Channel 

East Channel (Mackenzie River) at Inuvik - gauging station 
Kalinek Channel - N.W. of Inuvik 

Middle Channel (Mackenzie River) above Napoiak Channel 
Napoiak Channel 

Middle Channel (Mackenzie River) west of Williams Island 
East Channel (Mackenzie River) below Williams Island 


Marcus Channel 


Sc pie 


Appendix I Continued 


23. Reindeer Channel 
24. East Channel (Mackenzie River) below Tununuk Point 


25. Channel to west of, and below Tununuk Point 


Ae Channel below Indian Village 
B. East Channel (Mackenzie River) above Kalinek Channel 


C. Middle Channel (Mackenzie River) above Horseshoe Bend 
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PROGRESS REPORT ON 
GEOMORPHIC AND CLIMATIC STUDIES IN THE MACKENZIE DELTA AREA 


AIDED BY DENDROCHRONOLOGY 
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Glaciology Division 
Water Resources Branch 
Department of the Environment 


under the 


Environmental-Social Program 
Northern Pipelines 
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During the month of August 1973, a field party consisting of 
W.E.S. Henoch, M. L. Parker (Western Forest Products Laboratory, 
Vancouver, B.C.) and D. Outhet (graduate student, University of 
Alberta) visited the Mackenzie River Delta to study some aspects 
of geomorphic processes and climate in the region. Extensive tree 
sampling was carried out during the field season, particularly in 
the vicinity of Point Separation where sections of buried tree stumps 
were collected to study: 
(a) growth of adventitious roots and the rate of sedimentation; 
(b) age of the holocene deposits in the Delta near Point 
Separation; and 
(c) minimum age of a buried tree stump pierced by an ice wedge. 
In addition, three pingo-like humps were located between drained channels 
on the Peel River floodplain near a site known as Indian Village. There, 
preliminary work was carried out as follows: 
(a) five vegetation zones were recognized and mapped. 
(b) flood levels were surveyed; 
(c) the depth of active layer across the humps was surveyed 
by probing at one metre intervals; 
(d) the active layer depth was also probed under a windthrow. 
Tree samples were taken to determine the age of the windthrow 
and the rate of degradation of the permafrost under it; and 
(e) tree samples were collected in order to study the relation 
between tree rings and permafrost-soil conditions. 
Two new tree ring chronologies were built under Parker's supervision 


in the Western Forest Products Laborary from the samples collected in 1972. 


ae AON 


These chronologies of ring width and density in the Inuvik and 
Fort Simpson areas were compared with recipitation, temperature 
and river discharge records for the growing months from May to 
September. Some highly significant correlations were obtained. 
Analytical work is in progress and the results will be reported 


in the near future. 
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ABSTRACT 


In 1973, hydrologic data collection was continued in the Twisty 
Creek watershed. This small sub-arctic upland basin in the Mackenzie 
Mountains was chosen to examine the relationships among precipitation, 
runoff and sediment yield. Runoff from a heavy rainstorm on July 6th 
resulted in a new peak flow of 16.99 rye or 2.45 mein The percentage 
of runoff to precipitation was found to average 91 per cent over the field 
season with a resulting low evapotranspiration component. About 80 per cent 
of total sediment yield occurred during peak flood events. Such results 
should be of value in developing design criteria for highway and pipeline 


crossing of rivers and streams. 


= 202) 5 


HYDROLOGIC STUDIES AT "TWISTY CREEK" 
al 
MACKENZIE MOUNTAINS, N.W.T. 


- 1973 Progress Report - 
By 
J.N. Jasper 


"TWISTY CREEK" WATERSHED STUDY 

"Twisty Creek" is a small sub-tributary of the Arctic 
Red River and is located in the foothills of the Mackenzie Mountains, 
N.W.T. at latitude 65°23' N and longitude 131°16'W (Figure 1). 
The study was initiated in 1972 to understand the hydrologic response 
characteristics for small drainage basins in the area (Sellars, 1973). 
This location was selected with the intention of filling a gap in the 
present network of meteorological stations as well as for the size, 
accessibility and representativeness of the watershed. In addition, 
the Mackenzie Mountains have been identified as the source area for 
major floods in the MacKenzie River tributaries on the west side of 


the Mackenzie Valley (Mackay, Fogarasi, and Spitzer, 1973). 


The specific objectives are to collect and publish relevant 
data, and derive the relationships among precipitation, runoff, and 
erosion for a small upland watershed. Results are presented for 1972 and 
1973 and may be of use in determining design criteria for river crossing 
structures in watersheds with similar physiographic and precipitation 


characteristics. 


=F 203. = 


135° 130° 125° 120° << 


70° 


e 
© Herschel Island 


MACKENZIE VALLEY REGION SCALE 


50 100 miles 
CLIMATIC STATIONS e 
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Figure 1. Mackenzie Valley region snowing tne location of Twisty Creek 
and other meteorologic stations. 
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METHODS OF DATA COLLECTION 

As in 1972, six M.S.C. pattern tipping bucket rain 
gauges and three Lambrecht model 252C thermohygrographs in 
Stevenson screens were utilized to record units of 1/100 of an 
inch of rainfall, air temperature and humidity. A continuous 
record of stream stage was kept using a Stevens type F water level 
recorder at the main discharge station. Stream velocities were 
measured with an Ott current meter and rating curves developed 
from 12 measurements in 1972 and 32 in 1973. 

Dissolved load in the stream has been estimated from 
water chemistry analyses of stream waters, which were tested for 
pH, calcium and total hardness, and alkalinity. Dissolved loads 
were roughly calculated as the sum of the latter two, with the 
assumption that the major proportion of dissolved solids are of 
this composition. 

A U.S. DH-48 suspended sediment sampler was used to 
collect depth integrated suspended sediment samples at the main 
stream gauging site. Bed load moving in the stream at this site 
was sampled with a 43-inch mesh basket-type sampler. Discharge- 
sediment rating curves have been developed from 191 water chemistry 
analyses and 60 suspended sediment samples taken during the two 
field seasons. The 12 bed load measurements from 1973 have been 


utilized to give a preliminary rating curve for bed load. 


ee OO = 


The records obtained are for the period July 2 - August 25, 
1972 and June 29 - September 5, 1973. The 1973 network was basically 
the same as that established in 1972, however the rain gauge at 
site Rl was moved to R7 and a standard gauge substituted at Rl 
(Figure 2). Water level recorder W2 was moved to site W3 to allow 
development of a unit hydrograph for a smaller unit basin. However 
problems with stilling well intakes at both W2 and W3 have resulted 
in discontinuous records which have not been analyzed. 

A ground control survey carried out in 1973 has resulted 
in the drafting of a 1:10,000 scale contour map from 1972 aerial 
photographs. Drainage area of the basin above the main stream gauging 
site is 6.94km*, 6 percent greater than the original estimate of 
6.55km*, which was measured from an enlargement of the available 


1:250,000 NTS coverage (Figure 2). 


RESULTS AND DISCUSSION 
Precipitation and Runoff 

Tables 1 and 2 are summaries of precipitation, runoff, 
and sediment yields compiled on a daily basis. From these records 
the average precipitation over the catchment for the periods of 
record has been computed to be 234mm in 1972 and 367mm in 1973. 
Total runoff was 162mm or 69% of precipitation in 1972, 332mm or 
91% in 1973. Both ratios are high and indicate small evapotrans- 
piration losses, as would be expected from a watershed with less 


than 20% vegetation cover. As well, the soil in the basin is 
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LEGEND 131°15'00" 
STEVENSON SCREEN 
WATER-LEVEL RECORDER 
RAIN GAUGE 
WATERSHED ABOVE WI........ 
CONTOUR INTERVAL 


65°22 30° 


65°22'30" 


TWISTY CREEK 
DRAINAGE BASIN 
N.W.T. 


| kilometre 


Figure 2. Location of instrumentsjn "Twisty Creek" drainage 
basin 
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quite shallow except for a narrow strip immediately adjacent to 
the main channel, 

Base flow separation has been carried out for both years' 
data from observations on the response of groundwater flow to 
short intense rainfalls, as well as long duration light intensity 
events. Figure 3 gives several examples from 1973. A short 
duration high intensity rainstorm typically produces a direct runoff 
peak, which occurs as rainfall ceases. Groundwater flow begins to 
increase 6 to 8h later and produces a second peak some 12-15h 
after rainfall has stopped. Sellars (1973, p. 527) has presented 
Similar evidence for 1972. Long duration or light intensity storms 
produce a similar rise in groundwater, usually without the direct 
runoff peak. These response characteristics of groundwater flow, 
and the linear recession curves of groundwater flow obtained from 
semi-logarithmic plots of the streamflow data, have allowed a 
physically based graphical separation of direct runoff and baseflow. 

Further analysis of the ratio of runoff to precipitation 
has shown that the high value of 91% for 1973 can be largely 
explained in terms of the rapid runoff which follows high intensity 
storms. The most impressive example is for the maximum flood 
experienced over the 2 years of record, when 57.78mm of rain 
(44.68mm in a 2h period) produced a peak discharge of 16.99m3?/sec on 
6 July 1973 (Figure 4). At the time, baseflow was very low (0.023m3/sec) 
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due to a lengthy dry period; thus the major proportion of runoff 
in the several days afterwards was largely a result of this storm. 
The hydrograph separation technique described above,results in a 
prediction of 41.58mm for direct runoff and interflow, which 
corresponds to a direct runoff ratio of 0.72. A second major 
rainfall of 60.19mm on 4 August 1973 produced a peak discharge of 
8.14m2/sec and 36.53mm of direct runoff + interflow. 

If these two events had not occurred the runoff- 
precipitation ratio for 1973 decreases to 0.70, which is 
comparable to the 1972 value of 0.69. Therefore the ratio of 0.91 
no longer seems as unrealistic as it first appeared. The 
estimates of total direct runoff + interflow based on hydrograph 
separation for the two field seasons are 51mm in 1972 and 120mm 
in 1973, which correspond to 31.6 and 36.2% of total streamflow, 


21.8 and 32.9% of precipitation. 


Sediment Yield 

The results of the water chemistry analyses are presented 
in Figure 5, with regression lines fitted by eye. Although fewer 
data are available for 1973, the pattern of changing dissolved 
solids concentration under different flow conditions noted during 
1972 (Jasper, 1973, p. 654) is repeated. Transition from one 
rating curve to the next coincided with changes in flow conditions, 


often resulting from major floods, and is attributed to the 
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relative contribution of groundwater flow to runoff. 

Suspended sediment data (Figure 6) show good consistency between 
the two field seasons. As 1973 was a much wetter year (see Tables | 
and 2), with an estimated 367mm of rain compared to 234mm in. | Sies 
a greater proportion of samples are for higher discharges than LN OW Ze 
Also the 1973 suspended sediment concentrations are higher for similar 
discharges, possibly due to higher moisture conditions prevalent during 
the 1973 field season. Simple regression between discharge and suspended 
sediment yields a correlation coefficient of 0.88. It is hoped that the 
77% of variance in suspended sediment concentration explained by variation 
of discharge may be increased, perhaps by representing moisture conditions 
through an antecedent precipitation index. 

The 12 bed load samples have been analyzed and total yields calculated 
for the stream cross-section. The results have been interpreted in a 
somewhat arbitrary manner for the following reasons. Firstly, it has been 
assumed that the traps were 100% efficient. This is of course not in 
agreement with the laboratory tests of Gibbs and Neill (1973) and others 
too numerous to mention. Moreover, the flow conditions were so turbulent 
during sampling that the traps had little effect on movement of bed load. 
In addition, a least squares fit to all 12 data points gives an estimated yiel 
rate of 1,080,000kg/sec at the peak runoff rate of bee oom/ sac on 6 July 19738 


As the runoff rate corresponds to 17,000 kg/sec of water, this estimated 
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value is unrealistic. Alternatively, separate regression lines were fitted to 
all data points as well as to the 4 highest points. The resulting correlation 
coefficients were 0.91 and 0.95 respectively. The upper line (Figure 7) 
estimates a more realistic transport rate of 1,030 kg/sec at 16.99 m>/sec and 
has been used when discharges exceed 3.5 m>/sec. 
Computed sediment yields for 19/72 are as follows: bed 
load - 0.51 million kg; suspended load - 0.56 million kg; dissolved 
load - 0.22 million kg; for a total loss of 1.29 million kg of 
material. The values for 1973 are: bed load - 4.87 million kg; 
suspended load - 4.57 million kg; dissolved load - 0.45 million kg; 
for a total loss of 9.86 million kg. Assuming an average rock 
density of 2.71, these totals correspond to basin erosion rates of 
0.0069 cm/yr for 1972 and 0.0524 cm/yr for 1973. Both values are 
within an order of magnitude of the estimated range of 0.0399 - 
0.0998 cm/yr quoted for high mountain areas by McPherson Gey alias 
p. 232). In addition McPherson obtained a value of 0.0495 cm/yr of 
erosion for Two O'Clock Creek, Rocky Mountains. 
The significance of the two major floods of 1973 and one 
in 1972 may be appreciated when one considers that 4.59 million kg 
bed load and 3.10 million kg suspended load were transported past 


the gauging station in a total of 16 h in 1973. Almost all bed load 
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movement and 0.51 million kg of suspended load occurred during an 8 h period 
in 1972. Daily values of dissolved material transported varied over a smaller 
range due to the decrease of dissolved solid concentrations with increasing 
streamflow. Thus approximately 80 percent of total erosion occurred in time 
intervals accounting for 0.6% of the 1972 field season and 0.8% in 1973. 

In addition to the problems of bank erosion and channel 
shifting due to a highly mobile bed during peak floods, care must 
also be taken to provide protection for structures in the stream 
channel and to ensure passage of high bed load. As an example, 
a concrete and boulder weir constructed at the mainstream gauging 
site in late June 1973 was largely destroyed during the July 6 flood 
by boulders being swept along the bed and into the upper face of 
the weir. The weir was reconstructed after the flood subsided but 
continued to be slowly worn away as the rest of the field season 
passed. It is recommended that special care be taken in the 
design of structures crossing streams with similar bed load movement 
patterns. The size of bed material, channel slopes, watershed 


relief, etc., are obvious reference indicators. 


CONCLUSIONS AND RECOMMENDATIONS 

Although two summers of data do not warrant the drawing 
of definite conclusions, several hydrologic and sediment movement 
characteristics have been observed which are thought to be valid 


for similar watersheds in the region. The observations and 
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and tentative conclusions are as follows: 
1. Mean rainfall was 4.25mm/day in 1972 and 4.65mm/day in 
1973. 

2. Mean total runoff was 2.95mm/day in 1972 and 4.20mm/day 
in’ 1973. 

3. Major floods result from heavy rainstorms produced 
orographically by systems approaching from the NW, N 
and NE (Sellars, 1973). 

4. Runoff of excess water is rapid, with lag time from 
centroid of the mass of effective rainfall to centre of 
area of the unit hydrograph computed to be 2.5h 
(Sellars, 1973). 

Ds Groundwater flow, as estimated from a physically-based 
graphical separation technique, accounts for 60-/0 per- 
cent of total runoff. Groundwater responds more 
slowly to rainfall, beginning to rise 6-8h after rainfall 
begins and peaking about 12-15h after rain ceases. 

6. The peak discharge of 16.99m?/sec (2.45m?/sec/km* or 
229 cfs/miz) is thought to be the highest since at least 
June 1972 and possibly for several years previous, due 
to lack of similar evidence of recent major bank and 


gully erosion at the commencement of the study. 


7. As much as 80% of total sediment yield may occur during 


the peak flood event(s) in a year. Extremely high rates 


= 2/6 i= 


Duration 


of bed load movement (i.e. 1030kg/sec or 141kg/km* on 

6 July 1973) can be destructive to improperly designed 
structures in the stream channel. Consideration for 
passage of high bed load should be incorporated into 
their design. 

Peak at-a-station and basin averaged rainfall for several 
durations, based on the available data, are presented 
below. The 25-year return period regional intensity- 
frequency curves for the area have been described by 
Northwest Hydraulic Consultants Ltd. (1972) and are given 
for reference. It should be noted that the latter are 
obtained from records at Yellowknife, Whitehorse, Fort 
McPherson, Inuvik, and Aklavik, all of which are located 


at low elevations. 


Mackenzie Region At-a-Station Watershed Averaged 


25-Year Return Period Precipitation Precipitation 
mm mm mm 
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TABLE 1. DAILY RAINFALL, RUNOFF AND SEDIMENT YIELDS FOR 1972 FIELD SEASON 
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ABSTRACT 


This study examines the behavior of oil under ice. 
Cold room experiments are carried out to determine the physical 
properties of crude oil in ice-covered water bodies. 


Comparisons are made with observations taken under 
field conditions. 


PROGRESS REPORT ON-A STUDY 
OF OIL POLLUTION IN ICE-COVERED RIVERS 


by 


Benjamin E. Keevil and René 0. Ramseier 


INTRODUCTION 


The discovery of oil in commercial quantities in northern 
Canada could lead to the construction of oi] pipelines to southern 
markets. Such pipelines would cross many rivers and streams that 
are subject to ice conditions of varying seasonal length. Any 
river crossing is hazardous from the oil pollution standpoint and 
a pipeline break at such a location could cause a great amount of 
environmental damage. Problems associated with containment and 
clean-up could also be intensified by a spill at a river crossing 
site, particularly if ice conditions prevailed at the time of 
the event. 


Oil spills can be large such as those caused by a major 
pipeline fracture or they can be small and insidious such as those 
developing through pinholes of about 1 cm diameter in the pipe. 

The behavior of oi] spilled on and under ice at particular stages 

of river ice cover development is not well understood. This study 
attempts to define the behavior of oil under ice and its interactions 
with an ice-covered river. 


To date, laboratory experiments to determine the behavior 
of oil under ice conditions have mainly been carried out in small 
0.3 m square tanks of 0.6 m depths topped by a cold plate connected 
to a refrigeration unit. In this study, a larger tank measuring 
1.5 m in diameter and 0.7 m in depth was set up in a cold room. 
Although this technique does not duplicate spill conditions in a 
natural ice-covered water body, it nevertheless provides considerable 
insight into the behavior of crude oi] under ice. In addition, 

a comparison of the experimental results with those of real spills 
adds to the understanding of these events. The objectives of the 
study are, therefore, to: (1) carry out cold room experiments 

to determine the physical properties of crude oil under ice 
conditions; (2) correlate the behavior of crude 011 in cold room 
tests with accidental spills in the field; and (3) write a scenario 
on the behavior and fate of spilled crude oil in an ice-covered 
river. 


METHOD AND INSTRUMENTS 


A circular aluminum basin, 1.5 m diameter and 0.7 m 
deep is installed in a cold room capable of maintaining the air 


Ya 


temperature at -15 + 3°C. The cold room is equipped with a fan 
to ensure uniform air flow and heat transfer in the basin. A 
pressure relief overflow pipe releases and measures the volume of 
water displaced by the growing ice sheet. Heating tapes are glued 
to the outside of the tank and are covered with 10 cm of polystyrene 
insulation. The heating tapes are connected to variacs set at very 
low power to cancel out the bottom and side heat losses. This 
combination of heating tapes and insulation allows one-dimensional 
freezing in the basin and effectively simulates the growth of 
natural ice. A schematic diagram is shown in Figure 1. 


The ice thickness is measured with ice-thickness gauges 
(Syrinkov, 1963). The gauge is a small steel cylinder suspended 
by a wire under the ice. To measure the thickness ,the cylinder 
is raised to the ice-water interface and the distance from the 
ice-air interface measured and subtracted from the total wire length. 
The oi] injection system is a 25-litre plastic tank with a gravity 
fed electric pump. The system is portable and can easily be carried 
into the cold room and connected to the oil injection pipe. 


An air-ice-water thermistor probe is used to obtain 
temperature profiles with a relative accuracy of + 0.2°C. The probe 
contains 19precision thermistors (YSI No. 44033) mounted in a PVC rod 
at a separation of 1 cm in the ice and 5 cm in the air and water. Each 
thermistor, solder joint and wire is waterproofed with epoxy to ensure 
that no water can leak into the system. A multi-conductor plug is 
installed at the top of the probe for connection of the bridge circuit. 
A rotary.switch enables each thermistor to be switched into the bridge, 
balanced and the resistance recorded. 


The typical test procedure is as follows. At a cold 
room temperature of -15-C the basin is filled with water toa 
height of 60 cm. Two days later, a primary ice layer begins to 
form. Its growth is predominantly in the horizontal plane forming 
a very thin ice layer. The secondary ice grows parallel to the heat 
flow, perpendicular to the water surface and reaches the desired 
(6 to 8 cm) thickness in three days time. The temperature gradient 
in the ice (about 0.5°C/cm) is linear. Norman Wells crude oil 
(5 litres at 20°C) is pumped for about 70 seconds through a 1.3 cm 
pipe and released 30 cm below the ice-water interface in the center of 
the basin. A plastic grid (2 cm reference lines) and clock are placed 
on top of the clear ice and the oi] behavior is photographed through 
the ice. Clear ice is obtained in the cold room due to one-dimensional 
heat flow in the basin. After the oi] is injected the ice continues to 
grow and "sandwiches" the oj] between ice layers. The mechanism of oi] 
entrapment is shown in Figure 2. 
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The following variables are kept constant for each 
test. 


oil type (Norman Wells crude, pour point -51°C) 

oil injection temperature (+20 °C) 

cold room temperature (-15°C) 2 

force and injection rate of oil pump (0.4 kg/cm’; 
/0ce/ Sec) 

ice roughness (smooth) 


PHYSICAL PROPERTIES OF NORMAN WELLS CRUDE OIL 


The most significant properties of crude oil for the 
under-ice experiments are "pour point" and specific gravity. 
Pour point is the temperature below which the oil does not flow 
freely. Norman Wells crude oil has a pour point of -519C. As 
all temperatures under jhe ice are above 0°C and the pipeline oil 
temperature is about 60°C, the crude oil will not freeze or 
become semi-solid. The specific gravity of Norman Wells crude 
(0.8) is less than water or fresh-water ice and as evaporation, 
aging and emulsification are negligible in still water under 
ice, crude oi] will tend to flow over rather than under the ice. 
However, a natural ice cover floats with its ice surface above 
the water level thus preventing oil from spreading over the ice. 


OIL BEHAVIOR UNDER ICE 


The circular basin tests are designed to investigate the basic 
behavior of crude oil under fresh-water ice. Observations of oi] particle 
formation, oil layer thickness, spreading rate under ice, sandwich mechanism 
and effect of a circular current have been made. 


Particles 


The most significant observation (Figure 5) of the 
behavior of hot crude oil released under ice is the separation of 
the oi] into hundreds of small globs or particles 0.1 to 2.0 cm in 
diameter. When released from the injection pipe, the oil acts 
much in the same way as water falling from a tap and is separated 
into particles by surface tension. The crude oil particles are 
less dense than water and will rise to the water surface. However, 
the ice sheet is a boundary and the particles are flattened on 
the under-ice surface by reason of their buoyancy. While at the 
ice-water interface the smaller particles coalesce to form larger 
particles. This coalescing continues until several minutes later 
an oil layer has formed, 
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0i1 Layer Thickness 


The oil at the ice-water interface forms a layer about 
1.0 cm thick under smooth planar ice. 


Spreading Rate 


The spreading rate of oil under ice is affected by the 
oi] injection rate, properties, partitioning, separation into 
particles and the ice roughness. In the circular basin tests the 
oil injection rate, properties and ice roughness were constant, 
however, the separation into particles and partitioning cannot be 
controlled. As a result the spreading mechanism is complicated 
and there is no simple equation to predict spreading rate. The 
graph (Figure 3) shows the oil spreading radially at a rate of 
about 1.0 cm/s. The contrast between oi] spreading on and under 
ice is shown in Figure 4. 


Sandwich Mechanism 


The mechanism by which an oi] layer becomes sandwiched 
between fresh-water ice layers has been established. Ice grows 
around the oil layer and the oi] is incorporated as a lense into 
the ice-oil-ice sheet. Even with a slow growth rate (0.5 cm/day) 
the oi] is not pushed ahead by the growing ice-water interface. 
During one test 5.0 litre of oi] was sandwiched between two 10 cm 
ice layers. A small hole was drilled down to the oil layer and 
about 1.0 litre of oil rose up the hole and spread on top of the 
ice. This indicates that sandwiched oil is pressurized within an 
ice cover which could be caused by cooling and subsequent contracting 
of the ice surrounding the oil layer. Future experiments will examine 
this mechanism further as it has practical application for clean-up 
of oil under ice. 


Circular Current 


To study the effect of a circular current on the behavior 
of oi] under ice an electric outboard motor was mounted on the 
Side of the basin and a circular 011 boom installed to prevent the 
oil from rising up the tank walls and spreading over the ice 
surface. The outboard motor was turned on when the basin was 
filled with water and left running for the duration of the test. 

A detailed velocity profile across the basin was not obtained, 
however, a small current meter indicated an average velocity of 
about 0.1 m/s. The oil] when released rose to the ice-water inter- 
face and circled around the center of the basin (Figure 6), Globs 
of oil broke up and coalesced due to the rotating motion. The 
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0i1 did not adhere to the smooth under-ice surface. The 

Berti tionimgsenrect. where the’ ligntey and heavier fractions: 07 
oil separate, was observed. During one circular current test 
freezing of slush ice produced a very rough and uneven ice-water 
interface. Most of the oil after injection collected in holes 
and pockets and one day later all the oil] was sandwiched by the 
growing ice (2.0 cm/day). 


OIL BEHAVIOR REPORTED IN ACCIDENTAL SPILLS 


To compare the cold room behavior of crude 011 with oil 
behavior at accidental spills a summary of oi] spills in rivers 
is being prepared. Unfortunately accurate field observations of 
oil spills in rivers are rare and it is difficult to find useful 
reports. Most reports just mention that spilled oil disperses 
rapidly in fast flowing waters. During 1972 and 1973 there have 
been about 20 reported oil spills in rivers but only a few in 
ice-covered rivers. 


TABLE 1. SUMMARY OF OIL SPILLS IN RIVERS - JAN. 1972 - AUG. 1973 


CAUSE NO. OF INCIDENTS GALLONS OF SPILLED OIL AT EACH INCIDENT 
Pipelines 3 20,000 52,000 250,000" 
Storage Tanks a 800* ‘ioe O00 

Tank Trucks 2 1200 7,000 

Tests ] 270 

Others 2 900 1,500 

Unknown 4 Minor 


* 701) Spi lle 1ee-covered river 


The only major reported spill in an ice-covered river 
occurred in Sweden. In February 1972, an estimated 250,000 gallons 
of diesel oi] was spilled under the frozen river Ume. (Jerbo, 

1973). The oil spread underneath the ice at the ice-water inter- 
face and collected in the hollow surfaces beneath the ice. When 

the oi] was sandwiched, the thickness was estimated to be about 

U.5 cm (personal communication, Jerbo). In the initial stages of the 
011 spill 100,000 gallons were contained by the ice and did not move 
more than 7 km downstream. 
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Two conclusions can be drawn from this summary of oil 
Suilis i ierivers. First, -o1l@spiiled an ice free yiyers is quickly 
dispersed and emulsified downstream. Second, oi] spilled in ice- 
covered rivers quickly disappears from sight under the ice and 
many spills are not reported. 


SUMMARY OF PROJECT PROGRESS 


Cold room experiments designed to simulate an under-ice 
pipeline leak demonstrated the basic behavior of crude oil under 
ice. The following observations were made: 


(1) That crude oil when released under the ice separates 
into small globs or particles 0.1 to 2.0 cm in diameter. 


(2) The particles coalesce at the ice-water interface and 
form an oil layer about 1.0 cm thick under smooth fresh-water ice. 


(3) In still water the oil layer is not pushed ahead of the 
growing ice-water interface but becomes sandwiched by the growth 
of the ice around the oil. 


(4) In moving water the oil can be sandwiched in pockets 
and holes at the ice-water interface providing the current is 
FeSsethannw= 1 m/s: 


Observations at accidental 011 spills have confirmed 
the basic behavior of oil under ice. 


PLANS FOR FUTURE WORK 


Future cold room experiments will determine the 
insulating effect of a sandwiched oi] layer on heat transfer 
through an ice sheet. The circular basin will be divided in two 
halves with polystyrene, one-half will have oil under ice and the 
other will be a control. 


Data on the fate and behavior of spilled oil at 
accidental oi] spills will be collected. Arrangements have been 
made with Environment Protection Service to notify us of any oil 
spills in ice-covered rivers. Sampling equipment to measure oi] 
layer thickness and to sample for oil in water under ice is being 
developed. 


This combination of cold room experiments and 
accidental 0i1 spill documentation will provide enough data to 
write a scenario on the behavior and fate of spilled crude oil in 
an ice-covered river. 
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ABSTRACT 


Analysis of time lapse photography in the field and from the air, 
along with other data collected in the field, indicates that in the 
southern Mackenzie River Delta, the shape of an eroding bank is positively 
correlated with the erosion process and the rate and character of erosion. 
The primary factors of current velocity, channel orientation to the wind, 
and ice content of the bank sediment act along with several minor factors 
to cause the differences in shapes. There are five different easily-dis- 
tinguished bank shapes in the study area each with its own maximum and 
minimum erosion rates and manner of erosion. This information allows the 
prediction of eroding bank behaviour and the production of a map showing 
the behaviour category into which each bank fits. This map may be used 
in the planning of construction in the area to avoid rapidly eroding banks 


such as those that may erode up to 30 m/yr. 
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PROGRESS REPORT ON BANK EROSION STUDIES 
IN THE MACKENZIE RIVER DELTA, N.W.T.* 
by 
DAVID N. OUTHET 


Introduction 


Human activities in the Mackenzie River Delta have increased markedly 
over the last few years due to the Search for oil and gas. Construction of 
wharves, buildings, drilling platforms and storage tanks has occurred or is 
contemplated at a number of sites within the delta. The choice of these 
Sites is important because delta flooding is common, bank erosion is variable 
and permafrost is present. This study examines the bank erosion problem and 
its relation to current velocity, wind, break-up associated current velocities, 
permafrost and other parameters. 

Bank erosion can be examined in a number of ways. Air photographs taken 
over a large interval of years can provide information from which estimates of 
erosion rates can be made. Prediction of bank behaviour may also be possible 
if the effects on channel geometry caused by changes in flow, suspended sedi- 
ment, break-up and other physical factors are known. A third clue to estimation 
and prediction of bank stability lies in the analysis of ages of successional 


vegetation on point bars and slip-off slopes situated cross-channel from 


retreating cut banks. A fourth approach (used in this study) is to examine 
the processes of erosion at specific delta sites and relate them to the 


character, shape and evolution of the channel banks. 


Area and Site Selection 


A study area in the upper delta south of 68°N (Figure 1), was chosen 
for several reasons: 
Leis inenapexcotsthe delta asrthe site of «highest fluvial 
energy and therefore the site of highest banks and 
highest erosion rates. 
2. Extensive field reconnaissance during 1971 and 1972 
showed that a wide range of eroding bank types exist 


on all sizes of channels with a variety of orientations. 


* This report is a summary of a thesis submitted by the author for partial 
fulfillment of an M.Sc. degree at the University of Alberta under the title 
of: Bank Erosion in, the. Southern.Mackenzie. River Delta, NiW. i... Canada’ 


85 pages. 
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3. Most of the channels are deep enough for navigation by 
boats with a 1 m draft through the summer period until 
Eee Zemup. 

4. The study area is easily accessible from Inuvik. 

5. No previous studies of erosion have been done in the 
study area. 

The locations of specific study sites are shown in Figure 2. These 
sites are characteristic of the varying erosional behaviour seen along the 
channels that form the upper delta complex. Several small reversing channels 
and one large one, a section of Peel Channel above and below Indian Village 
(Figure 2), are included as are channel segments with large ice wedges (Plate 1) 
and ones with ice-rich sediments (Plate 2). The levees of the area are 
approximately 10 m above average late summer low water ( Mackay, , 1963). The 
details of the 15 sites chosen for intensive study are as follows: 

Site A - an overhanging bank with a deep thermo-erosional niche 

(hereafter referred to as a niche), located on a NW oriented 
reach; ice layers visible in the thawing bank face; high 
current velocities near the bank levee; surface covered by 
a mature white spruce forest; large stumps with well 
developed adventitious rooting* in the bank material. 

Site B - an overhanging bank with a shallow niche; a former 

point bar deposit; no visible ree in the bank; very 
high current velocities with boils and eddies; a 
reversing flow during break-up. 

Site C = a concave bank similar to Site A. with 2 difterent 

orientation; ice layers visible in the bank. 

Site D - situated on a concave bank of the reversing part of 

Peel River; bank appearance similar to that of Site C. 
Site E - typical of an eroding island with a willow-alder plant 

association; low ice-content sediment; very little 

vegetation mat overhang; current fairly strong; NW 


Orientation of the channel reach. 


* White spruce can send out lateral roots into new additions of alluvium 
enabling it to survive the associated upward movement of the permafrost 


surface occurring in such environments. 
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Site F - very similar to E although on a convex bank; 
erosion is by current coming out of the channel 
running past Indian Village; orientation 
Similar to E but shorter fetch in a NNW Direction. 

Site G - on a convex bank eroded by current from the Peel 
River entering the Peel Channel; bank much lower 
than others (4 m instead of 8 m); erosion of a 
former slip-off slope. 

Site H - on a short concave bank backed by a levee of the 
Middle Channel; this site has been cleared of 
trees by man. 

Site I - similar to site C; a concave bank of the Peel River. 

Site J - a short concave section in a convex bank being eroded 
by current from the Peel River. 

site sk rae concave devee Of the East Channel ;=cnosensaserepre— 
sentative for this particular size of channel. 

Site L - representative of many eroding levees with beaches 
along the Middle Channel; a NNW orientation with long 
open-water fetch and active wave erosion. 

Site M - similar to L except for a larger beach and some ice- 
rafted pebbles; no measurable current velocity at 
sites L and M when chosen. 

Site N - chosen as representative of eroding levees along 
channels the size of Peel Channel; on a concave 
bank, with some vegetation on slopes. 

Site O - similar to site N, although bank material more bonded 


by roots and driftwood. 


Only banks which appeared to be actively eroding were chosen for 
continued observation. Choices were made in August, 1972. Field work 
and air photo interpretation were carried out over the following year. 

Many banks in the area are relatively stable and did not warrant 
observation. Prograding point bars with temporary strandline erosion* 
(Plate 3) were omitted as were highly vegetated banks (Plate 4) and 
banks with sufficient vegetation (Plate 5) to minimize annual erosion. 


*Strandline erosion here refers to fluvial erosion at the water line 
along a channel. A succession of eroded strandlines can develop with 


a drop in water levels. 
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Methods and Data Coilection 


As mentioned previously, one method of determining bank erosion and 
retreat is to examine air photographs. Photos for the area taken in the 
years 1950, 1962, 1965, 1971 and 1972 were examined and estimates of bank 
retreat were made for the 23-year period of aerial photo coverage. 

Data collection in the field involved the establishment of control 
points at all sites in August, 1972, These points were either on trees 
or consisted of stakes placed sufficiently deep in the ground that frost- 
heaving would not significantly affect horizontal control. The distance 
from the control point to the bank edge was measured as in Figure 3. 

Stereo photos (Plate 6) were taken of the banks at the time of 
measurement. This provided a scaled visual record of: 1) the distance 
of the bank face to the horizontal control point; 2) the height of the 
bank above water level; 3) the shape of the bank; and 4) the extent of 
roots and driftwood incorporated in the sediment. Measurements were made 
every 2 weeks in 1973, covering the period from spring ice clearance 
through to the end of August. 

A number of factors considered important in bank erosion were 
measured during the field program. Current velocity measurements (Table 1) 
were made every two weeks 5 m out from the water's edge and 1 m below the 
water surface. Windspeed and direction were monitored continuously with a 
Lambrecht Woelfle type recorder located at the junction of the Middle and 
East channels (Figure 2). The results are presented in the form of a wind 
rose for the period of measurement in 1973 (Figure 4). The open-water fetch, 
measured in a NNW direction to coincide with the prevailing wind, is shown 
in Table 2 and is later used in multivariate analysis of the erosion problem. 
Wind blowing upstream over a long open-water fetch produces larger waves than 
a wind of the same velocity blowing downstream. 

Sediment samples taken at a minimum of 3 points along the bank profile 
were tested for moisture content, liquid limit and, in some cases, for particle 
Size distribution. Averaged results for moisture content and liquid limit are 


recorded in Table 3 and particle size distribution ranges are shown in Figure 5. 
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Erosional Processes 


The visible effects of erosional processes were examined and recorded 
at each bank site. Sloughing occurs at many bank faces. Sloughing is the 
falling away of small blocks (several ne from the vertical face of a bank. 
The size of these blocks, where they accumulated, and the type of surface 
left on the bank face were noted. 

Soil flow occurs at bank faces with ice-rich sediment. Thawing releases 
water from segregated ice layers which saturates surrounding thawed sediment 
layers causing an unstable condition. The movement, appearance, and location 
of flows were noted. 

Wave action undercuts some banks at the waterline and builds beaches of 
loosened material through swash and backwash. The extent and location of 
beaches developed by this process was noted along with the approximate heights 
of waves producing them. 

Depths and locations of thermo-erosional niches were noted although 
penetration of the larger niches (more than 3 m) was only estimated due to 
the hazardous Situation at such- locations.” As the niche develops, the weight 
of the overhanging block of soil increases the tension in the bank along a 
vertical line from the apex of the niche to the top of the bank. Eventually, 
the tension exceeds the tensile strength of the sediment, failure occurs, and 
the block falls into the channel. Notes were made on the size of such blocks, 
the dimensions and location of tension cracks, and the nature of features such 
as ice wedges which provide zones of weakness. 

A number of other erosional factors have been considered. Roots are 
important for the bonding of bank materials as indicated, for example, in 
Plate 7. Surface bank temperatures and water temperatures can be significant 
as frozen silt is harder to erode by fluvial action than unfrozen silt. 
Although break-up current velocities may be high, water temperatures close 
to freezing allow only loose material to be eroded. The action of ice alone 
in the delta channel break-up can sometimes be important, but its erosional 
effects are often attenuated by a near shore ice wreckage barrier. Plate 8 
demonstrates the protection afforded banks by near-shore ice. Mackenzie River 
ice also enters the delta and, with high water conditions, is swept onto beaches, 
banks and levees with apparent minor erosional degradation (plates 9, 10 and 11). 
Another factor, the surface vegetational mat retards thaw of bank materials and 


reduces erosion. 


SO 


Results and Conclusions 


Banks at research sites were categorized according to profiles or 
"shapes"! and placed into 5 classes. Idealized profiles considered repre- 
sentative of each class are shown in Figures 6, 7 and 8. Although the 
procedure is arbitrary for banks which are borderline cases, the groupings 
are realistic on the basis of initial field examination in August 1972. 
Plates 12 to 18 illustrate banks in the various classes. 

Figures 9 to 12 show the amount of 1973 linear erosion from the 
date of one measurement to the next. The initial measure of bank recession 
on June 7th, 1973 indicates the total amount from the date of selection in 
August 1972. These non-cumulative graphs show that much of this erosion for 
some banks (shapes 3, 4 and 5) took place during the 1973 break-up (Table 5). 
Other banks (shapes 1 and 2) eroded mostly during the two weeks after break-up. 

Mean rates of bank retreat determined from inspection of aerial photo- 
graphs are listed in Table 4. The August 1973 position of many lengths of 
bank were plotted on old photos and mean yearly rates determined by measurement 
of differences. One example of air photo use is given in Plate 19. 

The association of a number of physical factors likely to be important 
in bank recession was tested statistically using correlation and regression 
techniques. The dependent variable was linear erosion taken every 2 weeks or 
over the year. The independent variables were current, open water fetch and 
sediment ice content. All variables were tested in various combinations. The 
variables, in general terms, are presented in Table 5 for each site. The results 
of the statistical tests are shown in Tables 6 and 7. It appears from the multiplé 
regression analysis that high near-shore current velocities are associated with 
actively eroding banks. This suggests that delta and river channels in perma- 
frost environments shift as elsewhere with erosion occurring on the outside or 
cut bank side of meanders where currents are highest. Large scale sloughing, 
however, is somewhat greater than in more temperate delta environments due to 
the development of thermo-erosional niches and to structural weaknesses associated 
with ice wedges and other forms of ground ice. Banks along channels with favourab! 
orientation for wave production@ long open-water fetch for upstream winds) will 
have erosion accelerated by wave action. Wave erosion is also significant in 


affecting the appearance of banks to produce beaches or strandline erosion. 
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Bank shapes may evolve from one to another due to the primary 
factor of near-shore current velocities mentioned above. Hydraulic 
geometry changes caused by new bars, scouring beneath ice jams or changes 
in channel discharge may change near-shore current velocities. A higher 
velocity will lower the number of the shape and vice versa. For example, 
deposition of a new bar may be associated with higher near-shore current 
velocity resulting in a change from shape 3 to shape 2 or even shape l. 
There is no evolution between shapes 4 and 5 as their distinction is caused 
by differences in the orientation and open-water fetch of a channel reach 
for upstream winds. 

Bank sediments with high near-surface ice content are subject to slow 
attrition by melting and movement of flowing mud and, in addition, to rapid 
degradation by development of thermo-erosional niches at the waterline with 
subsequent large-scale sloughing. Banks made up of low-ice content sediments 
will be more susceptible to sloughing of materials rather than a flowage and 
are less likely to develop thermo-erosional niches. 

Tables 4 and 5 indicate that, within the study area, banks classed 


according to shape experience the following rates of recession: 


Shape 1 5-30 m/yr 
Shape 2 10-15 m/yr 
Shape 3: 5-10 m/yr 
Shape 4: 0-5 m/yr 


Shape 5: 0-5 m/yr 

Table 5 also indicates the character of erosion associated with each 
of the above shapes. Catastrophic erosion (type 1) is caused by the under- 
mining of a bank by a niche which may cause a block to fall. This may occur 
at both bank shapes 1 and 2 if niches are deep enough. Continuous removal of 
material by current and wave action (type 2) also happens at bank shapes 1 and 
2 which have moderate to high current velocities all summer. Continuous removal 
is why these bank shapes have higher erosion rates than other shapes. A com- 
bination of erosion types 1 and 2 may produce the extremely high rate of 28 m/yr 
found at Site B. Intermittent removal of material (type 3) is caused by variations 
in channel discharge or variations in wind velocities and directions. Type 3 
erosion is associated with bank shape 3. Because of soil flow, banks of shape 
4 or 5 with ice-rich sediment will have continously-retreating bank faces all 


summer (type 4). Banks with low ice content sediment will have very slowly- 
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retreating bank faces all summer. Mass wasting by sloughing is much 
slower than by soil flow. Both types 4 and 5 are characterized by 

the removal of sediment during break-up when unconsolidated slough slope 
material is carried away. 

The application of these rates shapes, and behaviour to part of the 
study area is demonstrated in Figure 13. Mapping of other delta areas could 
prove of benefit to agencies concerned with development of logistic networks, 
community services, oil and gas exploration and other activities. Further 
data collection will enhance the basis for making recommendations relevant 


to construction and development within the Mackenzie Delta. 
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Table 4 EROSION RATES MEASURED FROM AIR PHOTOGRAPHS AT LOCATIONS 
SHAPES DETERMINED BY BOAT 


CHECKED BY HELICOPTER SURVEY. 
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Figure 6 Diagrammnatic Side Views of Bank Shape Category 1, in 3 Stages 


ee ee 


— — —_—__----———- eer 


SHAPE 2 


! face 


SHAPE 2 aoe we: Fe ie 


face 


slough slope 


/ Slough slope 
SrA here 
slough slope 
0) 5 10 15 meters 


Figure 7. Diagrammatic Side Views of Bank Shepe Categories 2 and 3 


= §3285> 


SHAPE 4 levee (Picea) 
vegetation 
mat overhang 


face 


slough slope 


beach 


SHAPE 5 ae 
levee (Picea) 

vegetation , 

mat overhang 


face 


slough slope 


strandline erosion 


1) 5 10 15 meters 
ee) 


Figure 8 Diagrammatic Side Views of Bank Shape Categories 4 and 5 
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Figure 12 Portion: of the Study Area Showing Surveyed Bank Shape Categories 


= 3534 


124° 3G" 
' 


3 aeO 
Vor 
ae 
eee 
\ = \ 
Wain 
ee 
iS ‘ 
Tr Ve 
\ 
\ 
O ‘. 
\2 
AZo 


Plate 1 An ice wedge in a bank with ice-rich sediment. Alluviunm layers 
have been upturned by its formation. 


Plate 2 Taking a sample from a bank with ice-rich sediment, Note the 
wet appearance of the thawing sediment and ice layers. 


isos = 


Plate 3 A point-bar near Reindeer Station eroded by strandline erosion 
(By permission, Don Gill). 


Plate 4 A stable bank of a small channel. 


Plate 5 <A relatively stable bank of the Peel Channel. 
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Plate 6 A stereo pair taken at site G from which bank retreat was 
measured (see Figure 3). All scale markers have 0.5 m 
divisions. 
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Plate 7 


Roots retarding the perturbation of a 
photo). 
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niche roof (center 


Of 


Plate 8 Break - up. Moving ice on the left is separated from bank - 
protecting stationary ice on the rieht by a siear line. 


Plate 9 Massive Mackenzie River ice pushed onto a bank at the junction 
of East and Middle Channels. 
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Plate 10 Willow stems broken by the ice pushed up in Plate 16, 


Plate Il Ridges of alluvium formed by the ice shown in Plate 16. 
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Platel2 A bank with shape 1, at stage 1, within a few days of initial 
niche development (June 8, 1973). 


Plate I3: The same bankas in Plate 7, ae %stage'2, 13edays “ater 
(June: Zhe) se 
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Platel4 The same bank as in Plate 7, at stage 3, 40 days later 
Chilveto. 1970) 


Plate #45 A bank with shape 2 (site F). 
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Plate !6 A bank with shape 3 
(site Hye 


Plate 17 A bank with shape 4 (site L). 
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Plate 18 A bank with shape 5 (site Oy; 


Plate 19 A part of Government 
of Canada air photo 
Al12848-11 (1950, from 
6,100 m altitude). 
Erosion is approxi- 
mately 8m/yr at this 
location (near site E) 
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Abstract 


The geomorphic character of Mackenzie tributary watersheds 


is described by a number of morphometric parameters measured 


from available maps. Statistically determined hydrologic- 
geomorphic relationships based on multiple - regression 
analysis are used to estimate selected hydrologic characteris- 
tics of ungauged basins. Mathematical procedures of factor 
analysis are used in conjunction with regression analysis to 
formulate new prediction relationships for mean 50-year and 
100-year floods. On the basis of results obtained so far the 
east side of the Mackenzie River seems to offer a safer route 


for pipelines than the west side. 
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NOTATION 


Morpho logy 
AVE = Average elevation (ft) 
BF = Bifurcation ratio 
CCM = Constant of channel mainenance (ft° ) 
CLA = Latitude (degrees) 
CLO = Longitude (degrees) 
CP = Basin perimeter (miles) 
CSIN = Sinuosity 
DA = Drainage area Geese ) 
DD = Drainage density Giles ies ) 
Pe— ovrean Prequency 
H = Basin relief (ft) 
HCK = Hypsometric coefficient 
HYI = Hypsometric integral 
LB = Basin length (miles) 
MAE = Maximum elevation (ft) 
MCL = Main channel length (miles) 
MEE = Mean elevation (ft) 
MIE = Minimum elevation (ft) 
PAA = Percent area above average elevation 
PAFI = Percent area at elevation 0-2000 ft 
PAE2 = Percent area at elevation 2000-4000 ft 
PAE3 = Percent area at elevation 4000-6000 ft 
PAE4 = Percent area at elevation 6000 ft and over 
PF = Percent area of basin covered by forest 
PGL = Percent area of basin covered by glacier 
PLSW = Percent area of basin covered by lake and swamp 
RA = Area ratio 
RB = Weighted mean bifurcation ratio 
RCs= Circulatory ratio 
RE = Elongation ratio 
RH = Relief ratio x 100 
RL = Length ratio 
SLOP = Slope factor (ft/mile) 
Hydro logy 
MAF = Mean annual flow (rt? — ee oe vial 5 
MAFF = Mean annual flood with 50-yr. return period (ft 38- miles— 
MAFH = Mean annual flood with 100-yr. return period, (ft> s- miles- ) 
T = Mean flood with 2.33 - yr. return period (ft~ s— miles-) 


SUBS 


Introduction 


Thakur and Lindeijer (1973) have discussed a method of describing the complete 
morphology of Mackenzie tributary watersheds. Morphometric characteristics 
of 54 Mackenzie tributary basins with a total drainage area of approximately 


250,000 square miles have already been published. 


Selected hydrologic characteristics related to mean peak flow conditions 

were estimated by flood frequency analysis (Dalrymple, 1960) from hydrologic 
records of 11 gauged basins. Statistical procedures of multiple regression 
analysis were used to develop relationships between hydrologic characteristics 
(dependent-variables) and morphometric characteristics (independent 

variables) of the basins. Criteria of error estimation and reliability 

were defined. The following equations express the relationships between 


hydrologic and geomorphic parameters: 


MAF = 0.00007 (DA) + 0.00026 (MEE) - 0.9391 (RL) + 2.038 (1) 

T 2 =7. 508 (HCK) = 2.576 Cai) =*0-ie GOhCY yy ta10257 (2) 
MAFF = 0.575 (OLO) + 15.45 (DD) - 9.58 (Kb) —" 55.28 sp), 
MAFH = 0.620 (CLO) + 18.23 (DD) = 9290 .(Rb) = Ob2e0 (4) 


An alternative set of equations was developed for a few basins which did not fit @ 
broad homogeneous geomorphic character of other basins in the greup. Thue 

for watersheds with drainage area (DA) greater than 10,000 square miles, 

mean annual flow (MAF) is estimated from the following equation: 


MAF =. 00,0008 (MEE) — 22966 (hb) =) 0.05) (Hd eye (5) 


For basins with a value of the hypsometric integral (HYI) greater than 0.5, 


the mean annual flood is estimated as: 


T= 0.378 (CLO) ~ DUAG026 tiie = For eee (6) 


For basins with longitude greater than 137.0 or less than 123.0 degrees, the 


mean 50-year and 100-year floods were estimated from the following 


equations: 
MAFE =" O20004) “CDA) —98-640 CRE) 4+713.974) (DD) + 155802 (7) 
MAFH = 0.00045 (DA) = 8.777 (RL) + 16.642 (DD) + 14.888 (8) 


Hydrologic characteristics of ungauged basins were then estimated by 
substituting morphometric variables in the established hydrologic-geo- 


morphic relationships (equations 1 - 8.) 


Morphometric Analysis 

Following the general procedure outlined in our earlier report (Thakur and 
Lindeijer, 1973), morphometric characteristics of nine additional Mackenzie 
tributary watersheds were determined. The results are given in Table l. 


Hypsometric curves of the basins under study are shown in Figures 1 - 5. 


From the results it is apparent that watersheds on the eastern side of 
Mackenzie River in general have smaller drainage areas with lower mean 
elevations and gentler slopes than those on the western side. Consequently, 
the flood potential and erosional hazards to pipelines built on the eastern 


side of the Mackenzie River are of lesser magnitude. 


Hydrologic Characteristics 

Hydrologic characteristics of nine ungauged rivers were determined using 
the morphometric parameters in equations 1 - 8. Results obtained are shown 
in Table 2. Flood values are presented in cubic feet per second per square 
mile. The results compare favourably with those recorded in our earlier 


report. 


=O = 


Eactorinalysi>: 


Factor analysis is a technique used to derive a smaller number of 

independent factors from a larger number of inter-correlated variables. 

From these factors, variables which are close to being independent in 

the original set are identified and subsequently employed in the step- 

wise multiple regression analysis to obtain hydrologic-geomorphic relation- 
ships. Factor analysis is carried out utilizing the Factor Analysis Computer 
program BMD 08M from the set of Biomedical computer programs of the Univer- 


sity of California. 


General Procedure 


The analysis is carried out following the procedure suggested by T.V.A. 
(1966). First, factor analysis with varimax rotation is performed on the 
complete set of 30 morphometric variables as determined for the full sample 
of 63 basins. The number of factors obtained depends upon the percentage 
contribution made by each factor to the cumulative proportion of the total 
variance explained by the analysis. If the percentage contribution by 
inclusion of an additional factor is less than 1% of the total variance 
explained, it is not included in the search for independent morphometric 


variables. 


In the identification of independent variables, only those with factor loadings 
greater than 0.800 are selected. Thus, with screening by means of factor 
analysis, the number of variables obtained is reduced to 18. These 

variables are then used in multiple regression analysis to obtain hydrologic- 
geomorphic relationships for gauged basins. For mean 50-year and 100-year 


floods, the results obtained are shown in Table 3. 


Four variables are used in each equation giving a high correlation coefficient 
and a correspondly low standard error of estimate. 


The final equations are, 


MAPF = 7.27929 (CLO) = 0.00256 (MAE )>=— O-G401U CELEW) 4 164010 -( ap) 
-171.26998 (9) 


MAFH = 1.32268 (CLO) + 0.00283 (MAE) - 0.69055 (PLSW) + 14.01378 (DD) 
-187.29584 (1.0) 


2655 oe 


Results 


The foregoing equations could not be applied to all basins, especially 
those with a centre of gravity located at 126 degrees longitude 

or less and maximum elevation in the basin 5,000 feet or less. Most of 
these basins are located on the eastern side of the Mackenzie River. On 
the other hand, as our original sample of gauged basins from which 
hydrologic-geomorphic relationships were obtained comes exclusively from 
the western side of Mackenzie, the above two relationships probably yield 
better results for basins on the western side of Mackenzie, because of 
higher correlations and lower standard errors than equations based on 


multiple regression analysis alone. 


Values of 50-year and 100-year floods based on equations 9 and 10 
are shown in Table 4. The results are only for basins on the western 
side of Mackenzie. Table 5, shows the results of all basins so far 


studied. These results are based on equations 3, 4, 7 and 8. 


Conclusions 


Using multiple regression procedures, flood characteristics of ungauged 
Mackenzie tributary watersheds were estimated. On an equivalent area 
basis, the flood magnitudes of east side Mackenzie tributaries appear 
to be lower than those on the west side, making the eastern side of the 


Mackenzie a safer route for pipelines. 


Factor analysis was performed on morphometric data to screen inter- 
correlated variables. The resulting variables, when used in stepwise 
multiple regression analysis, yielded very stable statistical relation- 
ships with four significant morphometric parameters. However, the 
applicability of the equations is restricted to tributaries on the western 


Side of Mackenzie. 
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Table 2 


Computed Hydrologic Characteristics*® 
(Multiple-regression analysis) 


AREA 

WATERSHED NO. (MILES©) MAF ae MAFF MAFH 
1. Blackwater at (3) (4) 

MacKenzie (E) 1OHCL S15 / 0.63 4,86 OAT LO 
2. Ochre at MacKenzie 

(E) 1OHC2 es 0.69 3.34 dias s: pees! 

3. Big Snith at 

MacKenzie (E) 10HC3 367 0.48 2.26 eo 8.39 
H Lattle Bear at 

MacKenzie LOKAL COLT Oabz 3.54 ee a OS 
5. Hanna at 

MacKenzie (E) 1OKA2 260 Oued 3.05 14, 14 15.53 
6. Oscar at 

MacKenzie (E) 10KA3 306 0.28 3.86 132.36 14737 
7. Ramparts at (2 

MacKenzie LOKD1 2908 0.50 11.30 ke. O00 Lee 
8. Hume at MacKenzie 10KD2 1254 0.40 (2)44 99 Te SOT 17.43 
9. Donnelly at 

MacKenzie (E) 1OKD3 614 0.64 6.14 163 Tee 


(E) Indicates basins east of MacKenzie River. 
(2,3,4) These values have been estimated fron equation 5,6 and 7 respectively 


* All runoff values are presented in cubic feet per second per square mile. 
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Table 3 


Multiple-Regression Analysis - Significant Morphometric Variables 
Presented in Order of entry into the multiple regression equations 


Variables Estimating Mean 50 —- Year Flood 


MAFF = 10.27182, o = 2.52895, % 0 / MAFF = 24.62 


MAE; CLO PLSW DD 
Correlation coefficient R , . 6697 .7990 .8983 9819 
coefficient of Determination R 4484 . 6383 .8070 9641 
Stepwise increase in R¢ 4484 1899 1687 .1571 
Levelof Significance 5% 5% 1M 1% 
Standard Error of Estimate (0c) 1.9798 1.7 C04 12275 .6188 
7% o/MAFF 19, 27 16.55 12.93 6.02 


Variables Estimating Mean 100 - Year Flood 


MAFH = 11.20455, o = 2.82729, % o/ MAFH = 25.23 


MAE, Ciao PLSW DD 
Correlation coefficient R 3 .6764 ~7994 MOOTS 9850 
Cofficient of Determination R 4575 6390 ioow .9702 
Stepwise increase in R© A575 1815 1491 Misi 
Level of Significance py Of 1% Js 
Standard Error of Estimate (co) 19798 127004 Benya) .618 
% o / MAFH 19.59 16.95 13.88 5.62 
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Table 4 


Computed Hydrologic Characteristics of watersheds on 
the western side of Mackenzie based on Multi-regression 
analysis after screening by Factor Analysis. 


WATERSHED 


Liard at upper crossing 
Frances Near Watson Lake 
Dease at McDame 

Blue at Dease 

Hyland near Lower post 
Turnagain at gauging station 
Kechika at Liard 

Kechika above Boya Creek 
Coal at Liard 

Flat at S. Nahanni 

S. Nahanni above Virginia Falls 
Liard above Meister 
Meister at Liard 
Rancheria at Liard 

Dease at Liard 

Turnagain at Kechika 
Beaver at Liard 

Root at Mackenzie 

N. Nahanni at Mackenzie 
Ram at N. Nahanni 
Tetcela at N. Nahanni 
Keele at Mackenzie 
Dahadinni at Mackenzie 
Johnson at Mackenzie 
Wrigley at Mackenzie 
Redstone at Mackenzie 
Carcajou at Mackenzie 
Arctic Red at Mackenzie 
Ontaratue at Mackenzie 
Wind at Peel 

Blackstone at Peel 
Ogilvie at Peel 

Hart at Peel 

Peel above Canyon Creek 
Snake at Peel 

Bonnet Plume at Peel 
Snake above Iron Creek 
Peel at Mackenzie Delta 
Blow at Arctic Ocean 
Babbage at Arctic Oean 
Firth at Arctic Ocean 
Malcolm at Arctic Ocean 
Mountain at Mackenzie 
Little Bear at Mackenzie 
Ramparts at Mackenzie 
Hume at Mackenzie 


*Values estimated using equations 3 and 4. 
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DA 
13675 
4952 
2719 
682 
3424 
ZG 
8533 
4249 
36 20 
3326 
5653 
3639 
805 
2021 
5754 
2696 
4184 
3835 
Pa he A 
1149 
L503 
10467 
1046 
855 
5O2 
6080 
3527 
8284 
2646 
3783 
1oL2 
2791 
4672 
4709 
2931 
3724 
167%: 
30404 
112 
1756 
2425 
427 
o785 
877 
2908 
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Computed Hydrologic Characteristics of all Watersheds under Study 


Table 5 


(multiple-regression analysis) 


WATERSHED 


Liard above Meister 
Meister at Liard 


. Rancheria at Liard 


Francis near Watson Lake 
Liard at Upper Crossing 
Dease at McDame 

Blue At Dease 

Dease at Liard 

Hyland near Lower Post 
Turnagain at Kechika 


. Turnagain Gauging Station 
. Kechika above Boya Creek 
. Kechika at Liard 

. Coa aleiiamnd 

. Beaver at Liard 


Flat at S. Nahanni 


. S. Nahanni above Virginia 


Falls 


. Root at MacKenzie 
. Willow Lake River at 


MacKenzie (E) 


. River Between Two 


NO. 


LOAAL 
LOAA2 
1LOAA3 
LOABIL 
1OAAOOL 
1OAC2 
10AC4 
1OAC3 
1OAD1 
1OBAL 
LOBAOOL 
1OBA002 
LOBBOO1 
1OBCOO1 
LOBD1L 
10EA003 


10EBOOL 
10GA1 


10GB2 


Mountains at MacKenzie (E)10GB1 


. N. Nahanni at MacKenzie 
. Ram at N. Nahanni 

. Tetcela at N. Nahanni 

. Keele at MacKenzie 

. Dahadinni at MacKenzie 
. Johnson at MacKenzie 

. Wrigley at MacKenzie 

. Redstone at MacKenzie 

. Carcajou at MacKenzie 

~ Avetic Red at Mackenzie 
. Loon at MacKenzie (E) 

. Tieda at MacKenzie (E) 
. Travaillant at MacKenzie 


(E) 


. Ontaratue at MacKenzie 
. Rabbit Hay at MacKenzie 


(E) 


. Pierre at MacKenzie (E) 
. Rengleng at MacKenzie (E) 


Gull at MacKenzie (E) 


. Hare Indian at 


MacKenzie (E) 


. Blue Fish at Hare 


Indian (E) 


. Wind at Peel 

. Blackstone at Peel 

7 Ogilvie at Peel 

Sear sac: Peel 

. Peel above Canyon Creek 


10GD1 
10GD2 
10GD3 
1OHAL 
1LOHBL 
10HB2 
10HB3 
1OHBOOL 
LORBL 
1OLA1 
1OLB1 
10LB2 


10LB3 
10LB4 


10LC1 
10LC2 
10LC3 
10LC4 


10LD1 


10LD2 
1OMAL 
LOMA2 
1LOMA3 
10MA4 
LOMAOOL 


36:3 


(AREA SQ. MILES) 


3639 
805 
2021 
4952 
13075 
2710 
682 
bios 
3424 
2696 
2519 
yalg 
6533 
3620 
4184 
3326 


5653 
3835 


8179 


1305 
aioe 
1149 
1509 
10467 
1046 
855 
502 
6080 
Seri 
8284 
1168 
365 


1124 
2646 


159 
epg 
55 
1662 


1383 


1580 
3783 
1512 
2791 
4672 
9707 


able 5 (Cont'd. ) 


fF) Indicates basins east of MacKenzie River. 


OIRO pee (SP NO PO OS Oy OS) [ney |! (a Ne) ey (OD) 


(miltiple-regression analysis) 


WATERSHEDS 


Snake at Peel 

Bonnet Plume at Peel 
Snake above Iron Creek 
Peel at Mackenzie Delta 
Blow at Arctic Ocean 
Babbage at Arctic Ocean 
Hirth au Arctic Ocean 
Maclecolm at Arctic Ocean 
Mountain at MacKenzie 
Blackwater at MacKenzie (E) 
Ochre at MacKenzie (E) 

Big Smith at MacKenzie (E) 
Little Bear at MacKenzie 
Hanna at MacKenzie (E) 
Oscar at MacKenzie (E) 
Ramparts at MacKenzie 
Hume at MacKenzie 
Donnelly at MacKenzie (E) 


,4 These values have been estimated by alternative equations 7 and 8 respectively. 
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FIGURE 1. Hypsometric curve. 
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FIGURE 2. Hypsometric curve. 
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FIGURE 4. Hypsometric curve. 
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